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A. 


INTRODUCTION 


The  objectives  of  this  program  have  been  to  conduct 
research  on  processing  and  microstructure  relations  on 
materials  suitable  for  infrared  windows  and  to  develop  design 
procedures  appropriate  for  these  brittle  materials.  The 
program  has  been  specifically  concerned  with  producing  high 
strength,  forged  alkali,  or  alkaline  earth  halides,  or  oxides 
and  with  characterizing  the  mechanical  behavior  and  develop- 
ing appropriate  design,  proof  test  or  inspection  techniques 
for  these  or  other  window  materials. 

Specific  areas  of  research  have  included:  a)  investiga- 
tion and  characterization  of  the  substructure  formation 
during  hot  forging;  b)  determination  of  the  effects  of 
additives  or  impurities  on  grain  boundary  mobility;  c) 

' nvestigation  into  the  causes  and  mechanisms  of  recrystal- 
I-  r-iation  or  microstructural  instabilities  including  the 
inhibiting  effects  of  solutes;  d)  evaluation  of  the 
statistics  of  macrostrength  as  influenced  by  statistical 
variations  in  the  strength  of  microstructural  elements;  e) 
develop  flaw  density  curves  for  candidate  materials;  f) 
explore  various  methods  of  non-destructive  evaluation 
suitable  for  finding  small,  sharp  cracks  of  the  type  expected 
in  brittle  materials;  h)  consider  design  procedures 
incorporating  statistical  predictions,  proof  testing,  and 
non-destructive  evaluation  as  appropriate. 

The  work  on  grain  boundary  mobility  focused  primarily 
on  further  elaborating  the  solute  drag  effects  and  refining 
the  theoretical  models  for  solute  drag  and  for  the  mechanism 
by  which  a moving  boundary  can  breakaway  from  a solute 
cloud.  The  solute  drag  model  was  modified  to  include  strain 
interactions  between  aliovalent  solute  ions  and  the  boundary 
in  addition  to  the  primary  electrostatic  forces.  The 
method  by  which  a moving  boundary  can  breakaway  from  its 
solute  cloud  was  considered  as  this  may  be  important  during 
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recrystallization  of  forged  materials  and  also  in  breakaway 
of  boundaries  from  pores  or  precipitates.  This  latter 
process  is  a particularly  important  problem  during  sintering 
of  oxide  materials  to  high  density.  Finally  an  extensive 
analysis  of  grain  growth  data  in  ceramics  was  performed  to 
put  into  perspective  the  importance  of  solute  drag  relative 
to  other  mechanisms  which  can  control  boundary  mobility  in 
oxides . 

Forging  of  CaF2  was  studied  with  particular  emphasis 
on  the  substructure  formation  during  hot  forging.  Crystals 
of  various  orientations  were  forged  over  a wide  temperature 
range  to  find  conditions  under  which  crack  free  forgings 
could  be  accomplished  and  to  investigate  the  subgrain 
formation,  microstructural  uniformity  and  extent  of  re- 
crystallization. Extensive  effort  was  devoted  to  micro- 
structural  examination  using  optical  and  electron  microscopy. 

The  effort  on  designing  with  brittle  materials  focused 
on  the  problems  of  detecting  flaws  and  developing  flaw  density 
curves  and  on  the  effect  of  microcrack  coalescence  on  brittle 
fracture  statistics.  Analytic  and  numerical  models  were 
developed  to  evaluate  the  effects  of  microcrack  coalescence 
on  the  size  and  principal  stress  effects  on  fracture 
statistics.  These  models  also  indicate  the  extent  of  pre- 
cracking to  be  expected  before  failure.  Several  methods 
were  evaluated  for  determining  flaw  distribution  curves 
using  non-destructive  or  destructive  techniques.  These 
results  provide  important  elements  for  designing  with  brittle 
materials  and  for  monitoring  structural  parts  with  acoustical 
emission  to  anticipate  failure. 

The  research  on  this  program  has  been  performed  by  three 
research  groups  and  has  involved  a significant  contribution 
from  students.  Those  who  have  contributed  include  Professors 
H.  K.  Bowen,  R.  M.  Cannon,  F.  A.  McClintock  and  J.  B.  Vander 
Sande,  and  research  assistants  P.  T.  Chiu,  A.  M.  Glaeser, 
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A.  M.  Glae&er  was  partially  supported  by  a 
the  Fanny  and  John  Hertz  Foundation, 


f 


I 


and  M.  F.  Yan. 
fellowship  from 


7 


B.  COUPLED  SPACE  CHARGE  AND  SIZE  MISFIT  CONTRIBUTIONS 

TO  SOLUTE  SEGREGATION. 

B.l.  Introduction: 

Solute  segregation  near  grain  boundaries  has  been 
observed  in  ceramics  as  well  as  metals.  Kingery  (1974)  and 
Gleiter  and  Chalmers  (1972)  have  extensively  reviewed  the 
driving  forces  for  the  segregation  phenomena.  The  major 
factors  leading  to  equilibrium  segregation  in  ceramics  are 
probably  the  electrostatic  potential  associated  with  a grain 
boundary,  and  the  reduction  in  strain  energy  of  solute  ions 
at  and  near  a grain  boundary,  Kingery  (1974).  However,  since 
the  strain  field  affects  the  distribution  of  one  of  the 
cha,.ged  species  and  therefore  the  charge  density  in  the  space 
charge  region  near  the  grain  boundary,  the  electrostatic 
potential  and  the  strain  field  are  coupled.  The  coupling 
between  these  two  interaction  potentials  also  leads  to  coupling 
of  the  differential  segregations  of  the  aliovalent  dopants  of 
the  same  valence  but  with  a different  degree  of  size  misfit 
with  the  matrix  ions. 

In  a moving  grain  boundary,  the  models  of  Cahn  (1962)  and 
Lucke  and  Stuwe  (1962)  for  the  solute  drag  depend  strongly  on 
the  detailed  form  of  the  interaction  forces.  Thus  such 
coupling  can  be  expected  to  affect  the  solute  drag. 

This  section  treats  the  equilibrium  interaction  forces  and 
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solute  distributions.  Applications  to  the  case  of  a moving 
boundary  will  follow  in  section  C. 

B.2.  Theory: 

The  existence  of  space  charge  regions  near  lattice 
discontinuities  (free  surfaces,  dislocations  and  grain 
boundaries)  in  an  ionic  solid  was  first  postulated  by  Frenkel 
(1946) ; the  charge  distributions  have  been  formulated  by 
Frenkel  (1946),  Lehovec  (1953),  Eshelby,  et  al.  (1958)  and 
Kliewer  and  Koehler  (1965) . The  potential  difference  between 
the  surface  and  the  bulk  leads  to  a nonuniform  distribution 
of  axiovalent  solute  ions  in  a region  which  typically  extends 
2-20  lattice  spacings  in  from  the  surface. 

Gibbs  showed  that  solute  elements  which  decrease  the 
surface  energy  should  concentrate  at  the  surface.  Many  other 
thermodynamic  approaches  after  Gibbs  have  been  reported.  In 
particular,  McLean  (1951)  showed  that  for  an  impurity  with 
an  adsorption  free  energy,  F^,  and  a bulk  concentration,  C^, 
the  solute  concentration  at  the  boundary,  Cj^,  is  given  as 

C exp  (-F  AT)‘ 

C = (1) 

1 + exp  (-^) 

For  a solute  ion  which  has  a size  misfit  in  the  n.atrix 
the  adsorption  free  energy  can  arise  from  the  partial  reduction 
of  the  strain  energy  for  an  atom  segregated  to  the  already 
distorted  boundary  regions.  Furthermore,  due  to  the  spatial 
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dependence  of  the  " attice  distortion  in  the  grain  boundary 
region,  the  adsorption  free  energy  is  likely  to  have  a 
certain  spatial  dependence.  However,  the  details  of  the 
spatial  dependence  are  related  to  the  grain  boundary  structure 
andare  not  known  with  precision.  The  range  of  the  strain  field 
normal  to  the  boundary  likely  corresponds  to  the  periodicity 
of  the  structural  units  in  the . coincidence  model  of  the 
grain  boundary  structure  and  is  about  several  lattice  con- 
stants, Chalmers  and  Gleiter  (1971).  The  maximum  magnitude 
of  the  strain  energy,  U^,  can  be  approximately  estimated  from 
the  size  misfit,  Ar,  of  the  solute  ion.  Young’s  modulus,  Y, 
and  Poisson's  ratio,  v,  as  given  by  Lucke  and  Stuwe  (1962)  as 


U 


o 


47fY  r-^  n 
(1  + V) 


(2) 


A IT 

where  n ^ - and  r is  the  radius  of  the  lattice  ions  of  the 

sublattice  in  which  the  solute  ion  is  substitutionally 
located . 

In  an  ionic  solid,  the  strain  field  interaction  between 
a grain  boundary  and  solute  ions  affects  the  charge  distribu- 
tions and  therefore  the  electrostatic  potential  distribution 
in  the  grain  boundary  region.  For  illustrative  purposes,  we 

analyze  the  impurity  segregation  of  divalent  cations,  I, 

2+  2+  2+ 

e.g. , Sr  , Ca  , Ba  in  the  grain  boundary  region  of  an 
ionic  solid,  e.g.,  KCl. 

We  use  the  approach  of  minimizing  the  free  energy  of  the 
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crystal  including  the  lattice  defects  and  solute  ions  and 
the  interaction  between  the  charged  species  and  the  grain 
boundary,  following  Kliewer  and  Koehler  (1965).  In 
formulating  the  free  energy  of  the  system  of  defects  within 
a distance,  L,  from  a grain  boundary,  we  include  the  strain 
field  and  the  electrostatic  interaction  energies  between 
the  defects  and  the  grain  boundary.  For  generality,  we 
assume  strain  interaction  energies,  Ug(x),  between  the  defects 
(vacancies,  vacancy  pairs,  vacancy-solute  pairs,  and  unassociated 
solute  ions)  and  the  grain  boundary.  In  particular,  we  assume 
that  the  strain  interaction  energy,  Ug(x),  between  the  un- 
associated divalent  cation  dopant,  I,  and  the  grain  boundary 
has  the  functional  form. 


Uq  (1  - (^)")eV  X < 2a 


U(x)^  = 


(3) 


0 eV 


X > 2a 


where  a is  the  lattice  constant  of  the  matrix,  and  a choice 
of  the  index,  n,  affects  the  detail  of  the  strain  function. 
The  range  of  strain  field  energy  is  chosen  as  two  lattice 
constants  and  is  within  the  periodicity  of  the  structural 
units  in  the  coincidence  grain  boundary  model,  Chalmers  and 
Gleiter  (1971).  For  specific  calculations  the  magnitudes  of 
Uq  are  chosen  as  + 0.2  eV  and  +0.05  eV.  A strain  energy, 
Uq/  of  -0.2  eV  is  probably  appropriate  for  Ca  and  -0.05 
eV  for  Sr^^  in  KCl. 


11 


lurtherraorc,  we  assume  ^-hat  a certain  binding  energy 

3S  rcyaired  for  an  excess  quantity  of  the  divalent  cation, 

I,  a;id  the  m,  trix  ions,  K and  Cl,  to  be  associated  within 

the  grain  boundary  core.  This  binding  energy,  similar  to 

the  delta  function  discussed  in  at  ai..  (1965)  can 

also  be  c function  oi  the  particular  grain  boundary. 

Tee  free  energy,  F,  j.t.  a system  with  number  densities 

of  n_^(x)  V3cm„a,3s , n_{x)  anion  vacancies,  n (x) 

B 

vacancy  corap  . ^ xes,  nj,(x)  unassociated  solutes  and  n^j^(x) 
associated  vacancy-solute  complexes,  can  be  expressed  as 


F = 


1, 

A 

dxtn  (?  tU^)  + n (F  + U~)  + 
j ^ s - s 

o 


(F  + F - B + u;  ) + n,U 


s f^s  ^ib 


-n  (F'  -B,  +U^)  +ip^] 


"'^^C  ' 1 n'^f'^  + n f 

I s s s -S  s 


s s 


M) 


where  F and  F are  the  formation  energies  of  the  cation  and 
anion  vacancies  respectively;  3 and  are  the  binding 
energies  of  the  vacancy  complexes,  (V^-V^^)*,  and  the 
vacancy-impurity  complexes,  respectively;  ±p>p  is 

the  electrostatic  energy  term  for  the  charge  density,  p,  with 
the  space  charge  potential,  <t>;  is  the  configurational 
entropy  of  the  system  but  excluding  that  in  the  boundary  core; 

IS  the  electrostatic  energy  term  for  the  boundary  core 
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with  a "surface"  charge  density,  a,  and  a "surface"  potential, 

‘^1'  ^s'  ^s  "umber  per  unit  area  of  the  un- 

associated  I,  K and  Cl  ions  respectively  in  the  boundary  core 

with  their  respective  binding  energies  F^,  f'*'  and  F~  • and  u'*’ 

- B f b s s s'  s' 

^s'  ^s'  ^s'  ^s  ^ strain  field  interaction  energies 

between  the  grain  boundary  and  the  unassociated  cation 

vacancies,  unassociated  anion  vacancies,  cation-anion  vacancy 

pairs,  unassociated  solute  ions  and  solute-vacancy  complexes. 

The  charged  species  in  the  crystal  are  related  by 

Poisson's  equation. 


47re 

e 


(Hf  + n_  - n^) 


(5) 


where  e is  the  magnitude  of  an  electronic  charge,  and  e the 
permittivity. 

The  free  energy,  F,  must  be  minimized  while  subjected 
to  the  constraints  of  Poisson's  equation  and  the  conditions 
of  electrical  neutrality  requiring  the  charge  per  unit  area 
on  the  boundary,  a,  to  be 

(nf  + n_  - n_|^)  dx  (6) 

and  that  the  unassociated  I ions  in  the  boundary  core  and  in 
the  boundary  region  are  conserved, 

L 

dx  = constant  (7) 

o 


+ 

s 
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and  that  the  excess  matrix  ions  in  the  core  are  compensated 
by  the  vacancies  in  the  crystal, 

L 

n’  - n”  = r [n^(x)-n  (x) ] dx  (8) 

s s \ + 

•*o 

These  constraints  are  somewhat  more  restrictive  than 
necessary,  bat  allow  solutions  which  contain  the  important 
elements  of  che  problem. 

After  the  minimization  of  the  free  energy,  F,  with 
these  constranits  we  obtain 
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where  N is  the  derisity  of  cation  sites  and  is  the  Langrange 
multiplier  which  is  determined  from  the  condition  that  the 
total  number  of  impurity  ions  must  be  conserved,  is  the 
total  impurity  concentration  in  the  bulk  and  p^  is  the 
fraction  of  these  associated. 

The  bulk  potential,  e(|)oo,  can  be  determined  from  the 
electrical  neutrality  condition  in  the  bulk, 


n^  (“)  = n_  (~)  + n^  (“)  , 


exp 
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With  the  assumption  that  and  have  no  spatial 
dependence,  the  Poisson's  equation  can  be  reduced  to 
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and  the  approximation  holds  in  the  extrinsic  temperatures 
where,  n_  («>) /n_j_  (oo)  <<i. 

Numerical  integrations  were  performed  to  solve  Poisson's 
equation,  Eq.  (14),  for  the  strain  field  of  the  functional 
form,  Eq.  (3)  and  for  the  total  impurity  concentration  of 
50  ppm  molar.  The  values  of  defect  formation  and  association 
energy  used  for  these  calculations  are  listed  in  Table  I. 

At  temperatures  above  the  isoelectric  point,  the  I ions  are 
electrostatically  repelled  from  the  postively  charged  grain  boundary 
An  attractive  strain  field,  < o,  counter-balances  the 
electrostatic  space  charge  effect  and  Increases  the  space 
charge  potential  in  the  charge  cioud.  However,  a repulsive 

field,  > 0,  decreases  the  space  charge  potential  to 
a more  negative  value.  Figs.  B1  and  B2. 

At  temperatures  below  the  isoelectric  point,  which  is 
about  444"c  for  this  impurity  concentration  level,  the  attract- 
ive space  charge  potential  is  complemented  by  an  attractive 
strain  field  such  that  the  potential  near  the  boundary  may 
exceed  the  bulk  potential.  Figs.  B3  and  B4.  As  seen  from  Eq.  (13) 
and  the  Figs.  B1  - b4  the  value  of  e*„  is  not  affected  by 
the  strain  interactions;  this  is  to  be  expected  physically. 

The  distribution  profiles  of  the  K vacancies  and  the  I 


TABLE  I 


Data  of  KCl  Vacancy  Formation  Energies,  and  Solute-Vacancy 
Binding  Energies. 


Quantity 

K— vacancy  formation  energy^ 
Ql~vacancy  formation  energy^ 
Schottky  defect  formation  energy^ 
I-solute  K vacancy  binding  energy 


Enthalpy 

0.84 

1.38 

2.49 


(eV) 


Entropy  (k) 
3.2 
3.2 
7.64 
1.88 


1.  Dreyfus  and  Nowick  (1962). 

2.  Fuller,  et  al.  (1968) . 
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ig.  Potential  distributions  in  the  grain  boundary  region  due  to  the  coupled 

effects  of  the  electrostatic  space  charge  and  the  mechanical  strain  field 

interactions  at  400°C  (<  T ) . 


ions  in  the  grain  boundary  region,  for  an  attractive  impurity 
strain  field,  u = -o.2eV,  are  shown  in  Fig.  b5.  At  the  location 
where  the  concentration  of  K vacancies  and  I ions  are  equal,  there 
is  no  net  charge  density  (with  Cl  vacancies  neglected);  and  this 
location  corresponds  to  the  inflection  point,  in  the 

potential  distributions  shown  in  Figs. B1-B4 . The  segregation  of  I ions 
within  two  lattice  spacings  of  the  grain  boundary  is  due  in  large 
part  to  the  attractive  strain  field.  Beyond  the  range  of  the 

strain  field  interaction,  the  space  charge  potential  is  repulsive 
to  I ions  even  at  400°C,  Fig.  b5. 

Several  functional  forms  of  the  strain  field  interaction  energy 

are  assigned  by  the  designation  of  n = 2,  1,  and  1/4  in  eq.  (3).  Due  to 

the  coupling  with  the  strain  field,  the  electrostatic  space  charge 
potential  distributions  are  different  for  different  strain  functions. 
The  distributions  of  the  charged  species  due  to  the  coupled  effects  of 
the  electrostatic  space  charge  potential  and  the  different  strain 
functions  are  shown  in  Fig.B6  for  a temperature  just  above  the  iso- 
electric temperature,  444°C. 

For  the  case  of  a sample  with  two  divalent  cation  solutes 
of  the  same  ionic  charge  but  with  different  ionic  radii,  Poisson's 
equation,  Eq . (14),  in  the  grain  boundary  regions  becomes 


d^2  ^ 1 

ds2  " 2 


(18) 


where  denotes  the  strain  energy  function  of  solute 
of  solute  2;  and  A is  the  fraction  of  the  total 
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Fig.  B6.  Calculated  values  of  the  solute  distribution 
near  the  grain  boundary  for  three  different 
assumed  forms  of  the  strain  field  interaction 
energy. 


divalent  cation  which  is  solute  2.  In  the  limit  of  a 

negligible  amount  of  dopant  2,  the  Poisson's  equation  is 
reduced  to  Eq.  (14) . 

For  the  case  where  the  ion  with  the  greater  misfit 
IS  the  smaller  fraction  of  total  impurity,  Fig.  B7  shows  the 
concentrations  of  the  impurities  calculated  from  Eq.  (18) 
and  the  expressions  for  the  unassociated  solute  densities, 
'^fl  ^f2  solutes  1 and  2 respectively,  are 


Hf^ (x) 
N 


= (1  - A)C  (1 ~ p )e 

00  ' 00 


V 


f j ^ 


(x)-U^ (») 
"kr 


) 


(19) 


N 


- AQ„  (1  - p }e 


/ I 

J 


kT 


(x)-U2 (") 


) 


(20) 


where  the  impurity-vacancy  binding  energy  and  the  probability 

of  association  are  assumed  to  be  identical  for  solutes  1 and 

2. 


B.3.  Discussion: 

(1)  A coupling  between  the  electrostatic  potential  and 
the  mechanical  strain  interaction  between  an  aliovalent  solute 
and  the  grain  boundary  is  shown  to  exist  and  to  induce  a non- 
uniform  distribution  of  the  solutes  and  the  electrostatic 
potential  in  the  grain  boundary  region. 

(2)  The  coupling  is  visibly  strong  at  temperatures 
close  to  the  isoelectric  temperature  where  the  electrostatic 
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potential  is  small.  Under  these  conditions  the  electric  fields 
will  be  much  stronger  than  in  the  strain  free  case,  and  both  the 
potentials  and  fields  may  reverse  their  sign  within  the  charge  cloud 

(3)  The  choice  of  the  functional  forms  of  the  strain 
field  potential  is  rather  arbitrary.  However,  the  magnitudes 
of  are  chosen  in  that  they  are  reasonable  values  for  Sr 
and  Ca  dopants  in  a KCl  matrix.  The  range  of  the  strain 
field  are  chosen  as  two  lattice  spacings,  and  this  corresponds 
to  the  periodicity  of  the  structural  units  in  the  coincidence 
model  of  the  grain  boundary  structure. 

(4)  The  solute  profile  in  the  grain  boundary  region  is 
more  sensitive  to  the  functional  form  of  the  strain  field  at 
temperatures  close  to  the  isoelectric  point. 

(5)  If  the  strain  interaction  goes  to  zero  at  infinity, 
the  bulk  potential,  e<^^,  and  the  isoelectric  temperature  are 
unaffected  by  the  strain  interaction  and  are  the  same  as  for 
the  purely  electrostatic  case. 

(6)  Although  the  binding  energies  for  the  defects  in 
the  boundary  core  are  included  in  the  initial  formation  they 
do  not  affect  the  results  in  an  important  way  without  making 
additional  restrictive  assumptions  about  the  surface  or 
boundaries.  As  discussed  by  Lifshitz  et  al,  the  terms 
must  be  functions  of  the  particular  boundary;  in  tie  absence 
of  such  specific  information,  if  the  terms  are  taken  as 
constants  they  have  no  interesting  effect  and  the  quantities 
(f"  + F^)  become  the  defect  creation  energies  used  by  others. 

(7)  Similarly  the  binding  energy  of  cations  to  the  core 
does  not  affect  the  spatial  distribution  of  potential,  cations 


27 


or  vacancies  near  the  boundary  except  indirectly  through  the 
affect  on  the  bulk  concentration,  C^.  This  effect  on  is 
negligible  if  L is  large.  At  high  concentrations  and  low 
temperatures  when  the  impurity  concentration  at  the  boundary 
core  becomes  high  enough  to  saturate  the  available  sites,  then 
a more  explicit  effect  would  be  anticipated  which  would  be 
evident  if  the  surface  were  treated  as  having  a finite  number 
of  sites  in  the  analysis. 

(8)  The  analysis  of  the  free  energy  minimization  im- 
plicity  assumes  a continuous  charge  distribution.  Although 
this  is  probably  not  correct  as  pointed  out  in  Whitworth, 

(1972,  1975)  the  resulting  errors  do  not  obscure  the  present 
results  regarding  the  effects  of  strain  interactions  on 
segregation. 
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C.  SOLUTE  DRAG  THEORY  FOR  GRAIN  BOUNDARY  MIGRATION  IN  KCl . 

C.l.  Introduction: 

Alkali  halides  are  appreciably  strengthened  by  the 
deformation  structure  which  forms  during  hot  forging. 

These  microstructures  are  comprised  of  low  angle  subgrains, 
typically  of  a size  of  lOym  or  less,  eg..  Cannon  et  al . 

(1976),  Yan  et  al.  (1975a, c).  In  pure  materials  these 
forged  crystals  were  found  to  be  susceptible  to  strength 
degradation  caused  by  recrystallization  and  growth  of  grains 
of  the  order  of  1 mm  in  size  during  stress  relief  anneals 
and  even  at  room  temperature,  eg.,  Koepke  et  al . (1974). 

The  additions  of  solutes,  particularly  divalent  cations, 
were  found  to  be  effective  in  inhibiting  this  undesirable 
recrystallization.  Although  the  effects  C-  solutes  on 
nucleation  are  still  not  understood,  the  inhibition  of 
recrystallization  is  apparently  due,  at  least  in  part,  to 
the  reduction  in  grain  boundary  mobility  caused  by  the 
additives,  Yan  et  al.  (1975a, b,c).  Cannon  et  al.  (1976). 

In  the  impurity  drag  theories  of  Cahn  (1962) , Lucke  and 
Stuwe  (1962) , the  nature  of  the  interactions  between  the 
impurity  atoms  and  the  grain  boundary  is  an  important 
feature.  For  grain  boundary  migration  in  metals,  the  inter- 
action energy  has  generally  been  assumed  to  be  a mechanical 
or  strain  field  interaction.  However,  the  experimental 
results  in  ionic  solids,  showed  that  aliovalent  cation 
solutes  have  a much  stronger  effect  than  monovalent  solutes 
in  retarding  the  grain  boundary  migration,  Yan  et  al.  (1975c) . Yan 
(1976)  . Further , the  space  charge  theories  showed  that  ionic 
crystal  interfaces  at  thermal  equilibrium  with  the  bulk 
crystal  are  electrostatically  charged. 

Consequently,  we  proposed  that  the  impurity-boundary 
interaction  energy  is  mainly  due  to  the  electrostatic 
interactions  between  the  space  charge  cloud  and  the  moving 
grain  boundary.  A theory  was  developed  on  the  previous 
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program  for  the  solute  drag  caused  by  divalent  cations  in 

KCl;  It  was  based  exclusively  on  electrostatic  interactions. 

he  theory  was  an  extension  of  the  equilibrium  space  charge 

heory,  and  allowed  calculations  of  the  mobility  based 

upon  experimentally  determined  values  of  solute  diffusivities 

and  other  eefect  energies  and  mobilities,  Van  et  al.  (1975c) 

The  only  parameter  in  this  idealised  theory  which  is  not 

sufficiently  well  known  is  the  difference  in  energy  for  a 

cation  vacancy  located  at  a boundary  compared  to  one  in  the 
lattice . 

In  the  case  of  a monovalent  dopant,  the  strain  field 
interaction,  rather  than  the  space  charge  potential,  is 
more  likeiy  to  be  the  dominant  interaction  potential  between 
he  dopant  and  the  grain  boundary.  It  is  generally  agreed 
that  the  strain  field  interaction  potential  originates  from 
a site  mismatch  between  the  dopant  and  the  matrix  ions, 

Lucke  and  Stuwe  (1962)  . However,  for  an  aliovalent  dopant 
wi  a significant  ionic  size  mismatch  with  the  matrix  ion 
in  the  same  sublattice,  both  the  strain  field  and  the 
electrostatic  interactions  are  significant.  The  strain 
field  interaction  energy  affects  the  distribution  of  the 
charged  solute  ions  and  therefore,  the  electrostatic 
potential  and  defect  distribution  in  the  grain  boundary 
region.  In  section  B,  we  analyse  the  coupled  effects  of 
the  strain  field  and  the  electrostatic  interactions  on  th, 
solute  and  space  charge  distribution. 

This  section  includes  two  extensions  of  the  grain 
boundary  mobility  model.  First  we  incorporate  the  strain 
leld  interaction  into  the  solute  drag  theory.  This  is  of 
particular  importance  at  temperatures  near  the  Isoelectric 
temperature  where  the  electrostatic  forces  become  low  enough 
that  the  solute  cloud  may  not  remain  attached  to  the  boundary 
in  the  absence  of  a strain  field  interaction.  Secondly  we 
consider  the  mechanism  by  which  a moving  boundary  can 
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breakaway  from  the  solute  cloud  and  move  at  high  velocity. 
It  had  been  previously  suggested  that  this  is  an  important 
step  in  the  discontinuous  recrystallization  observed  at 
low  temperatures  in  KCl,  Yan  et  al.  (1975a),  Cannon  et  al. 
(1976)  . 

C.2.  Calculation  of  Solute  Drag: 

The  solute  drag  force  can  be  calculated  from  the 
asymmetric  solute  distribution  about  the  boundary  at  a 
given  velocity.  The  solute  profiles  are  determined  by  sol- 
ution of  the  flux  equation;  this  requires  explicit  expres- 
sion of  the  vacancy  concentrations  near  the  boundary  in 
order  to  calculate  the  dif f usiviry . Considering  only 
electrostatic  effects  the  solute  diffusivity  :an  be 
expressed  as,  Yan  et  al . (1975b),  Yan  (1976). 

D^(x)  = D^p(x)  = (1) 

p®+  (l-Po,)exp(-^i^j,-~  ^‘^°?) 

where  is  the  diffusion  coefficient  of  the  associated 
vacancy-solute  complexes,  and  p(x)  is  the  probability  of 
association  near  che  boundary 


p«,  + (l-poo)exp(-m.^j,~  ^^fi) 

The  electrostatic  field,  has  the  value  in  the  bulk. 

The  concentration  of  vacancy-solute  complexes  in  the  bulk  is 
Coopoo,  where  C<»  is  the  nominal  solute  concentration. 

Such  calculations  were  preformed  previously;  however, 
it  is  computationally  simpler  to  evaluate  the  solute  drag 
from  an  approximate  expression,  which  is  constructed  from 
Taylor's  expansions  of  the  integral  form  of  the  drag  force 
at  the  high  or  low  velocity  limits,  Cahn  (1962) . 
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consequently,  the  total  drag  force,  F,  experienced  by  a 

moving  grain  boundary  with  velocity,  V,  can  be  approximated 
as 


F 


g Cco  V 

2 2 
1+B  V 


(3) 


where  Mq  is  the  intrinsic  boundary  mobility,  is  the 

mobility  per  unit  impurity  concentration  in  the  solute 
controlled  low  velocity  region;  and  6“^  is  the  drift 
velocity  of  the  impurity  atoms  across  the  grain  boundary 
and  it  is  also  the  boundary  velocity  at  which  the  impurity 
drag  reaches  the  maximum.  Furthermore,  as  shown  in  Cahn 
(1962),  the  quantities  a and  6 are  related  to  the  impurity- 
boundary interaction  energy,  U(x),  and  the  impurity 
diffusivity,  D(x),  as 

r . 2 

g = 4NkT  sinh  U(x)/2kT 


-2-  ^ ' ,du.  ^ 

g2  kT  j 

^ J 


D(x)  dx 


(5) 


where  N is  the  number  of  cations  per  unit  volume,  and  we  take 


U(x)  = e<p  (x)  - e(po 


(6) 


This  procedure  can  be  directly  modified  to  include 
the  strain  interaction  energy  by  using 


U(x)  = (e(p  - e(poo)  + U 


strain 


D(x)  = 


D p 
0^00 


e(p  - e(j)oo  + U ^ 
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''^hsiTG  thiG  GlGctros tatic  potontial  is  modifisd  Hy  t".hr  strain 
interaction  as  discussed  in  section  B. 

Calculations  of  a and  3 ^ have  been  done  for  strain 
functions  of  the  form. 


r 


u = 

s 


(l-(x/2a)'^) 


0 


X < 2a 


X > 2a 


(9) 


where  a is  the  lattice  parameter,  n determines  the  rate  of 
decrease  of  the  strain  interaction  energy  in  moving  away 
from  the  boundary  a distance  x.  The  maximum  interaction 
force  Uq  is  calculated  from  the  elastic  strain  energy  for  a 
sphere  with  a radius  difference  Ar  compared  to  the  host  ion; 

3 

4rrEr  2 

^o  ~'(l+vT  ^r^^ 


where  E is  the  elastic  modulus,  and  v is  Poisson’s  ratio. 
Using  typical  values  for  KCl  with  a Ca"^^,  Sr^^,  or  Ba^^ 
dopant  gives  values  ot  of  -0.2,  -0.05,  and  -10"^eV 
respectively. 

Plots  of  a and  3 ^ are  shown  in  Figs.  (Cl)  and  (C2) 
These  were  calculated  for  n = 1,  and  = +0.2,  and  +0.05. 
The  values  for  Uq  =+  0.2  and  + 0.05  are  less  likely  to  be 
physical  (i.e.,  a strain  energy  interaction  will  always  be 
attractive,  < 0)  but  were  included  for  completeness.  At 
temperatures  near  the  isoelectric  point  the  drag  parameter, 
a / does  not  vanish  for  a finite  strain  field  of  either  sign. 
Fig. (Cl) . The  strain  field  interaction  tends  to  decrease 
the  breakaway  velocities,  3 r over  most  of  the  temperature 
range. 

the  drag  forces  must  be  balanced  by  the  driving 
force  in  order  to  maintain  a steady  state  boundary  motion, 

• (3)  , presents  a cubic  equation  of  the  boundary  velocity 

for  a given  driving  force.  The  solutions  of  boundary 
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Figured.  Calculated  values  of  a,  the  drag  force  per  unit 
velocity,  per  unit  dopant  concentration  on  a 
moving  grain  boundary,  considering  the  coupled 
interactions  of  the  electrostatic  and  mechanical 
strain  fields. 
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Figure  C2.  Calculated  values  of  the  b.reakav\7ay  velocity 

of  a grain  boundary  from  the  associated  solute 
cloud.  Both  the  electrostatic  and  mechanical 
strain  interactions  are  included. 
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5 2 

velocities  under  a driving  force  of  10  dynes/cm  , which 
corresponds  to  an  average  of  7.3ym  for  the  linear  grain  size 
intercepts,  are  shown  in  Figs.  (C3)  and  (C4) . Mobility 
calculations  alternatively  assuming  n = 1 and  n = 2 are 
shown  in  Fig.(C3);  the  purely  electrostatic  case  is  shown 
for  comparison  in  Fig. (C4).  The  values  of  defect  dif fusivities 
and  association  energies,  and  of  Mq  were  those  previously 
used,  Yan  et  al.  (1975b),  Yan  (1976).  The  experimental  data 
for  similar  Sr  concentration,  50  ppm,  and  driving  force  are 
also  shown. 

It  can  be  seen  that  away  from  the  isoelectric  tempera- 
ture, the  effect  of  the  coupled  strain  field  is  relatively 
small  , but  near  the  isoelectric  temperature  the  calculated 
mobility  is  substantially  modified.  For  the  assumed  value 
of  n = 1,  the  general  shape  of  the  curves  is  similar.  There 
is  a region  in  which  the  driving  force  exceeds  the  drag 
force  so  that  the  boundary  would  move  at  the  intrinsic 
velocity  in  this  temperature  range.  For  the  purely  electro- 
static interaction,  the  temperature  range  for  this  high 
velocity  behavior  is  centered  about  the  isoelectric.  For 
the  coupled  case  with  n = 1,  the  temperature  range  for  a 
single  high  velocity  solution  is  smaller  and  further  it  is 
displaced  to  a higher  temperature  range  rather  than  being 
centered  around  the  isoelectric  temperature.  If  a 
repulsive  force,  > 0,  is  used  the  gap  in  the  low  velocity 
drag  curves  shifts  to  lower  temperatures  and  is  centered 
below  the  isoelectric.  For  the  n = 2 assumption  the  detailed 
form  of  the  velocity  curve  is  changed;  the  cause  for  this  is 
not  presently  understood.  However  for  this  case  the  maximum 
in  the  velocity  is  also  shifted  to  a temperature  higher  than 
the  isoelectric.  Where  a low  velocity  solution  does  exist 
it  is  very  similar  for  both  assumptions. 

2 

The  formulation  of  the  expression  for  a and  a/3  , Eqs . 

(4)  and  (5)  were  based  on  the  assumption  that  the  solute- 
boundary interaction  energy  U is  independent  of  velocity. 
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Figure  C4.  Calculated  grain  boundary  mobility 
KCl  considering  only  electrostatic 

for  Sr  doped 
interactions 

This  may  be  so  for  a strain  interaction,  particularly  i i 
dilute  solution,  but  it  is  not  obviously  so  for  the  electro- 
static interactions.  The  origin  of  the  space  charge  cloud 
is  essentially  due  to  the  difference  in  the  distributions 
of  the  charged  species,  namely  the  unassociated  Sr  ions,  and 
K and  Cl  vacancies.  The  space  charge  potential,  cj) , must 
satisfy  Poisson's  equation,  ^ - ilP  . However,  p is 

composed  of  two  major  components,  namely  free  Sr  ions,  Sr^, 
and  free  K vacancies,  The  contribution  of  free  Cl 

vacancies,  , is  not  significant  in  Sr-doped  KCl  in  the 
extrinsic  temperature  region.  Furthermore,  the  relative 
proportion  of  Sr^  and  in  the  total  charge  density,  p, 
also  varies  with  temperature.  ^At  temperatures  above  the 
isoelectric  temperature,  the  constitutes  the  bulk  of  the 
charge  density  whereas  at  temperatures  below  the  isoelectric 
temperature,  the  charge  density  is  mainly  due  to  the  Sr' 
impurity  cloud.  Furthermore,  the  charge  distributions  are 
affected  by  the  strain  field  interaction  between  the  charged 
species  and  the  gr^.in  boundary. 

When  the  grain  boundary  migrates,  the  space  charge 
potential  will  be  affected  by  the  variation  of  the  impurity 
distribution  which  is  a function  of  the  boundary  velocity. 
This  will  introduce  a velocity  dependence  to  the  interaction 
energy,  U,  between  the  dopant  and  a moving  boundary.  We 
previously  showed  that  for  the  purely  electrostatic  case 
the  velocity  dependence  of  was  relatively  small  because 
the  distribution  of  the  more  mobile  vacancies  readjusts  in 
response  to  the  changes  in  solute  distribution;  Yan  et  al. 
(1975b)  . Consequently  the  static  values  of  <p  could  be  used 

in  Eqs.  (4)  and  (5)  or  in  other  calculations  of  solute  or 
defect  distributions. 

The  calculations  for  the  coupled  strain  and  electro- 
static forces  are  more  approximate  than  those  for  the  purely 
electrostatic  case  because  the  detailed  form  of  the  strain 
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field  is  rivTt  known,  but  must  be  assumed  and  treated 
parametrically.  In  addition,  the  assumption  of  a negligible 
velocity  dependence  of  the  electrostatic  field  is  not  as  good 
for  the  coupled  case  as  for  the  purely  electrostatic  case, 
especially  at  high  velocities.  Nevertheless,  the  calcula- 
tions, Fig.c3  indicate  that  the  low  velocity  drag  is 
relatively  insensitive  to  the  detailed  form  of  the  strain 
field  assumed. 

Especially  at  temperatures  near  the  isoelectric  temp- 
erature, the  existence  of  the  strain  field  interaction 
causes  a significant  change  in  the  electrostatic  field  near 
the  boundary,  (see  B,  Figs.  Bl  - B4).  Thus  as  the 

velocity  is  increased  and  the  excess  solute  concentration 
near  the  boundary  decreases  the  electrostatic  potential  will 
approach  that  for  a uniform  distribution  of  solute.  This 
is  much  more  like  that  for  the  static  case  with  Uq  = 0, 
than  with  the  finite  strain  interaction.  The  effect  on  a 
will  be  quite  small  since  it  is  representative  of  the  low 
velocity  limit  where  the  solute  profiles  are  changed  only 
very  slightly.  The  greater  effect  will  be  on  the  integral 
in  Eq.  (5) , and  so  is  more  likely  to  be  effected  than 

will  a.  This  means  the  greater  uncertainty  from  this 
approximation  will  be  on  the  conditions  for  breakaway  from 
the  'olute  cloud. 

Overall  the  agreement  between  the  experimental  data  and 
the  theoretical  calculations  is  improved  by  including  a 
strain  field  interaction.  This  suggests  the  basic  concept 
is  valid  and  warrants  further  development  of  the  theory. 

The  calculations  done  to  date  all  suggest  that  the  tempera- 
ture range  for  which  the  impurity  drag  would  be  effective 
would  increase  in  order  for  the  additives  Ba^"*",  Sr^'*'  and  Ca^'*’. 
Attempts  to  verify  this  experimentally  were  inconclusive  because 
of  unexpectedly  low  rates  of  nucleation  of  large  grains,  Bowen 
et  al  (1977) . 

C.3.  Breakaway  from  the  Solute  Cloud: 

From  the  steady  state  solution,  if  Mq  and  (1/aCoo)  are 
sufficiently  different,  the  drag  force  has  a local  maximum 


41 


4 


of  F 


Good  n“l  , ,u  , , , - 1 Zl  1/2 

p at  the  boundary  velocity,  V - 6 il+— — ^ 

max  ^ C 


max  " 2 

and  also  a local  minimum,  of  P 

V . (M  C a ) 

min  ' o «> 


min 


(aC^yM^)^/2  g-1 
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force,  F,  is  F^ 
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there  exist  three  possible 


ii^in  ‘max' 

velocity  solutions  for  Eq.  (3),  as  shown  in  Figs.  (C3)  and 

in  some  temperature  range.  The  intermediate  solution  is 
unstable  and  so  physically  uninteresti.ng . However,  in  this 
range  it  is  possible  for  a boundary  to  move  at  either  the 
slow,  drag  limited  velocity,  or  the  fast , intrinsic 
velocity.  The  time  dependent  problem  which  would  indicate 


the  conditions  for  achieving  either  mode  has  not  been  solved. 
It  is  anticipated  that  transitions  between  the  two  modes  may 
result  from  local  fluctuations  along  the  boundary.  If  the 
amount  of  solute  at  the  boundary  is  reduced  it  i.iay  lead  to 
a transition  from  the  slow  to  the  high  speed  mode.  Con- 
versely, if  a region  of  high  soeed  boundary  slows  down  it 
may  accumulate  impurity  and  then  shift  to  the  slow  velocity 
mode.  An  estimate  of  the  conditions  for  this  instability 
can  be  obtained  from  consideration  of  the  geometrical 
instabilities  which  can  result  during  the  motion  of  a planar 
boundary . 

If  a boundary  has  a waviness,  as  in  Fig.C5  the  local 
motion  of  an  impurity  as  it  is  dragged  by  the  boundary  will 
include  a component  perpendicular  to  the  average  direction 
of  motion  of  the  boundary.  This  lead  to  a tendency  to 
accumulate  impurity  in  the  slower  regions  of  the  boundary 
which  increases  the  drag  and  may  further  slow  such  regioriS . 
This  convective  effect  similarly  tends  to  remove  impurity 
from  the  more  rapid,  bulging  regions  leading  to  a local 
reduction  in  the  drag.  The  impurity  redistribution  (parallel 
to  the  boundary)  is  resisted  by  parallel  diffusion  of  the 
impurity, with  D',and  the  surface  tension  of  the  boundary,  y. 
From  a simplified  model  of  the  problem,  Roy  and  Bauer  (1975) 
showed  that  for  small  fluctuations  along  a boundary,  the 
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position  of  the  boundary  and  the  impurity  at  the  boundary 
can  be  expressed  in  terms  of  the  steady  state  values  for 
the  average  velocity,  Vq,  solute  concentration  at  the 
boundary,  Cq , plus  the  correction  terms  of  the  form: 


AX(y,t)  = [A^expiq^'^t)  +A2exp(q^  t)  ] cos 


^TTy 

A 


and 


AC  = [A  exp(q  ■^-c)  +A  exp(q  t)  ] cos  ^ 


(11) 
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where  are  amplitudes  of  the  pertubations.  By  modifying 
the  perturbation  analysis  to  include  the  nonlinear  drag 
force  and  excess  solute  concentration  at  the  boundary  it 
can  be  shown  that 
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and  in  the  high  velocity  limit, 

Y V„ 


Q = D + 


V ciC  V 
o °°  o 

M , . 2 5 C 


2„  2 


f 23 

f o o 


° I 


i+3^v  ^ r 

O J 


(15) 


It  is  obvious  from  Eq.  (13)  that  q is  always  negative 
and  cannot  contribute  any  instability  to  the  grain  boundary 
motion.  However,  where  q"*”  > 0,  i.e.. 


D < 


C a V ^ 

s , c °°  o 
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the  perturbations  in  both  the  shape  and  concentration  profile 
will  grow  and  instability  in  the  grain  boundary  motion 
will  result.  This  condition  of  instability.  Eg.  (16), 
holds  for  both  the  low  and  high  velocity  limits;  and  to  the 
irst  approximation,  it  may  be  correct  for  all  velocities. 

Experimentally,  it  is  easier  to  treat  the  driving 
force  F rather  than  the  dependent  variable  V . m the 
solute  limited  region  the  driving  force  is  balanced  by  the 
solute  drag  force.  Eg.  (3).  Consequently,  the  condition  of 
instability,  Eq . (16),  can  be  expressed  in  terms  of  the 
driving  force,  F,  as 


> a:  = 


(27t)  ^ a C y D 

<v»  • 


F^  (1  + 6^  V^^) 


(17) 


(2tt) 


~2  the  low  velocity  limit 

Essentially  for  long  wavelengths.  A,  at  which  q > 0,  any 
fluctuations  are  stable  and  will  grow  spontaneously.  For 
all  driving  forces,  F,  if  A > A^  dimensional  instability 
of  boundaries  can  result  leading  to  wavy  boundaries.  Further, 
^ min  ^ ^ < ^max'  geometrical  instability  may  lead 

to  pertubations  which  allow  the  transition  between  velocity 
modes.  The  explicit  inclusion  of  the  possibility  of 
evaporation  or  condensation  of  impurity  from  the  boundary  is 
expected  to  lead  to  a wider  range  of  instability,  i.e., 
smaller  A^.  if  a half  wave  segment  shifts  from  the  low  to 
high  velocity  regime,  it  may  spread  laterally,  Fig.cs  and 
carry  a large  segment  of  the  entire  bo  indary  into  the  high 
velocity  mode.  It  is  thought  that  for  these  intermediate 
riving  forces,  the  actual  velocities  of  boundaries  may  be 

distinctly  nonuniform  or  jerky  as  they  switch  from  one  mode 
to  the  other. 
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switch  to  high  velocities  could  lead  to  orders  of 
magnitude  higher  rates  o«  ornwei,  a ■ orders  ,f 

V,  ^ . growth  during  reorystalli zation 

V an  et  al.  ,1975a,.  Cannon  et  al..  ,i„e,.  xhe  Irregularly 

Shaped  grains  frequently  seen  during  low  temperature  re- 
crystallization  of  KCl  Fia  ra. 

five  of  th-  ' 9-  , re  thought  to  be  representa- 

ve  of  this  process.  Idiomorphic  grains  are  also  seen 
which  are  thought  to  resul-^  fT-r^Tr,  k j 

entire  ho  a boundaries  in  which  the 

entire  boundary  has  made  the  transition  to  intrinsic 

velocity,  .hus  for  a given  situation  Eg.  ,ly,  provides  a 

on  1 ion  on  the  amount  and  type  of  solute  necessary  to 

in  ibit  breakaway  and  rapid  recrystallization.  The  upper 

rmi  to  1 IS  of  course  the  physical  dimension  of  the  piece 

b«  as  a practical  matter  may  not  need  to  be  quite  L ' 
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ABSTRACT 

Conditions  under  which  four  possible  mechanisms  may 
control  grain  boundary  mobility  are  described.  Intrinsic 
mobility,  impurity  drag,  pore  or  particle  drag,  and  liquid 
pliase  control  are  considered  and  possible  mechanisias  for 
transitions  among  these  regimes  are  suggested.  The  particular 
mechanism  controlling  mobility  can  affect  the  kinetics  of  grain 
growth  and  the  result  ing  grain  morphologies.  A survey  of  the  data 
for  many  ceramics  indicates  qualitative  agreement  with  theory. 

INTRODUCTION 

The  importance  of  raicrostructural  control  in  achieving  the 
desired  mechanical,  electrical,  magnetic  and  optical  properties 
in  ceramic  materials  cannot  be  over-emphasized.  For  sintered 
materials  the  microstructural  evolution  depends  upon  the 
characteristics  of  the  starting  powder  (e.g.,  particle  size, 
size  distribution,  particle  shape  and  particle  agglomeration); 
the  green  compact  microstructure  (e.g.,  green  density,  binder 
material  and  the  state  of  particle  agglomerates);  and  finally 
the  sintering  and  coarsening  processes.  During  sintering,  pores 
are  removed  from  a ceramic  body  and  grain  bound  ries  develop 
between  particles.  If  a dopant  is  added,  the  dopant  ions  will 
further  react  with  the  major  constituents  and  may  bo  red istributed 
for  thermodynamic  or  kinetic  reasons.  The  grain  boundaries 
developed  during  sintering  act  as  vacancy  sinks  permitting  pore 
removal  and  densif ication.  However,  the  boundaries  tend  to 
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migrate  toward  their  center  of  curvature  to  reduce  the  total 
gram  boundary  area  and  energy.  Consequently,  the  kinetics  of 
grain  boundary  migration  have  a profound  effect  on  the  sinter- 
ing process,  on  the  achieveable  final  grain  sizes  and  densities, 
and  on  resultant  physical  properties.  Even  in  forged,  cast,  or 
CVD  ceramics  grain  boundary  motion  may  be  important. 

Although  the  boundary  mobility  and  the  microstructure  can 
be  significantly  affected  by  subtle  changes  in  the  dopant^  anc. 
porosity  levels,  the  mechanisms  of  grain  boundary  migration  are 
not  yet  well  understood.  A dramatic  demonstration  of  this 
variability  can  be  seen  in  the  collection  of  grain  boundary 
mobility  data  in  Fig.  1 for  many  ceramic  oxides  with  varying 
purities,  densities  and  amounts  of  second  phases.  For 

example,  the  measured  boundary  mobilities  in  MgO  differ  by  as 
much  as  10^  at  one  temperature. 


THEORY 


In  order  to  understand  data  (like  that  in  Fig.  Dl)  , the 
motion  of  grain  boundaries  is  first  considered  in  terms  of  four 
different  limiting  cases  which  are  shown  schematically  in  Fig.  E 
2,  and  are  referred  to  as:  intrinsic,  impurity  drag,  pore  drag, 
and  liquid  phase  controlled.  The  conditions  under  which  the 
various  mechanisms  may  be  controlling  are  then  examined, 
frequently  more  than  one  regime  may  be  possible  for  a given 
experimental  situation;  possible  mechanisms  for  transitions 
among  the  regimes  are  considered.  Finally  the  effects  of  the 
mechanism  controlling  mobility  on  the  kinetics  of  grain  growth 
and  the  resulting  grain  morphologies  and  pore  distributions 
are  discussed. 

^ The  intrinsic  grain  boundary  mobility  refers  to  motion 
limited  by  the  movement  of  ions  across  the  grain  boundary  of 
a width  of  a few  Angstroms  and  includes  the  processes  of 
detachment  from  one  side  of  the  boundary  and  attachment  to  the 
other  side.  Under  the  proper  conditions,  this  process  can 
control  the  bouncary  motion  even  in  typical  ceramic  materials 
which  have  impurity  or  nonstoichiometry  which  create 
atomic  defect  concentrations  greater  than  those  due  to  thermal 
energy  (e.g.,  Schottky  or  Frenkel).  Thus  the  use  of  the  term 
intrinsic  does  not  refer  to  a process  which  occurs  only  in 
pure  ceramics.  For  simplicity,  it  is  assumed  that  the 
impurities  within  the  grain  can  on  average  move  as  rapidly 
across  the  boundary  as  the  host  ions  and  that  they  do  not 
affect  the  grain  boundary  diffusivity.  Thus , although  it  need 
not  be,  the  intrinsic  mobility  is  taken  to  be  independent  of 
solute  species  and  concentration. 

The  major  effect  of  impurities  which  do  not  exist  as  a 
second  phase  is  shown  descriptively  in  Fig.D2.  If  the 
impurities  are  attracted  (or  repelled)  to  the  grain  boundary 
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concentratioiis  which  are  greater  (or  less)  than 
the  bulk  value  in  the  "near  grain  boundary  region",^  then 
these  ions  may  be  carried  along  by  a moving  boundary.  If 

' region,  26,  wider  thuii  the  boundary  core, 
a'vi  TV  controlling  their  motion  will  be  similar  to 

bulk  diffusivity  rather  than  a grain  boundary  diffusivity. 

The  motion  of  these  impurity  ions  may  limit  the  velocity  of 
the  boundary. 

tha  thickness,  £,  between  grains, 

the  transport  in  the  liquid  is  rapid  but  the  diffusion  distance 

will  be  ^arger  than  <d.  Two  limiting  processes  controlling 
boundary  motion  are  considered:  diffusion  through  the  liquid 
and  surface  reaction  kinetics. 

Finally,  the  fourth  situation  arises  when  a drag  force  is 
exerted  by  second  phases  which  may  be  crystalline  particles, 
iquid  droplets,  or  pores.  Tlien  the  motion  of  the  boundary  is 
frequently  limited  by  the  motion  of  the  second  phase  particles. 

Intrinsic  Grain  Boundary  Mobility 


From  absolute  reaction-rate  theory. 


. - . lurnhull  has  derived 

an  expression  for  the  velocity,  V,  of  a high  angle  grain 

width^’^w-^''*'^*^  driving  force,  F,  across  the  boundary 


V = 


kT  eJ)u)F  , . , Du  op 

IT  ’ “kT  ^^P(-AG./kT)=  (g) 


(1) 


where  AGi  is  the  free  energy  of  activation  for  the  migration 
process;  h,  Planck  s constant;  9.  is  the  ionic  volume;  Dk  is  the 
appropriate  grain  boundary  diffusion  coefficient;  kT  has  its 
usual  meaning;  and  e is  the  base  of  Naperlan  logarithms. 

In  a ceramic,  the  migration  of  a grain  boundary  requires 
the  simultaneous  diffusion  of  cations  and  anions.  If  the 
diffusivities  of  the  two  ions  are  different,  the  migration 
process  is  limited  by  the  diffusivity  of  the  slower  diffusing 
spec  es.  enerally,  the  cations  diffuse  much  slower  than  the 
anions  in  the  grain  boundaries;  then,  the  intrinsic  grain 

The  effects  of  grain  boundary  structure  may  be  important 
leading  to  mobilities  lower  than  those  given  by  the  Turnbull 
expression,  Eq.  (1),  to  anisotropic  mobilities  and  possibly 
to  impurity  effects.  The  structural  effects  envisaged  are 

atiin®co%J?  crystal  growth  models  in  which  an  accomod- 

ation coefficient  is  included  to  account  for  atomistic 
attachment  or  detachmentoccurring  at  only  a fraction  of  the 
atomic  Pltes;  e.g.,  screw  dislocation  growth. Gleiter^®  has 

-^-insic  grain  boundary  migration 
in  which  a specific  atomistic  model  was  assumed.  The  transport 
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of  an  atom  from  a kink 

shrinking  grain  to  a si-ilaV  sitp'n^^  “orface  of  a 

assumed  to  occur  by  several  -I  k?  ^ growing  grain  was 
dissociation  of  thl  s": "r'or  hT 

diffusion  along  the  sten-  atomic 

(d)  diffosio,,  on  tL  sSace  andT““‘°"  “'<='’• 

boundary  region;  and  (a)  drffusrof in"rh“°gra"“b'''b‘‘''‘‘‘" 
aasorption  on  the  surf^ipp  pP  -u  gram  boundary  and 

verse  attachment  st^  -- 

each  oJ'the'above  incorrect; 

relatively  slow  leading  to  ^ '^ns  assumed  to  always  be 

sum  of  those  for  all  tL  p^it’essL 

the  slowest  process(es)  which  ‘han  just  that  of 

across  the  boundary.  However  the  i t atomic  jump 

continued  in  the  expression  terestlng  features  are 

V = f ^ 

h bo  kT  ® (2) 

compared  to  Eq.  (irim'^in'^the  ^"^Portant  differences 

the  step  density  and  is  limited  trSTf 

Chat  f may  he  controlled  hv  ..  T ^ ~ i'leiter  argued 

under  certain  condit  oL  it  ^“"-^'^-^tions  and  that 

force  leadine  to  . n^T  ^ affected  by  the  driving 

driving  force.  The  activation^^^^^°"^'’J’’  between  velocity  and 

Slow  process(es).  but  there  may  b'"a^"dit?on\Vprf 

terms,  g,  to  account  for  the  othar-  Pre-exponential 

anisotropy  ray  result  fron  f such  thlrbr^^  ■* 
density,  f „ln  vanish  and  the 
Solute  Drag 

ions  and\;c:tat\‘old^a^™it\:t::  ttd“7" 

uniformly  distributed  in  the  near  emin  ^ a 
a stationary  grain  boundary,  the  dStJibuMn 
symmetrical  and  there  is  no  nar  a . ^ profile  is 

grain  boundary  and  the  solute  n k force  between  the 

As  the  boundary  migraLs  th"  boundary, 

asymmetric  and  this  resultf?  Oration  profile  becomes 

boundary.  The  derivatioi  of  a rate"' 

culation  of  the  imuuritv  r^^cf^  -r  ^ equation  entails  the  cal- 

uolute  telati.c  to  the  moving  J„L“?y""‘’lf'’th"°J'°"  ■ 

too  high,  the  boundary  will  brea"awav  fr  2”""’* 
cloud  and  the  velocitv  win  ^ from  the  impurity 

eiocity  will  approach  that  given  by  Eqs.  (1)  or  (2). 
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Impurity  drag  models  have  been  derived  by  Cahn^® , Lucke 
and  Stuwe  , Hillert  and  Sundman**^  and  others.  The  more  quan- 
titative and  concise  model  of  Cahn^^is  described  here.  This 
one  dimensional  model  is  based  on  the  following  assumptions: 

(1)  There  is  a dilute  solution  of  solute  in  the  matrix, 
thus  the  chemical  potential  for  the  impurity  is,p  = Pg  +kT£nC, 
where  C is  the  dopant  concentration.  2)  There  exists  an 
interaction  energy,  U,  between  th  ite  and  the  grain  bound- 

ary which  is  assumed  to  have  the  i ..  ^ing  properties;  (a)  it 
is  not  c function  of  velocity;  (b)  it  is  not  an  explicit 
function  of  local  dopant  concentration;  but  (c)  it  is  a 
function  of  the  spatial  coordinates  in  the  frame  of  reference 
of  the  moving  boundary.  (3)  The  Einstein  relation  between 
the  diffusion  coefficient  and  the  mobility  is  assumed.  (4) 

Only  the  steady  state  condition  is  considered;  thus  the  grain 
boundary  velocity,  V,  is  a constant. 

Because  of  the  convective  motion  of  the  impurity  relative 
to  the  boundary 


3t 


(3) 


The  chemical  potential  of  the  impurity  in  the  near  grain 
boundary  region  is 

U = kT£nC  + U(x)  +U  (4) 

o 

where  Uq  is  a constant.  The  impurity  flux,  J, as  observed 
from  the  reference  frame  of  the  moving  grain  boundary  is 


kT  dx 


(5) 


The  impurity  profile,  C(x),  can  now  be  calculated  from  the 
continuity  equation 


9x  3t 


0 


(6) 


and  the  boundary  conditions  at  x = «> : dC/dx  = 0,  dU/dx  = 0 and 
C(“)  = Coo.  Thus,  C(x)  must  satisfy. 


D 


^ dU 
kT  dx 


+ V(C  - C ) = 0 

oo' 


(7) 


The  drag  force  per  unit  of  solute  is  -dU/dx.  The  net 
impurity  drag  force  on  the  boundary,  F^,  due  to  C(x)  is. 


F 


i 


(C(x)  - C^] 


dx 


dx 


(8) 


where  Ny  is  the  matrix  atom  density.  The  C^x)  which 
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satisfies  Eq.  (7)  is  used  to  calculate 
(8).  A simpler  approximate  form  which 
low  velocity  extremes  is 

aC  V 

OO 

^'i  " 1 + g2"vT 


a drag  force  from  Eq. 
describes  both  high  and 


(9) 


wliere 


and 


u,  s kT 

V 


. ,2  .u(x)  ^ 

I sinh  (-L') 


_a  _ 

~ kT 


5 

r 


D(x) 


* dx 


d^ 

dx 


) D(x)  dx 


(10) 


(11) 


Physically,  a is  the  impurity  drag  per  unit  velocity  and 
per  unit  dopant  concentration  in  tlie  low  velocity  limit;  and 
B ^ is  the  drift  velocity  with  wliich  an  impurity  atom  moves 
across  the  grain  boundary  region.  From  the  form  of  a,  Eq.  (10), 
dopants  with  either  attractive  or  repulsive  interaction 
energies  of  equal  magnitude  will  exert  similar  drag  forces. 

The  Intrinsic  drag  on  the  moving  grain  boundary,  Fq,  is 
given  by  an  expression  such  as  that  from  Ti’rnbull,  e.g.,  Eq. 

(1)  or  (2). 


F 

o 


(12) 


where  Mq  is  the  intrinsic  boundary  mobility.  The  total  drag 
force,  Fig.  3,  is  the  sum  of  these  two  and  is  just  equaltothe 
driving  force  at  steady  state. 


F 


F 

o 


+ F. 
1 


(n) 


In  the  low  velocity  limit,  where  the  impurity  drift 
velocity  across  tliC  near  grain  boundary  region  is  high 
relative  to  the  grain  boundary  velocity,  6~'>>V,  there  is  a 
Linear  relationship  between  the  velocity  and  the  drag  force 


V = 


F 


f aC 


OO 


(14) 


This  is  the  experimentally  interesting  regime  considered  for 
the  impurity  drag  model  in  Fig.D2.  At  high  velocities  the 
amount  of  excess  impurity  near  the  boundary  is  very  small  and 
the  velocity  approaches  the  intrinsic  value.  Even  at  low 
velocities  the  amount  of  excess  solute  at  the  boundary  de- 
creases as  V increases.  From  Flg.D3,  it  can  be  seen  that  there 
is  a range  of  conditions  for  which  there  are  three  possible 
solutions  to  the  force-velocity  equation.  The  intermediate 
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.mi„porta„t.  Bowovep 

sip»  v'locu,  ’ 

must  b"  speclfiL”i.fitdor\“‘'  i'*!  ; they 

cular  sysL.  These  ?™c  a parti- 

metals.  In  ionic  materials 

impurities  and  charged  defec  s ar™L“?h“r  °/  '‘“<’''=>l™t 
trolled  by  electrosLtlr  ^ boundary  are  con- 

other  IntL-actJon^  f“  a b:  ? L‘,"a  “ 

forces  and  defect  disrrii  ^ halides  the  electrostatic 

calcularpH  • buttons  around  boundaries  can  be 

where  D(x)  L"u(xrcL^r®^  V^°T‘  P"°vides  one  case 

impurity  Lafca  cll'L  ‘ -P— d and  the 

include^secon^darJ^1?;:e^:,n.Lt  T,  \ , 

problem  has  not  been  solved  f^r-  ^ -.4  actions.  Although  the 
are  expected  to  e“"t  Interactions 

contln^uSItr’::  llZlZ  tsT  ’'t?“"“ 

mas  first  Postulated  by  F^^heCtand  Jhr  h 

have  been  formulated  by  Frenkel « f "llstributions 

and  Kliewer  and  Koehl„“  • Enholby  at  al“, 

because  of  the  differences  in  ^ field  develops 

:ijt“i:i-^aur‘r  t;r"  " 

:Subieia“L' 

Both  the  charged  vacincles  Ind  au“o“ale"t  “1^° 

charge  region  experience  an  eLtrostatic  We 

with  the  grain  boundarv  ind  interaction  force 

in  order  to 

diffusiuin,  n . • ....  energy.  Since  the  vacancy 

the  drag  f^ce  L r2“ng‘j  '"h‘"  J''®"  Blf fusivlty. 

3»aller*tha„  th:  imp:;u7dr7  - "->■ 

ions  :ay:L^ro“‘gLa^::i- - 

mismatch  betmeen  the  matrix  and  impurity  ^^1  e "he 

lattice.  For  aliovalent  ions  the  strain  F^  i j • ^ 

coupled  with  the  electric  field  h interaction  is 

excess  chartye  Th  . because  these  ions  have  an 

defects  within  the^^charL  charged  ions  and 

distribution  cafbe  s?gn"lfUa°narafL  t'^!'b“°r“‘‘^ 

or  Other  secondary  interactions.  y the  strain  field 
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For  a given  functional  form  of  the  qrrain  ^ 

terras  of  solnro  -na  a"r  ® equation  is  expressed  in 

f elds  are  shorn  i„  pig.Dd  for  KCl  contalnlog  50  ppm  Ca  Thi: 

er«i“;:rr[i:;ar-T,»n“<=‘> : 

between  iranurltv  in  strain  energy  difference 

:?e~  F 

tatic  potential  near  the  surface  as  corapared  to  the  bulk 
Ihe  maximum  value  of  the  tirmi,,  duik. 

be‘*°  strain  interaction  energy  was  taken  to 


U^(o)  = 


fArV 

L +V )l r 7 


radCs"  iL  Ar  "rhe‘“i''’  ’’“‘f'’"’*  ■•  the  matrix  ion 

enhanced  and  dip  eted  f;;“tr  "" 

cmmlr  i r-  . depleted  solute  concentration  because  of  these 
coupled  interactions.  Further,  even  near  the  isoelectric 
temperature  there  are  appreciable  electric  fields  near  tL 
boundaries  because  of  this  coupling. 

near  Ihl  diffusion  coefficient  has  a spacial  dependence 

cC  ra  lonr'  C"  h"  C ''--tion  in  the  defect  cL- 
ntrations.  Consider  Sr-doped  KCl;  the  Sr-ions  can  onlv 

move  by  associating  with  K vacancies.  The  probability  of 

association  in  the  bulk,  p^,  is  determined  largely  by  the 

between  the  differently  Carged 
solute  ion  and  uananrsir  TUi  ■ 1,  . 


solute  ion  and  vacancy 
coefficient  is  ‘ 


r.,,aigt.u 

i-he  specially  dependent  diffusion 


,.Sr  , , 

D U)  = 


^aiPcr,"^  C^exp  - 


e(^(x)  - ei|)(°o)  + [Us  (x)-Ug  (oo) 


Do  IS  the  diffusivity  of  Sr-vacancy  complexes, 
ivitv  thCdpf  literature  for  the  Sr-vacancy  diffus- 

“d  H^r:  i:Zf7„  If  T°" <“>■  “■'heevaluft'- 
derermlned  "f  fE,  ^ -e 

fix  Id  r "-5“™  -.perar  „re  for’l 

driving  force  of  105  dvne/C2  ^^‘'''  ‘'^‘^^^  is  shown  in  Fig.  D5.  The 

f ::s; 

energy  of  about  40  ergs/cra^.  Over  mucn  of  the  temperature  range 
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^<tOUCtO  DiSTrtNCt  ^ WOM  aOuNDANiji/ft) 


temperature,  °c 


Figure  D4.  Concentrations  of 
vacancies  and  Sr-impurity 
near  the  grain  boundary  of 
KCl.  The  range  of  the  attract- 
ive strain  field  was  asrsunied 
to  be  2 lattice  spacLngs. 

The  maximum  value  (-0.2eV) 
is  appropriate  to  Ca  in  KCl. 


Figure  D5.  Calculated  velocity 
of  grain  boundaries  in  Sr- 
doped  KCl  in  which  electro- 
static and  strain  energy 
coupling  are  both  included. 
The  intrinsic  and  impurity 
drag  regions  are  noted;  the 
dashed  line  is  an  unstable 
solution.  The  points  are  ex- 
perimental data.  For  these 
conditions  the  isoelectric 
temperature  is  444°C. 


motion  in  either  the  low  velocity,  impurity  drag  mode  or  the 
high  velocity,  intrinsic  mode  is  possible  for  this  driving 
force.  The  gap  in  the  impurity  drag  solution  results  from 
the  lower  electrostatic  attraction  near  the  isoelectric 
temperature . 

Recently,  Hillert  and  Sundman***  , developed  a solute  drag 
model  based  upon  the  evaluation  of  the  free  energy  dissipation 
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due  io  diffusion  when  a c?rair.  k 

volume  of  the  lattice.  Their  through  a unit 

solution  model  for  the  solute  activ^r"  ^ regular 

even  for  a high  solute  consent raM^  T therefore  valid 

may  be  required  in  the  more  compleToxijf  ’ formulation 

ferrites  and  titanates  in  which  soL  ^ 
components  commonly  occur.  At  ^ ^ solutions  of  several 
is  similar  to  the  Cahn  model.  ^ concentrations  this  model 

shows  that  Ler  a range^ord^Ling''fo°"  impurity  drag 

at  either  of  two  velocities  the  fi  boundary  may  mov£ 

would  indicate  the  conditions  .^/^“®/®Pondent  problem  which 
either  mode  has  not  been  sXd  remaining  in 

transitions  between  the  two  moH ' anticipated  that 

tions  along  the  boundary.  Tf  a^-irwlv'^m^^l  fluctua 

a region  where  the  driving  fsr  boundary  enters 

concentration  is  lowerth^n  t^r  " " 

locally  F>F^a^=aT/2B  ^'^lees  such  that 

the  hieher,  nearly  rntrinl^ryL^clty™''"," 

ition  it  would  shed  most  of  its  excess  i trans- 

if  a boundary  at  hieh  sneeH  . impurity.  Conversely, 

force  or  higher  solute  concentration  Inr°i'  driving 

accumulate  further  imnurltv  a u r ^ ii  ™ay 

These  transitions  need  not  o'"^  i’ll  t lo  the  low  velocity  mode , 
but  instead  may  star"  in  I oartl  T ^ boundary, 

ally  along  the  boundary.  The  degre^onr" • spread  later- 
in  F or  C.  sufficient  to  cauL  a transit 
IS  a tendency  for  imnuritv  rr*.  a ^ ff 

boundary.  Any  non-unifor^  i iscribute  laterally  along  a 
slow  boundary  would  r / redistribution  along  a 

lose  more  impurity  laterallv'"a  accelerate,  to  then 

and  perhaps  L maL  t^rfrsn  . surrounding  lattice. 

-le  OP ,,,  condi^^o;ff:^;h^"^°th::e^^t:i^r^^^• 

can  be  obtained  from  consideration  of  H can  occur 

ities  which  can  result  durine  geometrical  instabil- 

If  a boundarrhas  a l "V  boundary, 

perturbation,  as  in  Fig  D6  the^lo'^^^^  ^ 'velocity 

dragged  by  the  boundary  will  include"  ®°fote  ions 

to  the  average  direction  of' motion  of  th^r^^n'  cular 

impurity  toward  the  slower  regions  of  the  h moves 

accumulate.  Increase  the  drag  anH  f «here  it  may 

This  convective  effect  renin«  ’ • urther  slow  such  regions. 
hulSing  regions  lead'^g'^L 
impurity  redistribution  parallel  to 

by  diff  ,.ion  of  the  impurity  paJaUeJ  to  resisted 

surface  tension  of  the  bounL^v  v f the  boundary  and  the 

impurity  at  the  boundary  can  be  e^p^sseTirtel"  o^^h" 
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steady  state  values  for  the  average  velocity,  Vq,  plus  first 
order  correction  terms  of  the  form: 


and 


AX(y , t)  - exp  (q'^’t ) + A 2 exp(q~t)]  cos 

A 

AC  - [Aqexp(q+t)+  A4  exp(q"t)]  cos  — ^ 


(17) 

(18) 


The  time  exponents  of  the  fluctuations,  q'*'  and  q~,  are 
complex  functions  of  the  driving  force,  velocity,  drag  para- 
meters and  di  f f usiv  it  ies . Fluctuations  with  q>0  will  grow 
spontaneously  whereas  those  with  q<0  decay.  By  modifying  the 
perturbation  analysis  to  Include  the  nonlinear  drag  force,  Eq . 
(IJ),  it  was  shown  that  q>0  ‘^or  long  wavelengths,  X,  which 
satisfy: ' 


,2  (27r)^.tCooYD' 


(19) 


whi^re  D'  is  the  etfective  diffusivity  parallel  to  the  bound- 
ary . 

For  anv  drivimt  force  if  X dimensional  instability  of 
boundaries  can  result  leading  to  wavy  boundaries.  However,  if 
^min^  this  geometrical  instability  may  lead  to  a 

transition  betwr-n  velocity  mode.,.  A better  model  explicitly 
including  tiie  possibility  of  gain  or  loss  of  impurity  from 
the  boundary  trom  or  to  the  matrix  is  expected  to  lead  to  a 
wicer  range  of  instability,  i.e.,  smaller  . If  a half  wave 
segment  shifts  from  low  to  high  velocity,  it  may  spread 
laterally.  Fig. bo,  and  carry  a large  segment  of  the  entire 
boundary  into  the  high  velocity  mode.  For  intermediate 
driving  forces,  the  actual  velocities  of  boundaries  may  be 
distinctly  non-uniform  or  jerky  as  they  switch  from  one  mode 
to  the  other. 


Drag  from  Fores  ,md  Parti cles 

Porosity  is  t!u  mc.st  common  second  phase  in  ceramic 
materials.  Pore::  distributed  throughout  the  matrix  or  located 
at  the  grain  boundaries  restrict  boundary  motion.  The  dis- 
cre.te  nature  of  pores  leads  to  important  differences  between 
pore  drag  and  impurity  drag:  (1)  the  drag  force  due  to  pores 
on  the  moving  front  of  the  grain  boundary  is  non-uniform  as  - 
compared  to  the  more  nearly  uniform  impurity  drag;  (2)  the 
pore  mobility  is  frequently  much  slower  than  the  impurity 
diffusion.  Second  phase  particles  similarly  restrict  grain 
boundary  motion,  but  the  mobility  of  particles  may  be  lower 
than  that  for  pores. 

A pore  pinned  to  a grain  boundary  can  be  dragged  by  the 
boundary  at  a velocity  which  depends  upon  the  path  of  atom 
transport  around  or  across  the  pore.  The  expressions  for 
pore  mobility,  Mp,  are  generally  of  the  form 
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»V«  bV 


Vo  + bV 


Figure  D6.  Schematic'  representation  of  the  break-away  of  a 
grain  boundary  from  its  impurity  cloud  resulting  from  the 
growth  of  unstable  vciocity  pertubat  ions . 


/k'r) 


r 


The  exponent  ol  tlie  pore  radius,  r,  depends  upon  the  transport 
path  witli  n = 3 for  lattice  diffusion,  n = 4 tor  surface?  dis— 
fusion,  and  n = 2 or  3 for  vapor  transport,  and  AGp  is  approp- 
riate to  the  particular  process.  For  particles  vapor  transport 
is  impossible  and  interface  diffusion  may  be  slower  than 
surface  dif'^usionin  the  case  of  pores. 

Brook  t reated  the  kinetics  of  grain  boundary  motion 
under  the  influence  of  pore  interactions.  If  the  pores  are 
attached  to  the  grain  boundary,  the  boundary  velocity,  Vt,,  is 
equal  to  the  pore  velocity,  Vp,  giving 

Vp  = FpMp  = Vg  = Mb(F  - NFp)  (21) 
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from  wlilch 


V 


b 


= F 


N\+Mp 


(22) 


where  the  pores  are  oving  under  a driving  force,  F . The 
boundary  mobility,  includes  both  intrinsic  and  impurity 

effects,  Eq.  (13).  Part  of  the  driving  force  on  the  grain 
boundary,  F,  is  neutralized  by  the  drag  forces  from  the  pores 
which  have  an  areal  density,  N. 

Two  limiting  cases  can  be  observed: 


V 

V 

FM 

(23) 

FM 

F - NF  « F 

(2A) 

p nr 

The  driving  force  on  the  boundary  is  nearly  balanced  by  the 
drag  forces  from  the  pores,  and  the  boundary  motion  is  limit- 
ed by  the  pore  mobility. 


b) 


and 


NF 

P 


V =V. 


FMb 


F « 


F 


(25) 

(26) 


The  driving  force  for  grain  boundary  motion  exceeds  the  total 
drag  forces  due  to  all  pores  and  the  boundary  motion  is  limit- 
ed by  its  mobility. 

The  maximum  force  which  can  be  exerted  on  a spheri-al 
pore  by  an  attached  grain  bound:’. r ' is  ' rv,  VPie.r  t'  re.in 
boundary  moves  at  a velocity.  V,  > ^ ' Yj  '■>  .'eoar.ites 

from  the  pores.  The  uoper  Iniit  for  the  f^orce  necessary  ror 
pore-grain  boundary  par, ir  Lon  is,  from  Eq.  (21): 


F ■ • TTry 


(N  + ^) 

Mb 


(27) 


The  pore  density  on  the  boundary  may  be  much  higher  for  a 
slow  boundary  than  for  a fast  one  in  which  separation  has 
occurred  and  so  an  analogous  situation  as  t ^r  impurity  drag  is 
expected  where  for  a particular  driving  force  there  may  be  two 
stable  conditions  of  widely  different  V^,. 

Qualitatively,  three  different  modes  of  grain  boundary 
migration  can  be  expected.  For  larger,  less  mobile  pores  and 
low  driving  forces  the  pores  will  stay  attached  to  the  bound- 
ary and  control  the  motion,  ^or  higher  driving  forces  the 
boundary  may  separate  from  the  pores  and  its  motion  will  be 
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limited  by  Smaller,  iao’"e  mobile  pores  may  remain  attached 

to  the  boundary,  but  not  exert  a significant  drag  on  the  bound- 
ary; impurities  which  lower  Mb  extend  the  range  of  this  regime 
of  pore  attachment.  The  fact  that  Mj,  can  assume  two  different 
values  for  some  ranges  in  F further  complicates  delineation 
of  the  conditions  for  these  modes.  The  separation  from  pores 
and  break-away  from  impurities  may  be  correlated  in  some  cases. 


Grain  Boundary  Mobility  with  a Liquid  Film 

It  has  often  been  observed  that  the  rate  of  grain  growth 
is  rapid  when  a small  amount  of  liquid  wets  the  grain  boundary 
area.  The  rate  limiting  cases  which  have  been  considered  are: 
(1)  diffusion  of  the  bulk  ions  through  the  liquid;  and  (2) 
surface  reaction,  i.e.,  control  by  dissolution  or  precipita- 
tion. 

The  analysis  is  essentially  that  due  to  Greenwood ^ and 
Wagner  wlio  treated  a similar  problem,  the  coarsening  of 
particles.  Smaller  particles  have  a higher  free  energy  and 
thus  the  ions  dissolve  into  tlie  liquid  and  precipitate  on  the 
Surface  of  larger  particles  which  have  lower  free  energy  due 
to  lower  curvature.  In  the  coarsening  case  the  diffusion 
field  was  assumed  to  be  large,  i.e.,  a dilute  concentration  of 
particles.  Lay  corrected  the  expressions  for  the  finite 
*^ifftision  field  for  a liquid  film  around  grains. 

Wlien  diffusion  through  the  liquid  layer  of  thickness  Z 
(Fig.  1)  is  rate  limiting 

fiF 

" ~ Z W (28) 

where  the  motion  results  from  the  difference  in  solubility,  C£ 
:uc  to  curvature,  AG£ = C„ftF/kT.  Although  and  Db  may  be  of 
comparable  orders  of  magnitude,  if  the  film  thickness  is  in 
the  range  10- - 10^  A the  mobility  is  much  less  than  the 
intrinsic  mobility,  Eq.  (1). 

The  process  may  alternatively  be  controlled  by  the  sur- 
face reaction  rate,  Wagner. If  it  is  slow  relative  to  Eq . (28) 


The  reaction  constant, K, is  generally  exponentially  dependent 
on  temperature  but  also  depends  on  the  precise  kinetic  process 
and  may  be  dependent  upon  orientation  and  driving  force. 

Grain  Growth  Kinetics 

Boundary  motion  can  result  from  the  free  energy  difference 
between  phases  during  a solid-solid  phase  transformation,  from 
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strain  energy  during  recrystallization,  or  simply  from  grain 
boundary  curvature  during  grain  growth.  In  prLciple.  tie 
overall  gram  growth  kinetics  can  be  predicted  from  the 
forcp^i^  velocity  relations  and  specification  of  the  driving 
orce  in  terms  of  the  evolving  microstructure.  A statistic- 
ally averaged  boundary  velocity  is  assumed  to  be 


V = iiG 

2 dt 


(30) 


Tlie  total  driving  force  available  can  easily  be  obtained  from 
the  density  of  grain  boundary  energy  as  “‘stained  from 

= ll  , 3x 
L G 


F = yS 

t V 


(31) 


where  Sv  is  the  boundary  area  per  unit  volume,  L 
inear  Intercept,  and  G the  average  grain  size  ^ 
growth  rate  is  then 


the  mean 
The  maximum 


— = 6My 

dt  G (32) 


•V  X JL  J 


rtowever,  the  net  driving  force  for  growth  of  most  grains 
be  less  than  that  predicted  by  Eq.  (31)  because  the  actual 

fir  fhp  Shown  that  a reasonable  approximation 

curvature^is  which  takes  into  account  the  actua] 


(33) 


where  n is  a factor  which  is  %5/4*.  Grains  smaller  than  the 

shJinra  ^ negative  curvature  and 

shrink  and  eventually  disappear;  those  larger  than  G,  have  a 

net  positive  curvature  and  grow.  After  a transient  stagl 
during  which  abnormal  grain  growth  may  occur,  a steady  state 

f this  so  called  normal  growth  the  distribution 

relative  sizes,  G/G^.  remains  constant  and  there  are  no 
grains  with  size  G > 2Gc-^  The  average  size  G = 8G(,/9  increase 


as 


dG  _ ^ 
dt  G 


(34) 


The  average  driving  force  and  growth  rates  are  about  1/6  that 
predicted  from  the  simple  model,  Eq.  (32).  However,  if  the 
growth  IS  abnormal  with  only  a few  large  growing  grains  the 
^qwth  rate  of  ^these  large_grains  is  given  by  Kg.  (32) 

G«G  In  the  limit  of 

b«Gc,  n-1;  for  G»Gc,n=3/2  from  Eq.  (31). 
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Wlien  pores  or  particles  are  present  and  are  relatively 
immobile  (i.e.  , I>JMb»Mp)  they  can  pin  boundaries  and  reduce 
tie  driving  force  sufficiently  that  a limiting  grain  size 
results.  Zener  has  shown  that  this  occurs  when  the  de- 
crease in  grain  boundary  area  from  grain  growth  is  balanced  by 
the  increase  in  area  as  boundaries  separate  from  the  pores. 

For  a uniform  size  distribution  the  limiting  grain  size.  Gi, 

IS  8r/3v,  where  v is  the  volume  fraction  of  pores.  The 
pores  provide  an  effective  back  stress  resisting  the  motion 
of  any  boundary  so  that  if  there  is  a distribution  of  grain 
sizes  the  gi owth  rate  is  approximately^^ 

^ • G (a) 


4r)M^>  (- 


<-•  1 


dt 


G,+G 
1 c 

c.(; 

l r 

"l-^^c 
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G,G 
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I c 
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Y(fT  -p  - ),  G > 


(35) 


(c) 
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G,  G > 
c 


1 


(d) 


Ihe  powth  rate  is  essentially  zero  for  those  grains  on 
which  the  net  driving  force  is  insufficient  to  tear  the 
boundary  away  from  the  particles.  Hillert^"  showed  that  for 
small  initial  grain  sizes,  normal  growth  can  occur  but  the 
rate  will  decrease  as  the  limiting  size  is  approached.  This 
asymptotic  decay  may  be  approximated  by:” 


dG 


(; 


dt 


(1-7^) 


(36) 


If  the  size  distribution  is  narrow  the  final  grain  size  G 
may  actually  fall  1„  the  range  | Gf  to  G,  becfuse  the 
larger  grains  become  immobile  at  G<Gj^,  i.e.,  Eq.  (35b). 

Normal  growth  can  only  continue  if  Gi  increases  by 
dissolution  or  coalesence  of  the  pores  or  particles.  If 
coarsening  occurs  by  diffusion  between  particles  i e 
Ostwald  rinening,  there  .ill  be  a slow  growth  of  r, ’which  may 

acpendence,  . allowing  a similar 
growth  of  Gc  and  G^.  This  vill  occur  on  a much  longer  time 
scale  than  growth  below  Gi.  Alternatively  pores  or  particles 
can  be  dragged  by  the  boundaries  which  also  results  in  much 
slower  growth  controlled  by  M As  the  smaller  grains  dis- 
appear the  pores  on  these  boundaries  will  coalesece”  leading 
to  an  increase  in  r and  Gl.  Coalescence  by  particle  drag 
will  be  more  rapid  than  by  interparticle  diffusion  only  at 
small  concentrations,  v,  of  second  phase. 


66 


For  normal  grain  growr.h  integration  of  Eq.  (34)  gives 
the  familiar  parabolic  growth  relation 


G^(t)  - G^(o)  = 2Myt 


(37) 


This  results  only  if  the  controlling  mobility,  Mb,  is  equal 
for  all  boundaries,  as  may  occur  under  proper  conditions  for 
the  intrinsic  mobility  or  for  solute  drag  in  the  low  velocity 
limit. 


^ ^7  > grain  growth  data  are  found  to  fit  the 

relation  G (t)  - 

be  inferred  that 


^ ivi  t-  w J.XU  L.11C 

C^(o)  = Kt.  If  the  growth  is  normal  it  may 


V = 1 dG 

2 dt 


K 


2nG 


n-1 


(38) 


further,  if  the  lower  driving  force,  Eqs.  (33)  and  (34),  is 
appropriate,  it  may  also  be  Inferred  that 


K 

nY^~ 


2 


(39) 


For  some  mechanisms  the  magnitude  of  the  r.obility  will  change 
with  grain  size  as  is  implied  by  Eq.  (3^)  if  nf  2.  When 
growth  is  controlled  by  particle  drag  and  coalescence,  the 
increase  in  r as  G increases  will  cause  Mp  cc  decrease,  Eq. 
(20).  Similarly  for  a boundary  with  a thickening  liquid  film 
M may  change  with  G,‘"i.e.,  Eq.  (28).  Cubic  growth,  n=3,  can 
obtain  for  impurity  drag  control  if  the  amount  of  solute 
segregated  at  the  boundaries  is  much  greater  than  that  in  the 
grain  interiors.  This  is  most  likely  to  occur  at  lew 
temperatures  and  fine  grain  sizes  with  impurities  whica  are 
•Strongly  attracted  to  the  boundaries.  Transitions  between 
impurity  drag  limited  motion  and  the  intrinsic  velocity  will 
result  in  various  boundaries  moving  with  much  different 
mobilities  which  will  invalidate  the  assumptions  and  detailed 
results  of  Hillert's  theories,  such  as  Eqs.  (34)  and  (36). 
The  rate  and  steady  state  size  distribution  of  normal  growth 
will  be  altered  by  rapid  disappearance  of  the  smaller  grains 
if  the  boundaries  breakaway  from  tiie  impurity  cloud. 
Alternative!)' , the  susceptibility  to  rapid  abnormal  growth 
will  be  increased  if  the  boundaries  of  the  larger  grains  can 
breakaway  from  the  solute  and  grow  at  high  velocity.  Care 
must  be  taken  to  distinguish  steady  growth  with  n=3  or  4 
from  asymptotically  decaying  growth  of  the  type  predicted  by 
Eqs.  (35)  and  (36).  An  analogous  decay  will  result  from  a 
shift  with  increasing  grain  size  from  growth  at  the  intrinsic 
mobility  to  control  by  impurity  drag. 

Generally,  during  sintering  of  ceramics  all  the  pores 


67 


are  initially  on  the  boundaries  so  the  grain  size  is 
essentially  at  the  limiting  size,  i.e.,  G - Gi.  Normal 
growth  can  occur  by  pore  drag  or  by  growth  of  Gl  as  farther 
sintering  reduces  the  volume  fraction,  v,  of  pores. 

Coalescence  by  diffusion  between  isolated  pores  will  not  be  as 
important  as  further  densif Ication^^  unless  pore  shrinkage  is 
inhibited  by  entrapped  gas  or  a low  dihedral  angle.  If  G^  is 
somewhat  less  than  Gi,  normal  growth  may  be  slow  relative  to 
abnormal  growth  if  there  are  a few  grains  large  enough  to 
have  sufficient  driving  force  to  separate  from  the  pores, 

Eq.  (35c).  Abnormal  grains  can  grow  to  much  larger  sizes 
than  Gl  with  the  final  grain  size  determined  by  the  number  of 
these  large  grains  which  eventually  impinge  and  remove  all 
grains  smaller  than  Gl. 

The  tendency  for  abnormal  growth  and  pore  entrapment 
will  be  significantly  affected  by  solute  drag.  If  there  is 
sufficient  solute,  continued  dens:  i'ication  will  result  in  a 
shift  to  impurity  drag  control  when  v and  r become  sufficiently 
small  such  that  NMb  « Mp.  The  force  for  pore  separation, 

Eq.  (27),  will  then  be  larger  than  the  case  in  which  Mb  is 
high  as  assumed  in  Eq.  O'^c)  and  pores  will  tend  to  stay 
attached  to  the  migrating  boundaries.  Abnormal  growth  may 
then  involve,  but  not  necessarily  require,  breakaway  from  the 
impurity  cloud  and  separation  from  pores.  If  a portion  of 
oundary  breaks  away  from  the  impurity  cloud,  it  may  have 
sufficient  driving  force  to  tear  away  from  pores  since  Mb  will 
then  be  significantly  increased.  The  occurrence  of  an 
abnormal  grain  with  hounaaries  which  are  solute  free  but  which 
are  carrying  pores  is  possible.  However,  its  growth  is 

likely  to  be  limit;  by  further  accumulation  of  pores  as  it 
migrates . 

The  V ?locity  fluctuation  method  of  breakaway  from  the 
solute  cloud  can  provide  a path  for  nucleation  of  abnormal 
grain  growth.  This  mechanism  will  become  important  when  the 
critical  waveleng'h,  Eq.  (19),  becomes  comparable  to  or  less 
than  the  mean  pore  spacing  along  the  boundaries,  i.e.,  Ac<s. 
During  sintering  the  spacing  between  pores  grows  as  v 
decreases  and  as  d increases  by  normal  growth  with  pore 
coalescence.  At  the  critical  stage  that  the  spacing  exceeds 
Ac  the  microstructure  becomes  unstable  with  respect  to 
separation  of  the  grain  boundaries  from  the  impurity  cloud 
and  from  the  pores.  The  ensuing  secondary  grain  growth  and 
pore  entrapment  may  stop  further  sintering.  Increasing  the 
solute  concentration  will  increase  the  resistance  of  the 
boundary  to  breakaway  from  the  impurity  cloud,  Eq.  (19),  by 
increasing  Ac,  as  well  as  lowering  M|j  in  the  drag  regime,  Eq. 

U4).  Additives  in  excess  of  the  solubility  limit  may  give 
lower  values  of  s because  of  precipitate  pinning  as  well  as 
give  ma>imum  value's  of  A^  and  minimum  values  of  Mb  all  of 
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which  further  stabilize  the  body  to  boundary-pore 


separation. 


EXPERIMENTAL  DETERMINATION  OF  MOBILITY 

D8.  Mos^o?  Figs.Dl,  D7,  and 

qtMHivfa  were  calculated  from  grain  growth 

til  - n i^intered  or  hot  pressed  samples.  For  all  of 

within  ergs/cn,2  which  is  generally 

within  a factor  of  -.wo  of  the  true  values.  For  the  relativelv 
in  requent  cases  ol  parabolic  growth  Eq . (37)  was  used.  When 

from'^Er'^nSl"  r ^ ^ mobility  was  calculated 

To  faciiitar  implies  a grain  size  dependent  mobility. 

Ld  to  mobilities  were  normal- 

ize^d  to  that  for  a gram  size  of  90nm  (equivalent  to 

Ft  - 10^  dyne/cm2).  There  were  a few  reports^  abnormal 

gram  growth  from  whic^^^  ^^bility  could  be  calculated  from 

single  c^ystaL  ' measured  in  hot  forged 


discussion 


Grain  Boundary  Mobilities 


Many  grain  growth  studies  reported  in  the  v 

by-products  of  sintering  studies-  1^=  literature  are 

ara  tor  porous  oamplos.  Evof-!uh  Jullv  d!!!,®”"" 

-rorlals  aro  Ukely  to  Po  cop^”  „ “ ^ 

can  significantly  affact  boundary  .obllity  mZ,!  tw 

s -pLi-d 

comparisons  show  Important  differences  This  is  as 
affectively  by  comparing  the  moblluLs  as  a fu^cttror” 
reciprocal  homologous^tempe-ature  (T/Tm)-1,  Figs.  D7  and  D8. 
n,„u.  , definitive  measurements  of  the  intrinsic 

known  TM  which  Db  for  t he  slow  io  i is  also 

expression  Eq  (1)  f validity  of  the  Turnbull 

rZV  consistent  with  Eq  (ifusiL  a 

reasonable  estimate  of  Dfe.  In  Fig.o8  estimates  of  M„ , Eq.(l),are 

U02,  wDb  of  U was  measured  by  tracer  diffusion.'®  The  wDh 
MiO  co^L'„tra?L"fbut  to^o  S If  ““f  "" 

m nt  T,  »aeed’’on  recent  ^rg;. 

ments  which  discounted  wide  boundaries  and  space  charge 
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Figure  D7.  Grain  boundary  mobilities  for  oxides  and  halides  as 
a function  of  reciprocal  homologous  temperature  (Tm/  ^-1 
The  additive  concentrations  are  in  cation  fraction  unless 
otherwise  noted. 


contributions  to  wDb.‘**^  The  significant  differences  in  Mo  for 
different  compounds  with  different  crystal  structures  may  be 
eoretical  interest.  All  of  the  measured  mobilities  for 
oxides  are  well  below  these  estimates  for  Mq. 

In  general,  the  mobilities  for  porous  samples  are 

MgO  aircon  'h  • dense  ones,  e.g.,  103-10^  slower  for 

MgO  and  CoQ  which  is  consistent  with  particle  coarsening  or 

sintering  the  porosity  varied  both 

SamantranJ  ^ob/'^  T ^^e  topology.  Gupta“  and 

Samanta  and  Coble  have  shown  that  a linear  relationship 

frequently  exists  between  density  and  grain  size  during  inter- 

porosity  is  continuous  and 

the  densification  and  particle  coarsening  are  obviously 
y interrelated.  ->e  co'arsening  may  be  limited  by 
surface  diftusion  or  vapor  transport  to  the  free  surfaces  or 
y ensification  rather  than  by  grain  boundary  motion.  “ 

lat^H  85-95%  density  when  the  pores  become  iso- 

lated either  normal  or  abnormal  growth  may  be  observed 
depending  upon  the  size  and  distribution  of  pores  and  impuri- 

pirtl-J  generally  higher  than  the 

partieLe  coarsening  rates.  For  many  of  these  the  grain 

5»  be_3  but  there  are  reported 


CaO,^  MgAl204,^i 


P°^oes  samples  of  BeO,^^ 

NiO,  and  U02  in  which  the  mobilities  were  low!  The  value 
o n oes  not  necessarily  distinguish  between  grain  growth 
models.  In  UO2  grain  growth  with  n=  3 and  with  pore  growth 

ttributed  to  pore  drag  controlled  by  vapor  transport  with  the 

coaLL°"ce  8*='“  <P  ■ Il'/t)  »sU  « 

. , dense  samples  with  second  phase 

precipitates,  such  as  MgO  doped  AI2O3,  " the  low  mobilities 
are  consistent  with  particle  drag  control. 

It  IS  suggested  that  the  mobilities  inmost  dense  samples 
were  impurity  drag  1 imited,  although  small  amounts  of  porLity 
ware  son, atin.cs  important.  Forinatance,  in  MgO  paraboUc  ^ 
growth  was  found  in  some  undoped  samples however  the 
mo  ilit-ies  are  20-100  times  slower  than  Mq  and  the  samples 
contaxned  appreciable  impurity  (»  100  ppm)  suggesting  these 
are  impurity  drag  limited.  Additions  of  Fe2ofor  smfu 

reductions  of  M.  ^ 

theoretically  dense,  hot  pressed  MgO  high  initial  growth 
ates  were  found  which  decreased  rapidly  due  to  the  dLelop- 
ment  of  small  amounts  (<  1%)  of  po  rosity  ^ 

f compelling  test  of  the  impurity  drag  theory  is 

mobiUtv"'"  significantly  reduces  the 

drag  has  been  explicitly  calculated,*'®®  Fig.Db.  The  ordL  of 

i^tir^h  encouraging,  but  further  refinements 

in  the  theor>  are  necessary.  In  the  absence  of  sufficient 
information  to  calculate  a and  6 from  Eqs.  (10)  and  (11)  the 
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FiguroD9.  CoU  alter  various  amounts  of  grain  growth  showing 
1 he  precloiiii  naiil  ly  abnormal  mode  of  growth  with 
impingement  to  ultimately  give  a narrower  distribu- 
tion of  grains;  averagt*  densitv  99.3%  (courtesy  of 
LI.  Ciiowdhry). 

the  detailed  time  dependence  cannot  be  expected  to  follow  any 
simple  tl/>i  law.  The  pattern  of  alniormal  growth  may  be 
ultimately  related  to  diftercnues  in  porosity  distribution 
caused  by  the  agglomerate  history  of  the  starting  powder. 

The  grain  shapes  observed  in  KCl  during  recrystallization 
provide  an  indication  ot  itie  particular  growth  mode.  The 
grains  in  Fig.DlO(c)  are  thought  to  be  growing  in  the  impurity 
drag  limited  range  and  apparently  do  so  rather  uniformly.  In 
contrast  the  straight  sided  grains  seen  by  the  authors,  Fig.D 
10(a),  and  others  '®  in  pure  KCl  result  from  rapidly  moving 
boundaries  at  low  tempi ritures  which  have  broken  away  from  the 
impurity  cloud;  there  are  thought  to  be  moving  essentially  at 
the  intrinsic  velocity.  The  faiets  may  be  the ^^esult  of 
anisotropy  of  Mq  as  is  predicted  by  the  (.leitet  > model,  Eq. 

(2)  , where  the  atomic  step  density  in  the  boundary  and  the 
mobility  are  lower  for  bouiutai ics  which  are  parallel  to  a 
low  index  plane  in  one  '^rain  Finally,  the  very  irregular 
grain.  Fig.  DlO(b)  is  thought  to  have  occurred  because  parts  of 
the  boundary  had  broken  away  from  the  impurity  cloud  by  the 
fluctuation  mechanism.  Note  the  long  lingers  which  can 
result  from  growth  without  spreading  of  one  segment.  .Mso 
note  the  regions  of  unrecrystal  1 ized  material  which  is 
apparently  entrapped  between  adjacent  breakaway  segments. 

Straight  sided  grains  frequently  result  from  rapid 
secondary  grain  growth  in  the  presence  of  a liquid  phase, 
such  straight  grains,  instead  of  cusped  boundaries , resulting 


low  velocity  drag  can  only  je  estimated 


very  crudely  from^ 


OO 

ulTtli  JatMc  diffusivity  ^usually  more 

"Se?Sv  ^ft»«‘va'i„teracur 

gy.  There  are  only  a few  studies  in  which  the  expected 
decrease  in  mobility  with  higher  additive  concentratiorJf 

aSdi^  ’ ^203.”  For  instance  in  Y2O3  higher 

in  anot^  f°i^nd  to  reduce  growth  rates  J'\ut 

mohn?^  ’’  ^ additions  of  TIO2  or  Zr02  increased  the 

increo.?^  relative  to  that  in  pure  samples.  This  apparent 
on  L^n!?  h”  combined  effects  of  the' additive 

difr  'mobility  as  well  as  on  the  solute 

diirusivity.  In  certain  ranges  of  extrinsic  behavior  where 
I tne  drag  may  be  independent  of  the  actual 

..purity  concentration.  “ Frcnently.  various  additlvL  have 

Which  a^r^t 

More  complex  solid  solutions  such  as  the  ferrites  and 
s mobilities  even  though  in  some  cases  (PLZT)  the 

-samples  were  single  phase  and  theoretically  dense  (v«  0.n)^  " 

This  may  be  an  indication  of  a low  intrinsic  mobility  due\o 

clTTetTci  in  few  sites  from  which  atoms 

can  detach  and  jump  across  the  boundary,  e.g,  Eq.  (2). 

than  th  phases  are  present  the  mobilities  are  higher 

than  those  for  porous  samples,  but  often  not  much  higher  tLn 
dense  samples.  Lay  has  shown  satisfactory  agreement 
between  growth  rates  in  UO2  + 2.6%  AI2O3.  in  which  n=3,  Ld 
magnitude  predictions  of  the  d.^ffusivity  and 
nha^  assuming  control  by  diffusion  in  the  liquid 

phase.  When  the  liquid  has  a high  dihedral  angle  and  does 
not  fully  penetrate  the  grains  the  growth  rates  are  slower.'" 

Microstructure 

Very  few  grain  growth  studies  in  oxides  fit  simple 

InZ.f  Microstructural  examination  of  sintered  materials 
ndicates  that  one  reason  for  this  is  that  the  growth  is 
frequently  not  normal.  In  Fig.D9  the  microstructure  is 

foHhtc^r  CoO  samples 

for  which  the  as-hot  pressed  grain  size  was  l.Su,^®  The 

initially  relatively  uniform  raicrostructure  (not  shown) 

^r^nor^L^H-T'^^K^  impingement  approached  a 

est  Although  the  mobility  can  be 

estimated  from  such  microstructures  using  Eq.  (32)  or  Eq.  (39), 
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tiguro  DIO.  Various  morphologies  of  recrv;: 
forged  and  annealed  Kri  .single 
t ine  subgrains  from  the  deform 

from  abnorniai  gr<.)wtli  have  beei*  elaimed  to  b< 
presence  of  a licpiid  piiase.  These  idioiiuir 
pure  KCl  indicate  tliat  there  are  other  cau; 
sided  boundaries. 

The  abnormal  growth  mode  as  well  as  ( 
features  are  seen  in  t'ue  microstrticture  of 
annealed  Fe304/^  Fig.  Dll.  The  large  secoi 
left  contains  entra|)ped  pores  as  are  frequi 
also  has  an  entrapped  ttrain.  A.  nerbans 


SUMMARY 


Models  for  the  grain  boundary  mobility  in  ceramics  are 


Figure  U1  ] , Hof  lot.sjtH  ouH  mnoalt'd  FC3O4  which  shows 

>.  ufiMpiK' I |)va-('s  , ; nine  in  non-random  arrays,  an 
onlra|'|u‘d  grain,  and  ]iredoininantly  abnormal 
growl  li;  ample  ih-n.sitv  H9.6%  (courtesy  of  J. 

Ha  1 J oran) . 

dibcussi-cl  for  tlu  cases  of  couIidI  by  Che  intrinsic  mobility, 
by  drag  ol  solute  atoms,  Hv  tlrag  of  pores  or  particles,  and 
for  difinsion  tliriugh  .a  lupiid  fi]m.  Tlie  conditions  for 
whic'li  flu  diilercnt  m cli  ai  ■ im>  ait  oj'orative  have  been  indicat- 
ed. Ovei  wide  rangt  s drivii.g  foice  a boundary  may  be  drag 
limited  and  movt  at  a low  .apeed  <>r  if  sei>arated  from  the 
solute  (loud  or  poie  if  ina  ■ movi  nuu  h faster.  A simplified 
mechanisiii  for  breakaway  from  the  itii|5urity  cloud  is  presented 
and  the  pos.sibility  of  ( ciopt  r at  i vo  separation  from  impurities 
and  pores  i.s  discussed. 

(•rain  growth  kinet  i t ( ■' irv  lor  different  migration 
mechanisms.  The  tiiiil ' i v.a  lue^d  fspt  i t of  the  mobility  can  lead  to 
complets  growth  kinet  it  ineluuing  dor  tying  growth  rates 
caused  by  a transition  in  I lie  limiting  mechanism.  The  condit- 
ions for  abnormal  growth,  with  or  wMtliout  pore  entrappment, 
are  also  affected  by  such  transitions.  This  is  of  obvious 
importance  in  sintering  as  [trevention  of  separation  of 
boundaries  from  pores  i.s  important  ii  dense,  uniform  micro- 
structures  are  to  be  produced. 

There  are  very  few  measurements  of  the  boundary  mobility 
in  ceramics  wliich  unaml.igiiously  show  the  distinction  between 
solute  and  second  pliase  effects  or  wtiich  provide  a quantitative 
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confirmation  of  any  of  the  mobility  or  grain  growth  models 
Fron,  a comparison _ of  the  trends  in  the  data  wUh  the  Se^ 

“(IWhrr^r  “^."'^'■^“^tructures  it  is  concluded  that 
(1)  The  best  confirmation  of  the  solute  grag  theory 

Sr  dopefKCl^si''^''^^"'''f mobility  calculations  for 
^ incorporating  electrostatic  and 

experimental  measurement^ in 

usine  estimate  of  the  intrinsic  mobility 

mobUitJes  evl'nii  tha..  observed  bounjlr) 

velocities  .ay  be  c"°”e";y’'lSp;rltrdrS" 
appear  to  approacb  I„„I„sL  v^locILfrf^ 

suggesti  tha°r  ““ta  for  oxides 

(a)  pores  or  particles  generally  exert  a larger  drae 

irhiRh  as  particularly  when  the  porosity 

IS  high  as  during  sintering.  ^ 

hlgj’r«‘s‘“'^r““  ‘ “1"  Sive  relatively 

high  rates  of  grain  growth.  ^ 

(c)  the  boundary  mobility  in  materials  with  simple 

structures  and  compositions  is  higher  than  that  for 

more  complex  materials,  e.g.,  MgO  vs.  PLZT  or  ferrites 

(4)  The  shape  of  grains  can  be  affected  by  the  control- 

1 mg  mechanism,  e.g.,  faceted  grains  may  result  from  growth 

at  the  intrinsic  mobility  or  from  liquid  phase  co^troUer 

growth  particularly  if  it  is  interfa  ce  reaction  controlled 

causrve?rirrerr  intrinsic  control  can  ' 

cause  very  irregularly  shaped  grains. 

Is  lrrcSa^"‘?M®'‘‘^"  boundaries 

ant  s^Mal  var?  t"  ““  depead- 

ent,spaclal  variations  in  the  driving  force.  Thus,  the 

boundaries  of  small  and  large  grains  have  higher  velocities. 
boundary°motion“^^"  mobility  mechanisms  can  also  cause  jerky 
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E.  HOT  FORGING  OF  CaF^ 

E.l.  Introiuction : 

Or  the  various  classes  of  materials,  the  alkaline  earth 
halides  have  among  the  lowest  absorption  in  the  IR  spectrum. 
This  makes  them  a candidate  material  for  use  in  the  CO  laser. 
The  mechanical  properties  of  these  materials  has  become  a 
limiting  factor  in  their  utilization  in  high  power  laser 
systems . 

Hoc  forging  of  single  crystals  of  alkali  halides  has 
proved  to  be  successful  in  improving  mechanical  properties 
without  degrading  optical  properties.  As  a result,  deforma- 
tion, subgrain  formation,  and  recrystallization  behavior  has 
been  studied  in  such  alkali  halides  as  KCl , KBr,  AgCl , and 
NaCl.  Data  of  tiiis  nature  for  alkaline  earth  halides  is  at 

best  sparse.  It  was  the  oojective  of  this  work  to  determine 
this  information  for  CaF^ . 

E.2.  Materials: 

Single  crystals  of  CaF2  were  purchased  from  two  vendors. 
All  crystals  were  in  the  shape  of  right  circular  cylinders 
with  a height  to  diameter  ratio  of  1 to  1.  All  crystals 
purchased  from  Harshaw  Chemical  Co.  were  "pure"  and  un- 
oriented.  Two  sets  of  crystals  were  purchased  from  Optovac; 
a)  pure, oriented  (<100>,  <110>,  and  <111>)  and  b)  doped, 
oriented  (100  ppm  Gdy  <110>  ). 


E.3.  Experimental  Procedure: 

All  samples  were  heated  and  subsequently  deformed 
under  vacuum.  A graphite  die  was  used  for  all  runs. 

Grafoill  was  used  to  reduce  friction  between  the  sample  and 
die  piston.  An  internal  pressure  of  less  than  lOy  was 
maintained.  Samples  were  heated  at  a maximum  rate  of  200°C/h 
and  allowed  to  soak  at  forging  temperature  for  at  least  30 
minutes  to  minimize  temperature  gradients  within  the  sample. 
(1)  Union  Carbide  pyrolitic  graphite  foil,  o.Ol  inches  thick 
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Samples  were  deformed  at  a constant  true  strain  rate  of 
0.034  + .003  min  ^ for  a total  of  15  min.  Ram  displacement 
and  load  readings  were  taken  at  one  minute  intervals  and 
converted  to  true  strain  and  true  stress,  respectively, 
using  previously  reported  procedures,  Yan  et  al.  (1975a, b)  . 
Samples  were  either  hot  ejected  (v3  minutes)  or  the  sample 
and  die  were  allowed  to  cool  to  room  temperature  (5-10  hrs . ) . 

Samples  were  mounted  and  polished  to  a ly  finish  using 
diamond  paste.  This  was  followed  by  .3  and  .06y  alumina. 

A thirty  minute  hand  lapping  on  the  .OSy  alumina  was  found 
to  greatly  improve  the  surface  condition. 

Several  etches  were  tried  including  concentrated 
U^SO^,  and  HCIO^^  and  hot  aqueous  solutions  of  NH^Cl.  The 
best  etching  behavior  was  obtained  from  a solution  of  50  ml 
HCl  + 20g  NH^Cl  + 20g  AlCl^  in  100  ml  H2O  at  70°C  or  from 
hot  or  cold  concentrated  HCl.  Most  samples  were  etched  in 
concentrated  (fresh)  hydrochloric  acid  at  room  temperature 
for  20-30  minutes.  Microstructures  were  examined  and  the 
mean  linear  intercept  determined. 

E.4.  Results: 

E.4.1.  Stress  Strain  Behavior 

Stress  strain  curves  were  constructed  for  each  sample. 
Upon  comparison  several  trends  were  observed.  First,  after 
some  amount  of  strain,  the  stress  approached  a constant 
(steady  state)  value  varying  between  several  hundred  to 
over  eleven  thousand  psi  depending  upon  the  temperature  of 
the  run  and  orientation  of  the  crystal.  Both  the  amount  of 
strain  to  reach  steady  state  and  the  magnitude  of  the 
steady  state  flow  stress  were  found  to  decrease  with 
increasing  temperature.  At  temperatures  close  to  0.5T  a 
Strain  of  the  order  of  40%  was  necessary  to  achieve  a constant 
stress  condition.  At  T = 0.8T„,  a strain  of  20%  was  sufficient, 
the  necessary  strain  continuing  to  decrease  as  T increased. 
Figure  El  i] lustrates  these  trends. 
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TRUE  STRAIN 

Figure  El.  Representative  stress  strain 


curves  for  forged  CaF 


For  T >0.65T^the  variation  in  steady  state  flow  stress 

due  to  orientation  was  too  small  to  be  observable.  The  flow 

stress  behavior  of  undoped  oriented  and  unoriented  crystals 

above  0 . 65T^  could  easily  be  described  by  a single  line, 

Figure  E2.  At  T < 0.6  5T^^the  steady  state  flow  stress 

behavior  was  orientation  dependent.  For  any  temperature 

within  this  regime  it  was  found  that  o >a  >n 

„ . <100>  <110>  <111>' 
Extrapolation  would  suggest  that  the  flow  stress  behavior 

for  the  various  orientations  is  converging  to  a single  line 

at  T = 0.45  Tjn.  Experimental  verification  of  this  is 

considered  unlikely  at  the  strain  rates  used  in  our 

experiments.  It  was  found  on  several  occasions  that  samples 

forged  at  T < 0 . 5 shattered  or  cracked  during  forging. 

Prior  to  the  forging  of  oriented  crystals  a certain  amount 

of  scatter  in  forging  stress  had  been  observed  on  un- 

oriented  crystals  at  lower  temperatures,  Bowen  et  al(1976).  If  one 

anticipates  the  values  of  a for  unoriented  crystals  to  fall  somewhere 

between  the  two  extremes  of  the  oriented  crystals,  the 

"scatter"  is  easily  explained. 

A limited  amount  of  work  was  done  on  doped  (100  ppm  Gd) 
oriented  <110>  crystals.  The  flow  stress  was  found  to 


deviate  rather  considerably,  but  erratically  from  t values 
for  undoped  <110^  crystals.  At  some  temperatures  a was 
increased  by  as  much  as  a factor  of  3 at  others  there  was 
a much  smaller  effect.  All  doped  samples  tested  had  higher 
flow  stresses  than  equivalent  undoped  crystals  indicating 
that  there  was  some  amount  of  dopant  in  all  crystals.  Our 
supposition  is  that  the  inconsistency  in  the  flow  stress 
behavior  was  to  a large  part  controlled  by  an  inhomogeneous 
dopant  concentration  from  sample  to  sample. 


E.4.2.  Subgrains: 

Subgrain  formation  occurred  during  the  forging  of  all 
samples.  The  subgrain  size  decreased  with  increasing  stress 
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(decreasing  temperature) . The  results  of  TEM  work  (see 
section  F)  indicate  that  the  subgrain  boundaries  are  low 
angle.  The  subgrain  boundaries  appeared  to  become  more 
diffuse  as  the  temperature  decreased  (o  increased) . The 
combination  of  increasing  boundary  width  and  the  decreasing 
subgrain  size  made  optical  microscopy  of  some  samples 
extremely  difficult. 

In  the  course  of  our  investigation  many  samples  were 
polished,  etched,  and  examined.  In  our  work  on  unoriented 
crystals  there  seemed  to  exist  a sample  to  sample  variation 
in  the  details  of  the  microstructure f Bowen  et  al  (1976,1977). 

The  work  on  oriented  single  crystals  suggests  that  the  uniformity  of 
the  ir.icrostructure,  the  shape  of  the  subgrains  developed,  and  the 
susceptibility  to  cracking  are  orientation  dependent. 

Crystals  with  a <100>  orientation  were  forged  at  temp- 
eratures ranging  from  0.4  8 to  0.74  T^^.  All  samples  were 
removed  from  the  die  free  of  optically  detectable  cracks. 

Only  the  sample  forged  at  0.48  T^  cracked  during  subsequent 
polishing  and  handling.  The  fracture  was  transgranular  in 
nature.  As  the  forging  temperature  was  decreased  the 
sample  shape  tended  to  be  more  blocky.  Samples  forged  at 
0.74  (940°C)  , 0.G6  (800°C)  and0.60T^^  (700°C)  displayed 
equiaxed  grains  which  were  uniform  both  in  size  and  dis- 
tribution, Figure  E3.  There  was  some  evidence  of  recrystal- 
lization in  the  sample  forged  atu,-74Tj^.  Occasionally  a 
subgrain  was  observed  which  was  perhaps  a factor  of  ten 
larger  than  the  average.  The  interior  of  these  larger 
subgrains  were  free  of  etch  pits  whereas  the  matrix  was  not. 

The  samples  forged  atO.  54  (630°C)  and  0.48  T^^  (510°C)  showed  evidence  of 
localized  strain  and  a greater  variation  in  the  subgrain 
size  about  the  average.  In  the  bandy  regions  some  sub- 
grains appeared  to  be  nearly  rectangular.  Parallel  bands 
of  subgrains  with  subgrain  free  regions  between  them  were 
observed  in  both  samples.  The  fraction  of  the  sample  which 
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Figure  E4.  Inhomogeneous  microstructure  in  <11 
sample  forged  at  0.56  Tm  (640°C) . 
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was  subgrain-free  increased  as  the  temperature  decreased. 

Forgings  of  <11]>  crystals  were  conducted  at  tempera- 
tures ranging  fromO.49  (530°C)  to  0 . 66  ( 800  °C ' All  but 

the  sample  forged  at  0.66  Tm  underwent  fracture  either 

during  polishing  (0 . 56  aud  0 . 6lTm)  or  during  forging  (0.49  T ). 

m 

Crystals  of  <111>  orientation  v.'ere  found  to  be  the  most 
susceptible  to  cracking. 

The  microstructures  of  samples  forged  at  0.66  (800°C), 

0.61  (720“C)  and  0.56  (640°C)  Tj^  were  examined.  The  micro- 
structure in  each  case  was  inhomogeneous.  There  were 
limited  regions  in  which  subgrains  of  uniform  size  were 
observed.  Figure  E4.  Withiii  these  regions  there  was  a 
tendency  for  the  grains  to  be  elongated.  The  most 
prominent  microstructural  features  were  bands  within  which 
the  structure  was  poorly  defined  and  kink  bands.  As  the 
forging  temperature  was  aecreased  the  grains  tended  to 
become  more  elongated,  the  extent  of  "banding"  and  the  areal 
fraction  with  poorly  developed  subgrains  increased. 

Forgings  of  <110>  crystals  were  conducted  at  tempera- 
tures ranging  from  0.51  to  0.78  T . All  forgings  were  found 

m 

to  bQ  crack  free  upon  removal  from  the  die.  Microstructures 
were  not  examined. 

E.4.3.  Subgrain-Subcell-Flow  Stress  Correlation 

Generally  in  the  creep  literature  the  subgrain  size  is 
correlated  to  the  flow  stress.  In  many  metal  systems  it  is 
observed  that  the  subgrain  size  is  inversely  proportional 
to  o.  A similar  correlation  was  sought  in  the  CaF2  system. 

During  the  earlier  portion  of  our  program  the  majority 
of  our  work  was  on  unoriented  crystals.  It  was  observed 
that  the  "subgrain"  size  was  inversely  proportional  to  o^, 
Bowen,  et  al.  (1976  ), Cannon,  et  al.  (1976).  There  was  also 
an  unsatisfying  amount  of  scatter  in  the  data.  We  believe 
that  our  work  on  oriented  crystals  and  the  results  of  TEM 
work  offer  a plausible  explanation. 
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siz.  h"  ““  correlation  between  subgrain 

and  o rt  is  desirable  to  have  a uniforn,  microstructure 
When  confronted  with  inhomogeneous  miorostruc-uras  such  as  ' 
those  Observed  in  <111>  .nystals  of  CaP,  it  ia  not  clear 
the  average-  subgrain  size,  5,  should  be  defined,  if 
one  includes  the  Inhomogeneities  the  value  of  6 becomes 

the  inhomogeneities  and  correlates  5 to  o one  must  assume 

that  the  homogeneity  of  the  sample  does  not  significantly 

influence  the  magnitude  of  These  nrnhi 

inese  problems  are  made 

worse  If  appreciable  dynamic  recrystallization  occurs,  m 

ght  of  these  considerations  it  was  felt  that  the  most 

meaning  ul  correlation  would  be  obtained  from  measurements 
on  <100>  crystals. 

When  etch  pits  were  observed  they  were  generally 

aL'as"o"f  th  subgrains.  Although  in  some 

he  samples  the  arrangement  of  etch  pits  appeared 
o e random,  m a substantial  fraction  of  the  area  the  etch 
Pi  s ormed  cellular  networks.  Figure  E5.  There  were  two 
distinct  microstructUral  elements  present  rn  the  sampi;;. 

ce  u ar  units,  whose  boundaries  were  defined  by  dis- 
Crete  etch  pits  are  referred  to  m this  paper  as  subcells. 

by  d Th  ”"b  (=hbcell  Size)  designated 

Pits  and  hence  appeared  to  be  continuous  are  referred  to  as 
subgrain  boundaries.  The  subgrain  size  is  designated  by  5. 

to  the  3 to  be  inversely  proportional 

he  stress  the  magnitude  of  6 ranging  from  10  to  60m  The 

subcell  size  was  found  to  vary  between  4 and  28m.  In  t^e 

course  of  TEM  worh  on  our  unoriented  specimens.  Sherry  and 

in  l"d  f section  F measured  cell  sizes.  Their  data  is 
included  and  complements  our  measurements. 

in  several  samples  with  a <111>  orientation  subgrains 
and  suboells  were  observed.  The  subcells  were  not  as  clearly 
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defined  as  in  ^100>  forgings  and  an  accurate  measurement  of 
d was  not  obtainable.  The  subgrain  size  in  <111>  forgings 
was  smaller  than  those  in  <100>  forgings  of  comparable 
stress.  This  may  be  the  result  of  the  mj crostruotural 
inliomogeiei ty  observed  in  ■ 111>  forgings. 

The  subcell  size  and  subgrain  size  determine  a band  on 
a log  5(d)  vs.  log  a plot.  Figure  E6.  Nearly  all  our  earlier 
data  falls  within  this  band.  Cur  data  most  likely  represents 
a mixture  of  subcell  and  subgrain  size.  This  may  be  due  to 
ouv  use  of  an  etchant  which  did  not  preferentially  etch 
dislocations . 

Ihe  reasons  for  the  difference  in  subcell  and  subgrain 
etching  behavior  are  not  clear.  We  assume  that  there  must 
be  an  energy  difference  which  accounts  for  the  difference. 

The  subgrain  boundaries  perhaps  correspond  to  a higher  mis- 
orientation  angle  and  hence  a smaller  disloccition  spacing  in 
the  boundary.  TEM  measurements  on  misorientation  angles  by 
means  of  spot  splitting  have  shown  that  0 can  vary  substan- 
tially within  a gi\'’en  sample.  Within  some  samples  a nearly 
bimodal  distribution  was  observed.  Table  I.  A problem 
arises  with  this  interpretation  when  one  compares  the  TEM 
values  of  misorientation  angle  for  subcells  (presumably  the 
m0.3°  values)  with  those  determined  from  etch  pit  spacing. 

Based  on  etch  pit  spacing  one  predicts  a misorientation  of 
the  order  of  0.01°.  The  possibility  exists  that  the  pits 
correspond  to  dislocations  which  have  not  yet  been  accomodated 
into  the  boundary.  These  extrinsic  dislocations  may  have  a 
higher  energy  than  the  knittec!  boundary  dislocations.  As  the 
energy  of  the  boundary  increases  with  increasing  0,  subgrain 
boundaries  may  etch  without  the  presence  of  extrinsic  dis- 
locations. Alternatively  the  etch'’’i7  of  individual  dislocations  may 
depend  critically  on  th«  angl*  of  intersection  of  the  dis- 
location with  the  surface. 
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table  I 


Sub-boundary  Misorientation  Angle 


Temperature 

590“C 


Misorientation  Ang^l^^  Degrees 


0.374 

1.206 

0.695 

1.273 

0.  781 

1.302 

0.955 

1.507 

1.042 

1.507 

1.074 

Mean 


1.07 


702°c 


817°C 


907°C 


0.318 

0.465 

0.478 

0.786 

0.310 

0.337 

0.674 


0.419 

0.441 

0.559 

0.637 

0.800 

0.819 


1.145 

1.309 

3.55 

4.010 

1.528 

1.637 

1.763 

2.191 

1.055 
1.066 
1.146 
1.273 
1.302 
1.  375 
1.432 


1.042 


1.563 


0.955 
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E.4.4.  Recrystallization 

The  microstructure  of  samples  forged  at  temperatur-^s 
in  excess  of  0.6  5 showed  increasing  evidence  of  recrystal- 
lization as  the  temperature  was  increased.  At  temperatures 
close  t:)  0.65  Tm  the  evidence  was  limited.  One  occasionally 
observed  larger  than  average  subgrains  outlined  by  a boundary 
which  appeared  to  have  a higher  than  average  misorientation 
angle.  These  "nuclei"  became  more  frequent  as  the  forging 
temperature  increased.  This  behavior  is  illustrated  in 
Fig.  E7.  The  presence  of  poorly  defined  subgrain  boundaries 
within  the  large  grain  suggests  that  recrystallization 
occurred  during  forging.  At  T ^0.79Tj^  regions  with  well 
defined  subgrains  became  increasingly  sparse.  These  regions 
were  limited  to  the  edges  of  the  forged  samples,  Fig.ES.  It 
is  thought  that  the  tendency  at  higher  temperatures  to  reach 
apparent  steady  state  at  lower  strains  may  be  related  to  the 
considerable  amount  of  dynamic  recrystallization.  The  few 
samples  which  were  Gd"*"^  doped  did  not  appear  qualitatively 
different  in  the  tendency  for  dynamic  recrystallization. 

'I’^O.STm  samples  became  opaque  after  forging.  This 
is  thought  to  be  the  result  of  CaO  formation  on  cooling  as  a 
result  of  0 contamination.  As  a result  of  the  extensive 
recrystallization  at  high  temperature  no  attempt  was  made  to 
eliminate  this  problem. 

E.4.5.  Thermal  Shock 

One  of  the  problems  which  is  likely  to  be  encountered 

in  materials  which  have  poor  thermal  conductivity  coupled 

with  large  thermal  contraction  is  quench  cracking.  Unoriented 

samples  which  were  forged  and  hot  ejected  tended  to  crack. 

There  appeared  to  be  a temperature  range  ('v0.6BT  -'\^0.81T 

mm 

over  which  samples  did  not  crack  due  to  hot  ejection.  In 
light  of  the  variability  between  <100>  and  <111>  in  their 
tendency  to  crack  during  polishing  we  suspect  that  this  may 
have  been  caused  by  orientation  variations. 
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Recrystallizod  grains  in  CaF2  sample 
forged  at  0.79  Tj„  (1018°C)  . 


Sample  forged  at  0.79  Tm  (1018°C)  showing 
well  developed  subgrains  only  near  the 
sample  surface. 
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F.  Transmission  Electron  Microscopy  of  CaF2 
F.l.  Introduction 

An  investigation  to  characterize  the  microstructure 
developed  during  hot  forging  of  single  crystal  CaF2  has  been 
accomplished  using  transmission  electron  microscopy.  Single 
crystal  CaF2  samples  hot  forged  at  temperatures  of  590,  702, 
817,  and  907°C  were  examined  in  the  transmission  electron 
microscope  using  a double~tilt ,liguid~helium  cold  stage. 

F.2.  Experimental  Procedure 

Hot  forged  CaF2  specimens  were  cut  to  approximately  500  y 
thick  slices  using  a string  saw.  The  cutting  solution  is 
comprised  of  SiC  powder,  glycerin  and  water,  and  the  string 
saw  wire  is  375  y in  diameter.  The  500  y thxck  specimens 
were  then  mounted  on  a "polishing  jig"  using  a low  melting 
point  wax  and  polished  on  500  grit  emery  paper.  After 
polishing,  the  specimens  were  approximately  375  y in 
thickness.  The  slices  of  CaF2  were  next  cut,  either  using 
the  string  saw  or  a surgical  blade,  to  a size  that  will  fit 
into  a 2.3mm  specimen  holder.  These  "2.3  mm  diameter"  by  375  y 
thick  specimens  were  then  chemically  jet  thinned  for  tranS” 
mission  electron  microscopy. 

Chemically  jet  thinning  CaF2  for  transmission  electron 
microscopy  is  a difficult  process.  This  arises  for  two 
reasons.  First,  because  CaF2  is  transparent,  conventional 
means  for  detecting  the  appearance  of  a hole  in  the  specimen 
during  the  jet  thinning  process  cannot  be  employed?  and 
second,  because  CaF2  polishes  (etches)  nonuniformly . 

Several  solutions  have  been  tried  as  possible  candidates 
for  the  jet  chemical  thinning  process?"  The  chemical  solution 
which  has  yielded  the  most  uniform  polish  of  any  tried  to  date 
is  given  below. 
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A 


Chemical  Polish’ ng  Solution 
at  70°C 


400  ml  H^O 
100  ml  HCl 
70  gm  AlCl^ 

70  gm  NH^Cl 

Thin  foils  are  made  using  a dual  jet  thinner.  The 
specimen  is  observed  during  the  thinning  process  with  a 
light  microscope.  Transmitted  light  is  used  to  illuminate 
the  specimen  surface.  When  the  light  is  properly  aligned, 
and  "reflected"  light  is  used  to  observe  the  specimen, 
fine"  structure  on  the  surface  of  the  specimen  may  be 
observed.  When  a hole  appears  it  is  observable;  but  gener- 
ally not  until  it  has  become  0.5  to  1 mm  in  diameter.  Due 
to  the  size  of  the  hole,  a major  portion  of  the  thin  area 
is  lost,  making  the  transmission  electron  microscopy  more 
dif  f ici; . c . 

Before  a specimen  is  placed  in  the  electron  microscope, 
it  is  coated  with  a thin  layer  of  carbon.  The  carbon  is 
applied  by  vapor  deposition  and  prevents  specimen  charging. 

TEM  of  CaF2  3t  room  temperature  is  precluded  by  interference 
from  radiation  damage  from  the  electron  beam^.  Thus  TEM  observa- 
tions were  made  below  60°K  using  a liquid-helium  cold  stage  dev- 
eloped under  the  previous  contract.^ 

F.3.  Results  and  Discussion 

The  results  obtained  from  the  observation  of  the 
deformation  structures  of  hot  forged  CaF2  will  be  diviaed 
between  an  analysi s of  the  sub-grain  boundaries  and  the 
J^elation  of  the  sub-grains  to  the  observed  macroscopic 
behavior . 

F.3.1.  The  Sub-Grain  Boundary  Structure 

Observations  of  the  microstructure  resulting  from  hot 
forging  has  led  to  the  characterization  of  two  types  of 
sub-grain  boundaries  (SGB’s).  The  types  of  SGB's  developed 
during  the  forging  process  are  dependent  upon  forging 


temperature  (flow  stress)  and  are  classified  as  equilibrium 
or  nonequilibrium  boundaries.  The  characterization  of  the 
two  types  of  boundaries  has  been  done  by  considering  the 
defect  comprising  the  SGB  structure.  These  are; 

(1)  Intrinsic  sub-grain  boundary  dislocations  which  are  a 

requirement  of  the  boundary  geometry  and  are  associated 
with  equilibrium  SGB's. 

(2)  Extrinsic  sub-grain  boundary  dislocations  which  are 
those  dislocations  resulting  from  the  accommodation  of 
matrix  dislocations  into  the  boundary  and  are  not  a 
requirement  of  boundary  geometry  but  may  decrease  the 

energy  of  the  system  while  increasing  the  energy  of  the 
boundary . 

Figure  F.l.  presents  a series  of  micrographs  typical  of 
the  sub-gram  boundaries  observed  after  forging  at  590 °C. 

Some  sub-grain  boundaries  are  highly  disordered  and  appear 
to  have  a nonequilibrium  structure.  Figure  F.la,  b. 

Evidence  of  the  interaction  of  extrinsic  dislocations  with 
those  of  the  boundary  is  seen  in  Figure  F.la  and  b.  At  A 
extrinsic  (matrix)  dislocations  are  interacting  with  the 
boundary  but  are  not  accommodated  into  the  boundary.  These 
disordered  boundaries  should  have  stress  fields  of  signifi- 
cant proportion  and  may  affect  the  deformation  behavior.  In 
addition  to  the  highly  disordered  sub-grain  boundaries, 
several  cwist,  Figure  F.lc,  and  tilt  boundaries  were  observed 
however,  the  majority  of  the  observed  boundaries  were  highly 
disordered.  The  twist  boundary  of  Figure  F.lc  consists  of 
intrinsic  dislocations  and  is  highly  equilibrated.  Further 
evidence  of  a nonequilibrium  microstructure  was  manifest  by 
a high  matrix  dislocation  density.  The  matrix  dislocations 

generally  existed  in  gross  tangles  and  were  extensive  in 
nature . 

The  effect  of  increasing  the  forging  temperature  from 
590°C  to  702°C  may  be  seen  in  Figure  F.2.  The  boundaries 
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sub-grain  bound 
rrom  hot  forging  at  590»C:  (a)  a - 

5 = (022),  (c)  g = ,ui).  ' ® 


TEM  micorgraphs  of  sub-grain  boundaries  produced 
from  hot  forging  at  702‘’C;  (a~d)  g = (022)  . 


are  „K,re  ordered.  Both  tilt.  Figure  F.2a,  and  t„ist, 

Frgura  F.2o,  houndarias  ara  aean  as  „aU  as  tha  „ora  complex 
boundarras.  The  matrix  dislocation  density  is  lower  than 
that  at  590«c.  Tha  ™ora  ordered  boundariel.  as  T rgua  f 2a 
and  F.2C  should  have  short-range  stress  Fields  and  o d n" 

tha  . . presents  a simple  tilt  boundary  whe -a 

the  rntrrnsrc  dislocations  are  easily  discerned.  The 

presence  of  extrinsic  dislocations  in  the  boundary  can  be 

en  drrectly  and  also  by  the  discontinuities  in  the  intrinsic 

tinuities°"'  drslocations  causing  the  discon- 

may  be  seen  by  using  other  operating  reflections) 
ese  extrinsic  dislocations  have  been  accommodated  into  the 
boundary.  In  an  effort  to  determine  the  misorientation 

tween  the  two  grains  in  Figure  F.2a,  an  electron  diffraction 

experiment  measuring  the  shift  in  Kibuchi  lines  between  1 
grains  was  performed  Th«  mr.  tween  two 

tion  Of  0 36.  A ■ yibWed  a misorienta- 

■ of  0.36  . Assuming  that  the  boundary  in  Figure  F 2a  is 

agree'en  itTth:: 

the  n measurements.  The  nature  of 

he  Burgers  vectors  of  the  dislocations  comprising  the  sub- 

determined. 

The  reactions  which  occur  are; 

(1)  A + B = c 

[100) 


(2)  E + F = D 


7 1101]  + 1 tlOll 


2 11011  + i toil)  * I [110] 


Cherl'b  r"tr  io  Proportional  to 

- ^ Bor^ors  vector  of  the  dislocation.  The 

energy  requirements  of  (1)  and 

4-ho  u2  . . ^ ° calculated  usina 

:::i:  ettoT^rir^iir^^^ 

the  energy  of  the  reactants  and  products'of  ur^re 'eglT 
C oaer  attention  must  be  given  to  the  exact  nature  orthe 
dislocations.  The  Burgers  vector  of  a and  B are  inclined  by 
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approximately  45-  and  30-  to  their  dislocation  U„es 
respectiyely.  while  that  of  c is  nearly  parallel.  ;ns 

is  almosr^  “hile'c 

dislooa- 

(«  b2  a^  dislocation 

. y reaction  where  the  reactants  have  a significant 
edge  component  and  the  orodn-h  ^ ^icant 

favoiable.  This  is  tL  1 will  he 

reaction  ' e <2)  - and  the 

1 gureF2d  is  a boundary  which  lies 

icTor:;?:: ----  - -“ion 

(3)  G + H = I i n 1 It 

2 [101]  + - [loi]  = [001] 

Applying  the  b=  criteria  shows  that  the  elasHc  i • 

r>  •F  +-K^  S-l-SStic  StjTfiin  SilGycix/ 

reactants  and  product  dislocations  are  equal  Th 
Burgers  vector  of  r i o • i ■ equal.  The 

::iroTthV:ea“^ 

• action,  I,  ij  inclined  45°  to  thc^ 

line.  Since  the  product  dislocation  is  mixed  th! 
each  dislocation  must  be  calculated  r ■ ' 

to  determine  if  the  reaction  is 


E = 

47T 


^~l"-v^'  cos^e]  In  ^ 


(F.l) 


where  p i.s  the  shear  modulus,  a a constant,  R the  radia  di  st 
from  the  dislocation  coro  fl  a radia  u distance 

line  and  the  Burgers  vector  a d """"""""  dislocation 

urgers  vector,  and  v Poisson's  ratio. 

^products”''  however,  only  by  the  ratio  of 

E reactants  1 • fact  that  the  reaction  is  not  overly 

favorable  may  be  evidenced  by  the  observed  constricted  length 
of  the  product  dislocation. 

The  SGB  s developed  during  hot  forging  at  Sll^c  are 
presented  in  Figure  F.3.  These  boundaries  are  highly  oidered 
and  have  become  equilibrated.  The  matrix  dislocation  d -nsity 
IS  lower  than  that  at  590-C  and  702-c.  The  boundaries  should 
have  short-range  stress  fields.  Figure  F.3a  is  a compound 
boundary  made  up  of  twist  and  tilt  components:  F.3b  shows  the 
intersection  of  two  highly  ordered  twist  boundaries.  The 
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roqraphs  of  sub-q 
f forging  at  817'' 


upper  boundary  has  extrinsic  dislocations  which  have  inter- 
sected the  boundary  although  there  does  not  appear  to  be  any 
interaction  with  the  boundary.  Figure  F.3c  and  F.3d  also 
show  hiohly  ordered  subboundaries  produced  from  the  hot 
forging . 

The  SGB ' s observed  after  hot  forging  at  907°C  are 
presented  in  Figure  F.4.  Nearly  all  the  boundaries  observed 
at  907°C  consisted  entirely  of  intrinsic  dislocations  and 
were  highly  ordered.  The  density  of  extrinsic  dislocations 
was  very  low,  indicating  that  these  dislocations  were  either 
annihilated  by  interactions  within  the  grains  or  accommodated 
in  the  SBG's.  The  SGB  in  Figure  F.4b  is  composed  of  a highly 
ordered  structure  with  a number  of  boundary  defects,  labeled  A. 
These  defects  may  be  either  extrinsic  dislocations  or  grain 
boundary  ledges.  The  misorientation  of  the  SGB  in  Figure  F.4b 
was  found  to  be  1°  by  measuring  spot  splitting  in  the  electron 
diffraction  pattern.  The  apparent  migration  of  extrinsic 
dislocations  to  a SGB  is  shown  in  Figure  F.4d.  These  disloca- 
tions are  being  accommodated  into  the  boundary  at  the  expense 
of  increasing  both  the  angle  and  energy  of  the  boundary  while 
decreasing  the  energy  of  the  system.  - The  accommodation  of 
these  extrinsic  dislocations  into  the  boundary  may  be  rate  a 
controlling  mechanism  in  the  overall  recrystallization  process; 
but  more  importantly  this  may  be  a mechanism  by  which  a "high" 
angle  grain  boundary  is  created.  This  boundary  could  then 
produce  abnormal  or  secondary  grain  growth. 

The  effect  of  increasing  the  forging  temperature  from 
590°C  to  907°C  on  the  observed  microstructure  may  be  summari:;ed 
as  follows: 

(1)  The  matrix  dislocation  density  decreases. 

(2)  The  SGB's  become  increasingly  ordered. 

(3)  There  appears  to  be  increased  interaction  of  extrinsic 

dislocations  with  the  boundary. 
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F.3.2.  The  Relation  of  Sub-Grain  Structure  to  Observed 
Macroscopic  Behavior 

The  previous  sections  have  considered  the  effect  of  flow 
stress  or  forging  temperature  on  the  type  of  sub-grain  boundary 
developed  during  hot  forging  of  single  crystal  calcium  fluoride 
to  compressive  strains  of  approximately  50%.  In  this  section, 
the  results  of  sub-grain  size  observations  and  calculations 
made  with  the  aid  of  transmission  electron  microscopy  are 
presented. 

Micrographs  typical  of  those  used  in  determining  sub-grain 
size  are  presented  in  Figures  F.5  and  F.6.  Figure  F.5  shows 
a variety  of  sub-grain  boundary  structures.  There  appears  to 
be  a duplex  sub-grain  boundary  structure  present.  Within  a 
more  well  defined  and  ordered  subboundary  are  several  lower 
dislocation  density,  less  ordered  subboundaries.  The  larger, 
more  ordered  subboundary  has  dimensions  of  20-30  pm  in 
diameter,  whereas  the  less  ordered  subboundaries  are  'v  4-6  pm 
in  size.  Although  Figure  F.5  clearly  suggests  a duplex  sub- 
grain boundary  structure  it  is  not  typical  of  that  which  is 
generally  observed.  Figure  F.6  shows  sub-grain  boundary 
structures  more  typical  of  the  .subboundaries  observed  as  a 
result  of  forging  at  temperatures  590,  702,  795,  and  907°C. 
The  boundary  structure  is  much  more  uniform  although  there  is 
still  a suggestion  of  a duplex  sub-grain  boundary  structure. 
However,  any  suggestion  of  a duplex  structure  disappears  at 
the  highest  forging  temperatures. 

Table  F.l  below  presents  a compilation  of  the  sub-grain 
boundary  calculations  made  on  specimens  hot  forged  at  590, 

702,  795,  and  907°C. 
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Figure  F.6.  Transmission  electron  micrograph  of  the 

microstructure  developed  from  forging  at  590 °C. 


Table  F.l 


(°C)  d (ym)  g (psi) 


590 

5.41; 

3.69 

5300 

702 

4.79; 

5.54 

4000 

795 

6.63; 

5.54 

3300 

907 

6.82 

2560 

A plot  (Figure  F.7  of  In  a (flow  stress)  versus  In  d 
(sub-grain  size)  demonstrates  the  power  law  dependence  of 
flow  stress  on  sub-grain  size.  A least  squares  analysis  of 
the  points  plotted  in  Figure  F.7  yields  the  linear  In-ln 
relationship  presented  in  Equation  (F.2) 

In  a = -1.027  In  d + 9.978  (f.2) 

Equation  (F.2)  can  be  converted  to  the  power  law  dependence 
of  Equation  (F.3) . 


o 3:  2.155  X lo"^ 


^-1.027 


(F.3) 


Compilation  of  available  data  in  metallic  materials^,  e.g.  Cu, 
Al,  Fe  and  Fe  + 0.007%  C,  demonstrates  that  the  flow  stress  is 
inversely  proportional  to  the  subgrain  or  cell  size.  Similar 
results  are  seen  in  KCl^,  and  KBr^. 
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Figure  F.7. 


Plot  of  In  (flow  stress)  as  a function  of  In 
(sub-grain  size) . 
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PRINCIPAL  STRESS  EFFECTS  ON  BRITTLE  CRACK  STATISTICS 


G.l  INTRODUCTION 

Ideally,  the  initiation  of  a micro-crack  will  not 
imiucaiately  lead  to  macroscopic  cracking,  even  in  brittle 
material.  For  this  to  be  true,  the  first  micro-crack  must 
be  stabilized  by  friction,  plastic  flow,  or  inhomogeneities. 
Only  after  there  are  a number  of  micro-cracks  will  they 
localize  to  form  a raacrocrack  and  limit  the  load  to  a 
i^sximum.  Calculating  the  progre*"”  of  micro— cracking 
towards  the  configuration  giving  maximum  load  then  involves 
repeated  stress  analysis  of  the  developing  micro-crack 
distribution,  as  well  as  the  statistics  of  the  strengths 
of  the  elements  that  are  about  to  break.  This  problem 
has  been  treated  numerically  and  analytically  by  McClintock 
and  Zaverl  (1974),  with  the  following  general  assumptions: 

a)  A two-dimensional  array  of  regular  hexagonal 
grains . 

b)  Grain  boundary  strengths  against  micro-crack 
nucleation  or  propagation  from  an  adjoining 
grain  that  vary  according  to  an  extreme-value 
distribution. 

c)  Homogeneous  elastic  constants. 

d)  No  dynamic  effects. 
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e)  No  plastic  flow. 

f)  Biaxial  applied  tension. 

The  last  restriction  is  relaxed  here  so  as  to  find  the 
effects  of  changing  the  ratio  of  the  applied  tensile  stress 
components,  while  still  maintaining  enough  tension  to 
avoid  the  effects  of  crack  closure  and  friction.  While 
these  models  are  still  far  from  representing  actual 
brittle  materials  such  as  rocks  and  ceramics,  it  is  hoped 
that  they  provide  a worth-while  improvement  in  insight 
over  previous  approximations  that  consider  the  effects 
of  statistics  without  the  micro-stress  analysis,  or  stress 
analysis  of  regular  arrays  of  micro-cracks. 


G.2  STATISTICAL  DISTRIBUTION  OF  GRAIN 
BOUNDARY  STRENGTH 


The  grain  boundary  strengths  are  taken  to  follow 
an  extreme-value  distribution  of  the  second  kind,  with 
the  distribution  of  strengths  increasing  as  some  power  m 
of  the  amount  by  which  the  local  normal  stress  on  the 
center  of  the  grain  boundary  exce  5 so;  e lower  limit 
of  strength  . More  specifical  . , uhe  probability  that 
the  strength  S of  a grain  boundary  lies  below  is 

given  in  terms  of  the  parameters  S^^,  S^,  and  m by 


p(S<S^)  = 1 


exp 


4s,-s  ^ 

m 

.1  L 

1 s„-s. 

_ \ 0 

(G.  1) 
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Then  fracture  nucleates  at  the  center  of  a grain  segment 
when  the  normal  stress  reaches  a critical  value 

*^nn  ^1  * (G.2) 

For  a review  of  the  derivation,  and  tables  of  mean  and 
standard  deviation  in  terms  of  part  size  and  the  parameters, 
see  for  example  McClintock  and  Zaverl  (1977).  In  general, 
one  might  expect  the  shear  as  well  as  the  normal  component 
of  stress  to  affect  fracture  of  a grain  boundary,  but  any 
such  shear  effects  will  be  neglected  here. 

In  the  presence  of  an  adjacent  crack,  a grain 
boundary  segment  may  still  crack  by  nucleation,  as  described 
above,  or  it  may  crack  by  propagation.  In  this  latter  case, 
the  resistance  to  cracking  will  be  described  by  the  singu- 
larity of  normal  stress  along  the  grain  boundary.  This 
singularity  of  normal  stress  is  related  to  the  stress 
intensity  factor  at  the  crack  tip  and  the  orientation  0 
of  the  grain  boundary  relative  to  the  crack.  To  eliminate 
the  repetitious  rr  we  follow  Sih  and  Liebowitz  (1968) 
and  define  the  stress  intensity  in  terms  of  crack  half- 

length  c and  the  applied  normal  stress  at  infinity  S 
such  that 

^ )e=0 
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(The  more  familiar  K_  is 

X 


Kj  = /tT  .)  (G.4  ) 

Then  the  intensity  of  normal  stress  on  a boundary  at  an 
angle  0 to  the  crack  from-  Sih  and  Liebowitz  (1968)  is 

kg  = Cgg  y 2r  = k^  (1/2)  (1  + cos  6)  (cos  6/2)  . (3-5  ) 

The  critical  value  of  this  normal  boundary  intensity  factor 
is  again  assumed  to  have  an  extreme-value  distribution  of 
the  second  kind,  so  that  the  probability  of  its  value 
being  less  than  some  particular  value  k.,  is 


P = 1 - exp 


|^91~^eL 
^0O“’^0L 


A' 


(G  .6) 


For  insight,  the  statistical  parameters  ko^  and  k-_ 

UIj  uu 

may  he  normal  i in  terms  of  the  grain  boundary  half 
length  c^  , conversion  factors  ^)^0L  ^k0O  ' 

and  the  corresponding  parameters  for  crack  nucleation; 


^9L  ~ ^k9L  ^l/^1  ' ^90  ^k0O  ’ ( G.  7) 


Then  if  propagation  is  viewed  as  nucleation  by  the  normal 
stress  S at  a radial  distance  6/c^  from  the  crack  tip, 

“ ^k6L  , or  , (G.8) 
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and  similarly  for  f 


k0O  * 

G.3  NUMERICAL  MODEL 

A numerical  model  was  used  to  suggest  assumptions 
for  and  to  test  an  analytical  model. 

The  numerical  model  is  a form  of  boundary  integral 
method  in  which  the  traction-free  cracks  are  modeled  by 
dislocation  arrays  simultaneously  chosen  to  cancel  the 
effects  of  the  applied  stress.  For  simplicicy,  each  crack 
is  approximately  modeled  by  one  dislocation  pair,  so  the 
number  of  degrees  of  freedom  is  just  the  current  number 
of  crack  segments.  On.ce  the  necessary  dislocation 
strengths  have  been  found  for  any  given  crack  configuration 
the  stress  is  calculated  at  each  grain  boundary,  and  the 
stress  intensities  are  found  at  each  crack  tip.  The 
weakest  element  relative  to  its  applied  stress  is  then 
found,  assumed  to  be  cracked,  and  the  calculation  is 
repeated.  Fracture  is  determined  by  the  weakest  segment 
relative  to  the  current  stress  on  it,  or  if  a crack  tip 
is  more  critical,  by  the  lowest  segment  toughness  relative 
to  the  current  local  intensity  of  normal  stress  across  a 
grain  boundary.  After  the  most  critical  segment  is  found, 
it  is  assumed  cracked.  Then  a new  elastic  solution  is 
calculated,  a new  most  critical  segment  is  found,  and  so  on 
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Some  typical  results  are  shown  in  Figs.  G-1  and 
G-2.  The  cracked  segments  are  numbered  in  order  of  cracking 
Those  that  cracked  after  a maximum  load  had  been  attained 
are  denoted  by  dashed  lines.  Two  different  orientations 
were  considered  for  the  hexagonal  array  relative  to  the 
direction  x^^  of  the  maximum  principal  component  of 
applied  stress,  . in  the  "vertex"  orientation,  one 

of  the  vertices  of  the  hexagons  is  in  the  direction  of 
maximim  principal  stress  component;  in  the  "face"  orienta- 
tion one  of  the  faces  of  the  hexagon  is  in  that  direction. 

With  the  vertex  orientation  of  Fig.  g-1,  the  cracks 
tend  to  form  along  the  set  of  zig-zag  lines  at  ± 60°  to 
the  maximum  principal  stress  direction.  Under  combined 
stress,  it  is  then  reasonable  to  build  an  analytical  model 
on  this  pattern,  taking  into  account  the  combined  effect 
of  the  two  applied  stress  components  on  the  normal  stress 
required  to  initiate  these  60°  cracks,  and  the  critical 

normal  grain  boundary  intensity,  kg  , required  to  propagate 
them. 


For  the  face  orientation  of  Fig.  G.2,  as  the 
transverse  stress  is  reduced  the  cracks  tend  to  form  along 
the  90°  segments  rather  than  on  those  at  30°  to  the  maximum 
principal  stress  direction.  Here  the  model  must  be  more 
complex.  First,  the  critical  event 


seems  to  be  not  the 


Fig«  G,i  Sample  cracking  pattern  for  a vertex 
orientation. 


Fig.  G.2  Sample  cracking  pattern  for  a face 
orientation. 
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formation  of  the  normal  cracks,  so  much  as  the  subsequent 
cracking  of  the  30'*  cracks.  These  latter  are  thus  taken 
to  be  tne  unit  building  blocks  out  of  which  an  approximate 
statistical  theory  for  macro-cracking  will  be  built. 

Further  discussion  of  the  results  from  the  numerical 
theory  will  be  deferred  until  the  analytical  model  has 
been  developed  for  comparison. 

G.4  ANALYTICAL  MODEL 

The  essence  of  the  analytical  models  is  to  estimate 
the  maximum  strength  by  first  finding  the  density  of  unit 
cracks  as  a function  of  applied  stress  level.  The  prob- 
ability of  aggregations  of  length  2c  of  these  cracks  is 
then  estimated.  When  the  crack  lengths  and  applied  stresses 
are  high  enough  to  cause  crack  extension  even  with  median 
strengths  at  the  crack  tip,  a maximum  strength  is  assumed 
to  have  been  attained. 

Vertex  Orienteition 

For  the  vertex  orientation,  the  cracking  is  observed 
to  be  predominantly  in  the  directions  at  ± 60°  to  the  direc- 
tion of  maximum  tensile  stress.  For  nucleation,  the  normal 
stress  on  these  segments  is  found  by  Mohr's  circle  from  the 
maximum  and  minimum  principal  stress  components  to  be 

'^06  “ ^ ^2^^  ■ (g.9) 
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The  fraction  of  ± 60°  boundaries  that  are  cracked  at 


principal  applied  stress  components  and  S2  is  now 

found  by  substituting  back  into  Eq.  9.1,  and  noting  that 
segments  of  this  orientation  comprise  2/3  of  the  total 
number  of  grain  segments: 


p = 1 - exp 


I - "l 


m 


(G.  10) 


The  total  number  of  cracked  grains  in  a part  of  size  N is 


= (2/3)  N p 


(G.ll) 


With  this  fraction  of  cracked  grains,  we  ask  the 
probability  that  a contiguous  chain  of  has  formed. 

Following  McClintock  (1974),  approximate  the  actual  condi- 
tions by  a linear  array  of  (2/3)  N grain  segments  of  which 
the  fraction  p are  cracked.  The  longest  contiguous  crack 
which  has  a probability  <I>  of  forming  is  approximately 


N = 


'n  ^ (2/3) 

V.  7. 5-l/fcn(p))(-£n  (1-<I))  ) / 

(P) 


(G.12) 


The  strength  of  the  crack  of  length  is  estimated 

from  fracture  mechanics  as  the  applied  stress  required  to 
crack  the  next  segment  ahead,  assuming  no  other  cracks  in 
the  neighborhood  of  the  tip.  (This  is  an  overestimate  if 
any  neighboring  cracks  lie:  ahead  of  the  crack,  but  an 


underestimate  if  any  such  cracks  lie  off  to  the  side  of 
the  tip,  tending  to  shield  it  as  a plastic  zone  would.) 
Furthenr.ore,  assume  the  average  direction  of  the  crack 
IS  normal  to  the  maximum  principal  applied  stress  and 

the  length  is  the  projected  length  of  the  segments, 

2c^cos  30°  . Then 


^Ic  ~ ^*^1  30°)/2 


(G.13) 


For  a long  crack,  k^^  is  found  by  the  lower  of  the  two 
median  strength,  for  nucleation  and  propagation.  Approxi- 
mating the  prior  zig-zag  crack  segments  by  a straight  line 
as  shown  in  Fig.  G,3,  the  normal  grain  boundary  intensity, 
for  propagation,  is  reached  when 

= ^ic  ^ cos  30°)  (cos  15»)/2j  = 0.901  k^^  . (a. 14) 

Fo^  a median  value  of  the  normal  grain  boundary 
intensity  and  noting  that  there  are  two  possible  crack 
sites,  one  at  each  crack  tip,  from  Eqs . G.6  and  G.7,  and 
combining  Eqs.  G.  13  and  G.  14  , 


0.5  = 1 - 


exp 


JO.901  s^/tTco. 


m 


(G.15) 


If  the  next  segment  cracks  by  nucleation,  use  Eq.  G.l  for 
the  distribution  of  strengths  and  find  the  stress  on  the 
next  grain  from  the  stress  intensity  factor: 
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Fig  . G . 3 


Approximating  the  stress  intensity  with 
applied  stress  in  the  vertex  direction. 
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(d  + cos  30°)  (cos  15°)/2]/72E^  = 0 . 901k^^//27^.  (G 


a 


06 


k 


Ic 


Combining  Eqs.  G.l,  G.13,  and  G.16, 


0.5  = 1 


exp 


0 

(0.901  sj.5^  cos  30°  - S ^ 
1C  u 

m 

— z 

S„  - 

V "^0  L 

- 

. (G  .17) 


Since  we  assume  k^  and  Ogg  to  be  calculated  from  the 
same  random  sample,  the  choice  between  crack  growth  by 
propagation  and  re-nucleation , Eqs.  G . 15  and  G , 17 , is  set 

by  the  parameters  ' ^kOO  ' ^0  ’ 

maximum  stress  is  now  the  value  of  (for  a given 

S_/S.)  found  by  eliminating  p and  N from  Eqs.  G.,  10, 
G.12,  and  G.15  or  G,17. 


Face  Orientation 

For  the  face  orientation,  the  boundaries  normal  to 
the  principal  stress  direction  crack  relatively  easily. 

The  critical  events  are  cracking  of  the  oblique  boundaries, 
at  ± 30°  to  the  maximum,  tensile  direction.  As  shown  in 
Fig.  G.4,  these  oblique  cracks  may  form  after  either  or 
both  of  their  neighbors.  First  consider  their  formation 
by  propagation.  If  the  interaction  of  the  two  normal 
cracks  is  neglected,  the  grain  boundary  intensity  is  set 
purely  by  does  not  depend  on  whether  one  or  both 

neighbors  is  cracked.  The  probability  of  an  oblique  crack 


Fig.  G.4  Forming  oblique  cracks  with  applied 
stress  normal  to  a face. 


Fig.  G.5  Approximating  the  stress  intensity 

with  applied  stress  normal  to  a face. 
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is  then  the  sura  of  the  probabilities  that  either  neighbor 
or  both  are  cracked,  times  the  probability  of  cracking  of 
the  oblique  crack  given  the  cracked  neighbors:  • 

Denoting  the  probability  of  a crack  nucleating  at  a 
neighbor  by  p^  , the  probability  that  an  oblique 
boundary  is  cracked  is 


p = (2Pn  (l-Pn)  + P„^)pe/n  • (G.18) 

The  probability  of  a crack  nucleating  on  a normal  boundary 
is  found  from  Eq.  G.l: 


The  probability  of  propagation  of  the  oblique  segment  from 
a normal  boundary  is  now  found,  noting  there  are  four 
possible  crack  sites,  two  at  each  tip,  from  an  equation 
for  the  grain  boundary  intensity  in  terras  of  the  applied 
stress  intensity  factor  and  the  orientation,  and  the  dis- 
tribution function  for  the  grain  boundary  intensity; 


kg  = k^  (1  + cos  60°) (cos  30°)/2  = 0.649  , (G.20) 


9/n 


= 1 - exp 


-4 


^0.649 

V ^k90  ^O'^l  “ ^k9L  ^L'/^l 


ml 


(G.21) 
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If  the  oblique  crack  forms  by  re-nucleation , the 
stress  at  its  center  will  be  influenced  by  the  two  adjacent 
cracks.  Taking  the  stress  field  around  a crack  of  finite 
length  from  Green  and  Zerna  (1963)  we  find  for  one  neighbor- 
ing normal  crack, 

= 0.301  S.  + 0.75  S , (G.22) 

0 0 1 2 

whereas  for  two  cracks,  with  only  one  possible  cracking 
site,  allowing  30%  interaction  effects, 

Oqq  = 0.370  + 0.75  S2  . (G.23) 
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uous  units,  counting  the  projected  distance  between  oblique 
cracks  as  one  unit,  is  found  as  before  from  Eq.  G.12. 

The  strength  of  a crack  of  length  is  found 

by  first  finding  the  stress  intensity  due  to  a chain  of 
cracks  with  projected  length  3c^  per  oblique  crack,  as 
Shown  in  Fig.  g.5.  Then 

^Ic  “ . (G.27) 

Next,  the  value  of  for  propagation  is  approximated 

as  shown  in  Fig.  G.s- 

k,  = k^^  (^(1  + cos  60»)(cos  30»)/2)  = 0.649  k^^  . (o.28) 

For  a median  value  of  this  normal  grain  boundary  intensity 
with  four  possible  crack  sites, 


0.5  = 1 - 


exp 


•4 1 ^i/"c'^v^>  - 

'k0L  ^L'/^l 


m 


^k0O 


(G.29) 


If  the  next  segment  cracks  by  nucleation,  use 
Eq.  G.l  for  the  distribution  of  strengths  and  find  the 
stress  in  the  next  grain  from  the  stress  intensity  factor 

^e0  “ ^ic  + cos  60“)  (cos  15“)/2  j = 

« 0.649  . 


(G.30) 


The  median  strength  against  nucleation  is  found  from 


0.5  = 1 


f"0.649  S^JW^'4)  - 

m 

-4  j 

^ =0  - =L  J 

(G.31) 


Since  again  kg  and  '^00  sre  assumed  to  be  determinate 
functions  of  the  same  grain  boundary  strength,  the  choice 
between  nucleation  and  propagation,  Eqs.  G.29  and  G.31,  is 
governed  by  the  parameters  , f^^g^  , , Sg  , and 

32/5^  . The  maximum  stress  is  now  the  value  of  (for 

the  given  ^ found  by  eliminating  p and  from 

Eqs.  G . 12  and  G.18,  G . 19 , G.21,  and  G.,29  for  propagation, 
or  g.19,  g.24,  g.25,  g.26  and  &«31  for  nucleation. 


g.5  formation  of  the  first  crack 


The  first  crack  is  especially  susceptible  to  the 
toe  of  the  distribution  function,  and  to  variables  in  the 
forming  processes  which  give  rise  to  it.  Nonetheless,  it 
is  of  interest  to  get  some  idea  of  the  increase  in  stress 
between  the  first  nucleation  of  cracks  and  the  final 
instability,  since  this  indicates  the  warning  signals, 
such  as  acoustic  emission,  w'hich  might  be  detected  prior 
to  fracture. 

The  probability  of  first  fracture  is  best  found 
from  the  product  of  the  probabilities  that  no  fractures 
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are  found  on  any  of  the  families  of  grain  boundaries 

having  a specific  orientation  relative  to  the  direction 
of  maximum  principal  stress. 


Vertex  Orientation 

In  a total  of  n grain  boundaries,  2N/3  have 
an  orientation  of  * 60"  to  the  direction  of  maximum 
principal  stress.  The  normal  stress  on  them  is  found 
from  Mohr's  circle: 


*^00  ^S^/4  + S./4 


(G.32) 


There  are  also  N/3  segments  normal  to  the 
stress  component,  for  v;hich 


second  principal 


^^00  “ ^2  • 


(G,33) 


The  probability  of  first  fracture  at 
stress  and  secondary  principal 

the  complement  of  the  probability  of 
of  these  families  of  planes; 


maximum  principal 
stress  S2  = C is 
no  fracture  on  either 


* « 1 - exp  [-  I [2(S^(3  t 0/4)"'  t (C  . 


(c>0)  (G.34) 


Face  Orientation 


For  the  face  orientation, 
subject  to  a normal  stress  of 


there  are  N/3  segments 
and  2N/3  segments 


jJ 


13C 


subject  to  a normal  stress  of 


^ee  “ ^1^^  30/4  . (G.35) 

Introducing  tl.ese  stress  components  into  the  extreme  value 
distribution  we  get  the  probability  of  first  fracture: 


f \ m' 

r 

b 3 ' 2 ( 

^S^{1  + 30 /4j^ 

)J  , (0-1/3)  . (G.36) 

With  either  orientation,  the  stress  level  S^.  for 
a given  probability  of  fracture  and  stress  ratio 
C = found  by  solving  Eq.  G.34  or  G.36. 

G.6  DISCUSSION  OF  RESULTS 

To  illustrate  the  above  equations,  they  were  eval- 
uated for  the  following  values  of  the  parameters : 

^k0O  ^ ^k0L  " 
m = 3,  10 
N = 10  - 100,000 
$ = 0.50. 

Figures  G.6  and  G.7  show  the  size  effect  for  two  different 
variabilities  (m  — 3,  10)  and  for  uniaxial  applied  stress. 
Qualitatively,  the  curves  are  similar  to  those  ojstained 

by  McClintock  and  Zaverl  (197  7)  for  biaxial  stress, 
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Fig.  G.6  Statistical  strength  as  a function  of  part  size 
for  a strength  distribution  with  m = 3,  S /Sf.=0, 

^keO  “ ^keL  “ 0*50.  Samples  of  5,  giving 
93.75%  confidence  limits  for  the  median. 
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Fig.  g.7  Statistical  strength  as  a function  of  part  size 
for  a strength  distribution  with  m = 10, 

Sr/Sn  = 0,  f,„,  = f, 


but  there  is  a smaller  stress  range  between  crack  initiation 
and  instability.  Again  an  extreme-value  distribution,  valid 
for  first  fracture,  would  bo  a poor  descriptor  of  maximum 
strength. 

The  fracture  loci  are  shown  in  Fig.  g.8  for  an 
exponent  m = 3 aud  a part  size  of  N = 100  grain  boundary 
segments.  Note  that  there  is  an  increase  of  the  order  of 
50%  in  the  stress  for  first  cracking  on  going  from  biaxial 
to  uniaxial  applied  stress.  Crack  closure  begins  a 
transverse  compressive  stress  of  about  half  the  tensile  stress. 

The  numerical  results  are  in  reasonably  good  agree- 
ment with  the  analytical  model  although  there  is  a slight 
discrepancy  due  to  neglecting  any  interaction  effects  in 
calculating  the  density  p of  unit  cracks  in  the  analytical 
formulations.  Larger  numbers  of  boundaries  would  have  to 
be  studied  to  show  whether  the  localization  of  cracking 
around  a crack  tip  as  it  formed  would  tend  to  promote  crack- 
ing along  the  line  of  the  crack,  or  inhibit  it  by  branch 
cracking  that  reduces  the  local  stress  concentration. 

More  accurate  modeling  should  be  done'  when  specific 
data  on  flaw  distributions  are  available  for  cases  of  practical 
interest.  Modeling  of  compression  cracks  would  require  a 
computer  program  that  would  treat  crack  closure  and  friction. 


Fig.G.8  Principal  stress  fracture  loci  for  first 

cracking  and  instability  for  an  extreme-value 

distribution  with  m = 3,  S_/S„  = 0,  f,  « = 

L'  0 ' keO 

fkeL  “ 0*50,  part  size  N >=  100. 
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H.  THE  EFFECTS  OF  SURFACE  INITIATED  FRACTURE  ON  BRITTLE 
CRACK  STATISTICS 

H.l  INTRODUCTION 

The  problem  of  c^^lcalating  the  progress  of  micro- 
cracking in  brittle  materials  towards  the  configuration 
giving  maximum  load, by  repeated  stress  analysis  of  the 
developing  micro-crack  distribution,  has  been  treated 
numerically  by  McClintock  and  Zaverl  (1975).  They  con- 
sidered grain  boundary  cracking  in  a two-dimensional  array 
of  regular  hexagon  grains  embedded  within  an  infinite  body. 
However,  brittle  materials  often  fracture  by  the  propaga- 
tion of  surface  cracks,  whether  they  are  initiated  at  sub- 
critical  stress  levels  or  are  pre-existing.  With  this  in 
mind  the  numerical  model  is  extended,  using  near 
surface  dislocation  dipoles,  to  study  the  effects  on 
brittle  crack  statistics  of  surface  or  near  surface 
initiated  fracture. 

H.2  A NUMERICAL  MODEL  FOR  SURFACE  INITIATED  FRACTURE 

The  above  numerical  model  of  the  mechanics  of  brittle 
crack  initiation  within  an  infinite  medium  is  a form  of 
boundary  integral  method.  Crack  segments  are  approximately 
modeled  by  one  dislocation  pair,  by  inserting  a relative 


displacement  dart  t^cross  the  m ' segment,  see{Fig. 

H.l.)  The  number  of  degrees  of  freedom  is  then  just  the 
current  nuxnbor  u:  c:  a<:K  segments.  The  dart  is  defined  as 
the  Burger's  vect'r  at  the  second  end  of  the  segment.  The 
resulting  tract i,;'  on  the  ii  segment  due  to  the 

til 

m cracked  so  - fmen ; i given  by  the  influence  coefficient 
T^^  , which  is  anaiagous  with  a Green's  function  (see  for 

example,  Hildebrcinu,  ,1952).  Accounting  for  the  applied 

»u  . , . u 

stress  at  infinity  t gives  the  traction  t as 


t^  - 


„;jm 


m ® 
D + t 


(H.l) 


All  segments  are  assigned  traction  f'"ee,  so 
the  darts  must  be  chosen  to  cancel  the  effects  of  the 
applied  stress.  Having  solved  for  all  darts,  the  stress 
can  be  found  ai.,'.  t.  .re  in  the  body  by  superimposing  the 
stress  fields  of  ail  cracked  segments.  The  stress  is 
calculated  at  each  gram  boundary  and  the  stress  inten- 
sities are  found  at  each  crack  tip.  Grain  boundary 
strengths  and  toaghnesscs  are  assumed  to  follow  an  extreme 
value  distribution.  Fracture  is  determined  by  the  weakest 
segment  relative  to  the  current  stress  on  it,  or  if  a 
crack  tip  is  more  critical,  by  the  lowest  segment  tough- 
ness relative  to  the  current  local  intensity  of  normal 
stress  across  a grain  boundary.  After  the  most  critical 
segment  is  found,  it  is  assumed  cracked.  Then  a new 
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I 

I 


jf'ig • H.l  The  m cracked  segment,  modeled  as  a 
dislocation  dipole,  and  the  resulting 
traction  on  the  segment. 
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elastic  solution  is  calculated,  a new  most  critical  segment 
is  found,  and  so  on. 

A model  of  surface  brittle  crack  initiation  requires 
the  normal  and  shear  stresses  at  the  free  surface  to  be 
zero.  This  condition  is  satisfied  by  using  the  stress 
field  cf  a disl ocaticn  near  a free  surface  as  given  by 
Head  (1953) . Tne  equations  giving  the  components  of 
stress  at  a point  (x,y)  for  the  more  general  case  of  a 
dislocation  at  d^^,y^)  near  an  x-axis  free  surface  are 
given  in  Appendix  „ These  equations,  with  f,  g and  h 
defined  in  Appendix  i.i  for  components  of  the  Burger's 
vector  b^  paralaol  to  and  b^,  normal  to  the  free  surface, 
take  the  abbreviated  form 


x: 


f b + f b 

XX  y y 


“9  b + q b 
yy  X X ^y  y 

^ ~ ^ ^ t h b 

X X y y 


(H.2) 


The  mean  and  deviatoric  stress  components  at  z 
due  to  a dislocation  at  z^  are,  in  complex 
form,  given  by 


)b  +(f  +g  )b 

a = YX  = X ^x  X w V 

m 2 2 

2 "^^*^xy“  2 ^ ^ (h  b +h  b ) 

z X X y y 


(H.  3) 
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The  deviatoric  part  of  the  stress  referred  to  the  normal 
and  shear  component  along  a segment  Az^  at  0 to  the  x- 
axis, using  Mohr  circle  for  the  complex  plane,  is  given  by: 


^d-  = ^d 


e-2i(e-V2)  ^ ^ 


{ Az^ 


(H.4) 


Giving  the  components  of  stress  in  n,s  coordinates  as 


a 

a 

a 


nn 

ss 

ns 


°m  " ) 


Im(a^, ) 


(H.5) 


The  traction  t^  on  the  segment  due  to  the 

t h 

m cracked  segment  in  n,s  components,  and  accounting 
for  the  Burger's  vector  at  each  end  of  the  dart,  is  given, 
after  some  manipulation  of  Eqs.  H. 3 to  H,5,  by 


+ It 


m 

' m=l 


+ a 


m 


+ CT 


(- 


(H.6) 


where  the  influence  coefficients  are: 
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n=l 
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“ cID^L 


<V’x'  - i'V^y'  - (<V9x> 


• ^<V^y')  sf  - “‘Vii'y! 


Az 
y'  Az 


ym 


Y"  1 ‘>!="  r „ , / 

- ^'fy%)  - ( 

n=  1 


(H  .7) 


(f  ~g  ) 

X 


+ i (f  -g  )1 

y y 7 


^ - i2(h  +ih  ) ^ 
Az  X y'  Az 


The  Burger's  vectors  are  related  to  the  darts  through  the 
parameter  db^/dD^^  which  takes  on  a value  of  -1  for  n=l 
and  +1  for  n=2. 

The  influence  coefficients  of  Eq.  H.7  are  now  used 
in  the  Brittle  Statistics  Program  BRTST3  (McClintock  and 
Zaverl,  1975)  to  find  the  Burger's  vectors  for  a given 
crack  distribution  and  then  Eqs.  h.3  are  used  to  find 
the  stress  for  next  cracking. 

To  minimize  the  errors  in  modeling  a surface  crack 
by  a relative  displacement  dart  its  ends  are  recessed  in 
from  the  ends  of  the  segment  by  a fraction  Az^Az"’ 
of  the  segment  length.  The  accuracy  of  the  model  compared 
with  theoretical  and  the  optimum  value  of  Az^/Az^  to 
choose  are  considered  in  detail  in  Appendix  H.2. 


For  a 


value  of  0.10  an  error  of  loss  than  10%  was  achieved 
in  the  mid-segment  stresses  on  neighboring  segments  and 
those  far  from  the  or  nek;  tip.  The  computed  intensity 
would  be  low  by  a factor  of  20%  . However,  this  can  be 
compensated  for  by  adjustment  of  the  computational  para- 
meter that  relates  intensities  to  the  grain  boundary 
strengths.  Furthermore,  the  requirements  of  a stress- 
free  surface  are  satisfied.  Also,  if  the  crack  is  away 
from  the  surface  by  a distance  very  much  greater  than  its 
length,  the  computed  stress  field  becomes  that  of  an 
embedded  crack  I'/ithin  an  infinite  medium. 

The  computed  effects  on  the  statistics  of  brittle 
materials  for  surface  or  near  surface  fracture  are  discussed 
after  the  analytical  theory  is  extended,  for  comparison, 
in  the  next  section, 

H,.  3 ANALYTICAL  MODEL 

The  extension  to  study  surface  effects  of  the 
analytical  formulation  of  Section  G.  4 is  straight-forward. 

As  before  the  models  estimate  the  maximum  strength  by 
finding  the  density  of  unit  cracks  as  a function  of  the 
applied  stress.  The  probability  of  aggregation  of  these 
cracks  is  estimated.  The  maximum  strength  being  achieved 
when  crack  length  and  stresses  are  sufficiently  high  to 
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Cause  crack  propagation  even  with  median  strengths  at  the 
crack  tip. 

The  density  p of  unit  cracks  under  a given 
applied  stress  and  corresponding  longest  contiguous  chain 
of  , having  probability  of  forming,  are  unaffected 

by  the  surface.  However,  the  intensities  due  to  surface 
cracks  and  the  number  of  possible  cracking  sites  at  each 
tip  are  different  from  cracks  embedded  within  the  material 
and  with  this  in  mind,  the  models  are  extended. 

Vertex  Orientation 

A surface  crack,  for  this  orientation,  is  shown  in 
Fig.  H.2,  comprising  of  cracked  grains.  The  average 

direction  of  the  crack  is  assumed  normal  to  the  applied 
stress  and  nas  a projected  length  of  N 2c, cos  30° 

The  intensity  of  the  crack  shown,  estimated  from  fracture 
mechanics  (Wigglesworth, 1957) , is 

kic  = 1.12  y2c^N^cos30°  . (H. 8) 

The  normal  grain  boundary  intensity  k^  , for  propagation 
of  the  crack,  is  reached  when 

" ^Ic  ((1+COS30')  (cosl^)/2j  = 0.901  k^^  . (H.9) 


For  a median  value  of  intensity,  noting  that  th-re  is  only 


Fig.  H.2  Approximating  the  stress  intensity  with  applied 
stress  for  a surface  crack  in  the  vertex  orient- 
ation. 


Fig.H.3  Approximating  the  stress  intensity  with  applied 

stress  for  a surface  crack  in  the  face  orientation. 
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one  possible  cracking  site  at  the  tip  of  a surface  crack, 
propagation  occurs  when 

.901x1.12  S-  y2c-N  cos30"-f,  ^ s, 

^ X C IC  D Ij  Ij 

*"k90^0  "^'^l  " ^kOL^L  '^l 

For  the  next  segment  to  crack  by  nucleation  the  midsegment 
stress  on  the  next  grain  from  the  normal  grain  boundary 
intensity  is 


0.5  = 1 - exp 


r 


So  for  a median  strength  of  the  next  grain,  crack  growth 
by  renucleation  occurs  when 


0.5  = 1 - exp 


/'0. 901x1. 12  S,  y U cos30°-S, 

I 1 c L 

I =o-"l 


m 


(H.12) 


The  maximum  stress  is  obtained  as  in  Section  G.4 
but  using  Eqs . H.IO  instead  of  G.15  and  H.12  instead  of 
G.  17  for  crack  growth  by  propagation  and  renucleation 
respectively. 


Face  Orientation 

For  this  orientation  the  projected  crack  length 
is  3c^N^  as  shown  in  Fig.  H.3.  The  corresponding  value 
of  normal  grain  boundary  intensity  is  then 


.10) 
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(H.13) 


kg  = ^(l+cos60)  (cos3(f)/2^  ~ 0.649x1.12  S 3 c ^ • 


In  this  case  there  are  two  possible  cracking  sites  at  the 
tip  of  a surface  crack.  For  a median  intensity  kg  , 
propagation  occurs  when 


0.5 


1 - exp 


-2 


'0.649x1.12 


^keo^o 


Sj  '/^‘^l^c  ~ ^kOL^L^l'^ 
•^1  “ ^kOL^L  '^1  ■'  _ 


(H.14) 


For  a median  strength  of  the  next  grain  crack  growth  by 
renucleation  occurs  when 


0.5  = 1 - exp 


-2 


'0.649x1.12  y3/2 


. (H.15) 


The  maximum  stress  is  obtained  as  in  Section  G.4 
but  using  Eqs.  h , 14  instead  of  G.29  and  H.15  instead  of 
G,31  for  crack  growth  by  propagation  and  renucleation 
respectively. 


H.4  DISCUSSION  OF  RESULTS 

Part  sizes  up  to  1000  grain  segments  were  studied 
using  the  numerical  model  for  two  different  variabilities 
(m=3,10).  The  results  are  shown  for  the  two  grain 
orientations  in  Figs.  H.4  and  H.5  and  are  in  fairly 
good  agreement  with  the  values  predicted  by  the  analytical 
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1.0 


Fig.H.4  Statistical  strength  for  the  face  orienta- 
tion as  a function  of  part  size  for  surface 
initiated  fracture.  Strength  distributions 
of  m = 3 and  10  with  S^^/Sq  = 0.  Samples  of 
5,  giving  93.75%  confidence  limits. for  the 
median. 
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Fig.  H.5  Statistical  strength  for  the  vertex  orienta- 
tion as  a function  of  part  size  for  surface 
initiated  fracture.  Strength  distributions 
of  m = 3 and  10  with  S^^/Sq  = 0. 
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models.  Maximum  strengths  for  surface  initiated  fracture 
are  10-20%,  depending  on  orientation  and  variability, 
below  those  for  embedded  fracture.  Consequently,  first 
cracking  governs  fracture  for  part  sizes  larger  than  embed- 
ded fracture,  but  such  part  sizes  are  still  small  compared 
with  even  laboratory  size  test  specimens.  The  relatively 
small  drop  in  maximum  strengths  suggests  that  the 
statistical  treatment  of  brittle  materials  is  somewhat 
independent  of  whether  fracture  is  initiated  in  the  bulk 
or  at  the  surface  of  the  material  as  long  as  the  relevant 
distribution  of  grain  boundary  strengths  is  used. 
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appendix  h.i. 


STRESS  FIELD  OF  A DISLOCATION  NEAR  A 
FREE  SURFACE 


The  numerical  surface  model,  developed  in  Section 
H,  requires  the  normal  and  shear  stresses  at  the  free 
surface  to  be  zero.  This  requirement  is  satisfied  by 
using  the  stress  field,  as  derived  by  Head  (1953),  of  a 
dislocation  near  a free  surface.  These  equations  des- 
cribing the  stresses  at  a point  (x,y)  for  the  more  general 
case  of  a dislocation  at  (x^,y^)  near  an  x-axis  free 
surface  are  given  below. 


For  a Burger's  vector  parallel  to  the  free 

surface : 


o = b 

XX  X 


E ' 

f 

— Y 

1 

+ W I ■ 

+ 2y 


(w^+u^) 


vw^-6vwu^+u^ 


n 


(w^+u^)  ^ 


II 


a = b 

yy  X 


E' 


- V 


[ v^-u^  \ ^ ( w^-u“  \ 

l(vW)2j  ■ 

f (3y+y^)w^-6ywu^-u‘^ 


(H.I.l) 


- 2y 


(H.I. 2) 
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f 2_  2 

“xy  = ^ - ,"2^“2,2 


2 2 
w -u 

, 2^  2>2 
(w  +u  ) 


^ 2 2 V 

I 3w  -u 


(H.I.3) 


where  u = ^ “ ^"^n'  ^ plane  strain 

with  Young's  modulus  E and  Poisson's  ratio  v 


4tt{1-v  ) 


(H.I.4) 


For  a xlurger's  vector  perpendicular  to  the  free 

surface, 


°xx  ==  - 


2 2 
V -u 

,2  2 '> 

(v  +u  ) 


2 2 
w -u 

/ 2 , 2,2 

(w  +u  ) 


~ 4y^ 


{2Vj^-y)w^+  (3y+2y^)  u^ 

, 2^  2,3 
(w  +u  ) 


(H.I. 5) 
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(H.I, 6) 
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(H.I.7) 


xy 


= b. 


f 

E' 

V 
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2 2 
V -u 
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/ 2 2 
w -u 

,2  2,2 

^(w  +u  ) 


- 2y  f-v«^+6ywu^-uM 


The  above  equations  for  cr  , a and  a have 

XX  yy  xy 

been  written  so  that  the  first  term  corresponds  to  the 
stress  field  of  a dislocation  with  Burger's  vector  b , 
and  the  second  term  to  the  field  of  a dislocation  with 
Burger's  vector  . -b  at  the  image  point  (x^,-y^).  It  can 
be  seen  that  the  stress  distribution  is  not  just  the 
superposition  of  the  stress  fields  of  the  dislocation 
and  an  image  but  the  third  term  in  each  of  the  above 
equations  is  necessary  to  give  the  required  zero  normal 

stress  on  the  boundary. 


The  above  Eqs.  H.I.l  to  H.1.3  and  Eqs.  H.1.4  to 

H.I.7  riay  be  superimposed  for  a dislocation  having 

components  of  Burger's  vector  b and  b and  for 

X y 

convenience  are  written  in  abbreviated  form  as 


a 

a 

a 


- b f + b f 


XX  XX 

y y 

= b g + 

b jg„ 

yy  x^x 

y y 

= b h + 

b h 

xy  XX 

y y 

(H.I.8) 
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APPENDIX  H.II. 


ACCURACY  OF  THE  SURFACE  MODEL 


In  this  Section  the  accuracy  achieved  in  modelling 
Mirface  cracks  by  dislocation  dipoles  is  considered  by 
comparing  computed  values  of  stresses  and  intensities  with 
th-,  theoretical. 


A crack  of  length  Az^  is  modelled  as  a near 
surface  dislocation  dipole,  by  inserting  a relative  dis- 
placement dart  D , with  ends  recessed  an  amount  Az'^  , 
as  shown  in  Fig.H.Il.l  Then  the  opening  mode  stress 
c^dip  V ' found  from  Eq.H.I.l  , is 


c 


2 (Az^-Az^‘)  _ 

y*^-  (Az"^-Az")^ 

- 2Az^ 

2 , A n. 

y -(Az  ) 


2 (Az”^-Az^)  (y-  (Az^-Az^) 
(y+  (Az^-Az^)  ) ^ 

, 2Az’^(y-Az’^) 

2 "TTprrj— 

(y+Az  ) 


) 


(H.II  .1) 


For  convenience  Eq.H.ii.l  is  writt»^n  as 


dip' 


=E'D 


F (y  ,Az^,Az’^) 


(H.II .2) 


The  midsegment  stress  due  to  the  dipole  must  cancel  the 

applied  stress  a so  that 

00 


a 


00 


y=0 . 5Az 


m 


= E'D  F(y=0.5Az^)  . 


(H.II. 3) 


Figure  H 




SURFACE  CRACK 


DIPOLE  MODEL 


II.  1 Surface  crack  and  the  corresponding 
dislocation  dipole  model. 
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Solving  for  the  dart  D 


in  Eg.  H.II.3  gives 


D 


o 

CO 


E'  F(y=0.5Az"') 


(H.II.4) 


Adding  in  the  applied  stress 
computed  opening  mode  stress 
from  Eqs.H.ii.2  and  H.II.4 


at  infinity,  the  total 
'^comp  y is  found 


a 

comp 


F(y,A2^,Az^~^)  ^ \ 

FCy^O.SAz"')  ; 


(H. ri. 5) 


S_tress  at  a Crack  Kalf-length  ahead  of  the  Crack  Tip 

Evaluation  of  the  series  solution  derived  by 
Wigglesworth  (1957),  gives  the  opening  mode  stress  a 

theo 

a distance  0.5Az"^  (y^l.SAz'^)  ahead  of  a surface  crack  as 

^theo  ^ 1-378  . (H.II.6) 

•The  ratio  of  the  computed  stress  given  by  the 
Qipole  model  to  the  theoretical  at  a crack  half-length 
ahead  of  the  crack  tip,  from  Egs.H,Ii.5  and  H.ii.6  is  then 


a 

comp  _ 1 

^theo  " 


F(y=1.5Az^) 
F(y=0.5Az"^)  j 


(H.II.7) 


Eq.  H.II.7  is  plotted  in  Fig.  H.II.2  for  different  values 

A n m 
of  Az  /Az  . 


Stress  at  Large  Distance  from  the  Crack  Tip 

The  theoretical  opening  mode  stress  given  by 
Wigglesworth  (1957)  ahead  of  a surface  crack,  length  Az^, 
at  large  distance  from  the  crack  tip,  approaches  the  value 
given  by  Green  and  Zerna  (1963)  for  an  embedded  crack, 
length  2Az"',  and  is 


°^theo 


a 

00 




/A  HI  A 2 

(Az  ) 


1/2 


(H.II.8) 


Subtracting  out  the  applied  stress,  the  ratio  of 
computed  stress  to  theoretical,  from  Eqs.  H.II.5  and  H.II.8 
at  large  distance  from  the  crack  tip  is  then 


f ° 

comp 


[*^theo 


T-.  / HI  n . 

F(y,  z , z ) 
F(y=0.5  z^") 


(Az"')  ^) 


- 1 


(H. II .9) 


The  error  in  computed  stresses  for  y>>Az"'  is  given  by 
Eq.  H.II.9  and  is  plotted  in  Fig.  H.II.3. 


Intensities  of  the  Dipole  Model 

The  theoretical  opening  mode  intensity  of 

a surface  crack  is  gi'^en  by  Wigglesworth  (1957)  as 
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Fig.  HII.2Ratio  of  the  computed  stress  to  theoretical  at  a 
half  segment  length  ahead  of  the  crack  tip. 


Fig.  HII . 3 Ratio  of  the  computed  stress  to  theoretical  far 
from  the  crack  tip. 


Fig. HI I. 4 Ratio  of  the  computed  intensity  to  theoretical. 


(H.II.IO) 


k.  . = l.±2 

cheo 


a /Az 
00 


m 


The  intensity  k^  found  from  the  crack  opening 
displacement  2u  at  a distance  6r  from  the  tip,  given  by 
Rice  (1972) , for  plane  strain  is 


k,  = (2u)  . (H.II.ll) 

/T^ 

The  computed  intensity  k is  found  from  Eg.  H.II.ll  by 

approximating  the  opening  displacement  by  the  dart  at  a 
crack  haif-length  back  from  the  tip.  Dividirg  Eq.  H.II.ll 
by  Eq.  H.II.IO,  and  substituting  Eq.  H.II.4  for  the  dart, 
gives  the  ratio  of  computed  intensity  to  theoretical: 


comp  _ 

^theo  1.12Az^‘  F (y=0 . 5Az'^') 


^ fi  ■ j. . I . 


The  error  in  computed  intensities  is  given  by  Eq.  H.II.12 
and  is  plotted  in  Fig.  H.II.4. 


Discussion  of  the  Numerical  Surface  Model 

Figs.  H.II.2,  H.II.3,  and  H.II.4  compare  the  computed 
values  of  stresses  and  intensities  with  theoretical  for 
different  values  of  Az”/Az"'.  For  a value  of  0.10,  an  error 
of  less  than  10%  was  found  in  the  midsegment  stresses  on 
neighboring  segments  and  those  at  'infinity'.  The  computed 
intensity  v'ould  be  low  by  a factor  of  20%.  However  this  can 
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be  compensated  for  by  adjustment  of  the  computational  para- 
meter that  relates  intensities  with  the  grain  boundary 
strengths.  Note  that  for  an  embedded  crack  the  optimum  value 
of  Az^/Az”^  to  choose  is  0.14. 
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I.  STRENGTH  DISTRIBUTION  FROM  EXPERIMENTAL  DATA 


I.l  INTRODUCTION 

Accurate  strength  distribution  curves  assist  the  designer 
in  predicting  the  behavior  of  brittle  components.  They  can  also 
give  insight  into  ways  in  which  manufacturing  techniques  can  be 

altered  to  improve  the  ultimate  performance  of  the  finished 
product. 

The  flaw  density  distributions  of  brittle  materials  can  be 
obtained  from  many  different  kinds  of  tests,  such  as  the  hardness 
indentation  test  and  the  three  or  four  point  bend  test.  Recently 
growing  interest  in  acoustic  emission  has  led  to  the  suggestion 

that  strength  distribution  can  also  be  obtained  by  acoustic 
emission . 

1.2  FLAW  DISTRIBUTION  FROM  HARDNESS  INDENTATION 
State  of  the  Art 

Hardness  indentation,  especially  spherical  indentation  (also 
known  as  Hertzian  indentation),  has  been  favored  by  many  resear- 
chers in  their  investigations  of  brittle  materials.  Because  of 
Its  ease  of  application,  and  because  only  a small  area  of  the 

specimen  is  affected,  numerous  data  can  be  obtained  within  a 
short  period  of  time. 

Many  workers  (Argon,  1959,  Suoov,  1962,  Oh  and  Finnic,  1967, 
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Tsai  and  Kolsky,  1967,  Hamilton  and  Rawson,  1970,  and  Poloniecki 
and  Wilshaw,  1971)  have  attempted  to  calculate  the  flaw  distri- 
bution curves  from  the  experimental  data,  but  their  methods  were 
either  too  complicated  or  of  an  empirical  nature. 

Earlier  in  this  research,  Matthews,  McClintock,  and  Shack 
(1976)  developed  a simple  and  straight-forward  method  which  made 
use  of  only  the  fracture  load  from  either  the  three  point  bend 
test  or  the  spherical  indentation  test.  In  the  spherical  inden- 
tation test,  if  the  specimen  with  elastic  properties  and  , 

is  indented  with  a sphere  of  radius  R and  with  elastic  proper- 
ties E2  and  V2  , the  contact  radius  a which  is  a function 
of  the  applied  load  P , is 


3PR 

4 

El  E2  JJ 

Outside  the  contact  circle  r = a the  radial  stress  at  the 

surface  is  tensile: 


0 


rr 


(r  > a)  . 


(1.2) 


The  maximum  tensile  stress  is  at  the  contact  circle: 


0 


m. 


(1.3) 


If  the  flaw  density  g(a)  is  defined  such  that  g(o)  do 
:.s  the  number  of  flav/s  per  unit  area  of  strength  between  o and 
o + do  , and  if  the  probability  of  failure  for  a stress  less  than 
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'^m  ^ the  spherical  indentation  test, 


g (a  ) = 

^ m , 

TTa^ 


$"(0^^)  [l-$  (a^^)  ] + [$'  (a^)]"' 


[1  - <!'  (Oj^)  ] ^ 


(1.4) 


4>(aj^)  can  be  ordered  from  the  ::xperimental  data,  and  by  dif- 
ferentiating it  twice,  g(%)  is  obtained. 

In  order  to  record  the  load  at  which  ring  cracking  starts, 
the  indented  area  must  be  observed  continuousl;y  during  the  test. 
If  the  specimen  is  transparent  (e.g.  glass) , observation  can  be 
made  with  an  optical  microscope  from  the  side  or  the  bottom  of 
the  specimen.  If  the  specimen  is  opaque  (e.g.  silicon  nitride) , 
Evans  and  Wilshaw  (1976)  suggested  that  the  initiation  of  cracks 
can  be  detected  by  acoustic  emission.  The  following  alternate 
method  was  proposed  by  Lawn  (1968).  With  the  maximum  load  as 
variable,  a systematic  series  of  indentations  was  performed. 

The  specimen  was  then  viewed  to  find  out  which  indentations  had 
produced  ring  cracks.  The  critical  fracture  load  could  be  de- 
termined to  within  about  5%. 

Mathews  et  al  applied  their  analysis  to  the  test  results  of 
Argon,  Hori,  and  Orowan  (1960)  to  obtain  the  flaw  density 
curves.  Their  method  has  also  been  applied  successfully  by 
Evans  and  Wilshaw  (1977) . Evans  and  Jones  (1977)  have  even 
extended  their  analysis  to  the  tensile  test,  the  four  point 
bend  test,  the  expanded  ring  test,  and  the  concentric  ring 
test . 
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The  Triction  Effect  wit-h  Spherical  Indente r 


It  was  pointed  out  by  Rawson  (1977)  that  the  analysis  of 
Matthews  et  al  overlooked  the  frictional  forces  at  the  interface 
between  the  indenter  and  the  substrate/  when  these  components 
have  different  elastic  properties.  These  forces  will  change  the 
stress  distribution  on  the  surface  of  the  specimen,  and  the  ap- 
parent flaw  density  curve  will  be  affected.  Johnson,  O'Connor, 
and  Woodward  (1973)  have  calculated  the  stress  distribution  for 
two  extreme  cases  of  com.plete  slip  and  no  slip.  In  reality,  the 
indenter  will  slip  in  an  annulus,  the  exact  size  and  location  of 
which  depends  on  the  elastic  properties  and  the  friction  coef- 
ficient. The  maximum  stress  when  a rigid  ball  is  indented  on  a 
less  rigid  surface  (e.g.  steel  indenter  on  glass  surface)  occurs 
not  at  the  contact  circle  as  predicted  by  the  Hertzian  theory, 
but  at  a short  diotar.;,e  outside  it.  Thus,  the  ring  crack  should 
start  outside  the  contact  circle.  This  finding  is  in  qualita- 
tive agreement  with  the  experimental  results  of  Argon  et  al  (1960) 
as  well  as  that  of  other  workers,  although  some  of  the  observed 
effect  is  associated  w^th  the  statistical  variability  of  frac- 
ture (Oh  and  Finnie,  1970,  Hamilton  and  Rawson,  1970). 

Johnson  et  al  have  prepared  two  curves  which  are  reproduced 
here  in  Figs.  I.l  and  1.2.  From  these  two  curves,  once  the  con- 
stant < in  terms  of  the  elastic  properties: 


(1+v^)  (l-2v^)  ^^"^^2^ 
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Fig.  I-l-  Maximum  values  of  the  radial  tensile  stress  (a)  as  a 
function  of  k (no-slip  theory)  and  (b)  as  a function 
of  p (complete  slip  theory).  An  estimated  curve 
(chain  line)  is  interploated  between  the  complete  slip 
theory  and  no-slip  theory  for  k = 0.4. 


Fig*  1.2  Position  of  the  maximum  radial  tensile  stress  (a)  as 

a function  of  k (no-slip  theory)  and  (b)  as  a function 
of  p (complete  slip  theory).  An  estimated  curve  (chain 
line)  is  interpolated  for  K = 0.4. 


i 
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and  the  friction  coefficient  li  are  known,  a reasonable  in- 
terpolation between  the  complete  slip  theory  and  the  no  slip 
theory  is  possible.  E’or  example, with  a steel  indenter  and  a 


clean  glass  surface,  Johnson  et  al  measured  an  initial  friction 
coefficient  of  only  0.12,  ana  assuming  = 2.07  x lO^MN/m^, 


^gl ass 
0.364. 


^ 6.90  X 10'*MN/m^  , ^g^eel  ^ 0*^0'  ^glass  ^ k is 

From  Figs.  I.l  and  1.2,  it  can  be  said  that  the  com- 


plete slip  theory  would  be  a better  approximation.  The  stress 


distribution  assuming  complete  slip  is  plotted  in  Fig.  1.3. 

•3r 


radial  distance,  r/a 

Fi;;.  i.l.  Surface  radial  GtrcGG  distribution  of  Gtecl  ba.ll  iadtntin'-,  on 

Sln.GG  specimen.  (s)  Hertzian  distribution,  friction  neglecti.d. 
(b)  Friction  modified  ntrcGG  distributinn  ar.GU'in^  conplete  : liii, 
!■  = 0.16',  and  p = 0.12. 


The  maximum  tensile  stress  is  ^ocated  around  r/a  = 1.12. 

To  minimize  the  friction  effect,  the  indenter  and  the  speci- 
men should  have  similar  elastic  properties,  or  the  interface 
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should  be  lubrioated  with  high  pressure  greases  sueh  as  those 
that  contain  graphite,  niolybdenua,  disulphide,  or  fatty  acids. 

On  the  other  hand,  the  analysis  developed  by  Matthews  et  al 
would  be  .uch  n.ore  useful  if  it  could  be  modified  to  account  for 
the  friction  by  assuming  a stress  distribution  that  reached  a 
um  at  r - r^(rj^>a).  Such  work  is  being  carried  out  as 


part  of  a subsequent  project. 


Although  a spherical  indenter  is  ideal  for  studying  surface 
flaws  in  brittle  materials,  a different  shape  of  indenter  would 
be  needed  to  study  sub-surface  flaws.  Such  an  indenter  should 
give  no  stress  on  the  specimen  surface  so  that  fracture  will  not 
originate  at  the  surface.  A cylindrical  indenter  with  its  axis 
parallel  to  the  surface  (Fig.  1.4)  does  not  generate  any  stress 


y 

Fig.  1.4.  Plane  strain  cylindrical  indentation. 
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on  the  surface  outside  the  contact  area  and  produce  only  com- 
pressive stresses  in  the  interior.  A cylindrical  indenter  is 
also  useful  in  modelling  the  contact  between  a car  wheel  and  a 
rail,  or  between  a roller  bearing  and  the  inner  or  outer  ring. 

The  elastic  stress  distribution  of  piano  strain  cylindrical 
indentation  can  be  found  in  McClintock  (1977)  or  Poritsky  (1950) 
in  complex  form: 

0 t 0 r 

^ 2P  z 

2 Tra 


(1.6) 


(1.7) 


where  z = x + iy  , 

and  a - contact  half-width 

2 

(1.8) 


^ 4PR 

1 - 1 - 

TT 

When  a brittle  material  is  under  biaxial  stresses,  Griffith 
(1924)  showed  that  when  the  smaller  stress  is  in  the  range 

— *^2  — *^1  ’ f^^^^ture  takes  place  when  the  algebraically 
larger  stress  reaches  the  tensile  strength  of  the  ma- 
terial. When  O2  is  in  the  range  02<  > fracture  takes 

place  when 
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(Oi  - 02)^  + Sa^  (a^  + 02)  = 0 


(1.9) 


Solving  for  o. 


^^1  - ^2^' 
8 (oj^  + 02) 


(I. 10) 


With  ^he  principal  stresses  and  02  calculated  from  the 

stress  distribution  functions  (Eqs.  1.6  and  1.7),  the  maximiam 
value  for  0^  is  found  to  occur  along  the  axis  at  a depth  of 
y/a  = 1.25.  McClintocl<  and  Walsh  (1962)  modified  Griffith's 
criterion  by  considering  the  effect  of  craclc  closure  with  fric- 
tion. If  0^  is  the  normal  stress  to  close  a craclc,  and  y is 
the  coefficient  of  friction,  they  showed  that  the  condition  for 
fracture  after  the  cracks  have  closed  is 


y (02  + 0^  - 2c^)  + (oj^  - 0,)  ,/TT 


y 2 = 4 


(I. 11) 

By  assuming  that  the  first  crack  caused  instability,  they  were 
able  to  fit  past  experimental  data  well  with  values  of  y of 
1.0,  and  0^  of  -3a^  . Using  these  values,  the  maximum  value 
of  0^  IS  found  to  be  along  the  axis  at  y/a  = 2.05  . Therefore, 
either  the  Griffith's  criterion  or  the  modified  criterion  pre- 
dicts fracture  in  compression  to  start  along  the  axis  in  the  in- 
terior of  the  specimen. 

Plane  strain  cylindrical  tests  were  performed  by  Thomas  and 
Hoersch  (1930)  who  tested  on  steel  plates,  and  Kolsky  (1952)  who 
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tested  on  both  glass  and  "Perspex"  (polymethylmethacrylate) . 
Kolsky,  using  steel  indenters  of  radius  of  curvature  of  12  in- 
ches, found  that  the  glass  specimens  fractured  along  the  outer 
edges  of  the  area  of  contact.  A similar  problem  of  compressing 
glass  cylinders  between  two  steel  plates  was  investigated  by 
Hooper  (1971).  Hondros  (1959)  showed  that  both  principal  stres- 
ses become  comipressive  in  the  loading  zones,  with  a surface 
value  equal  to  the  applied  normal  stress.  The  maximum  tensile 
stress  occurs  at  the  center  of  the  cylinder.  Therefore,  frac- 
ture is  expected  to  initiate  from  the  center  of  the  cylinder. 
However,  Hooper  (1971)  found  from  experiments  that  fracture 
started  from  the  surface  of  the  cylinder. 

Hooper  suggested  that  the  frictional  forces  at  the  interface 
due  to  the  mismatch  of  elastic  properties  induced  tensile  stres- 
ses at  the  edges  of  the  contact  zone.  He  solved  the  surface 
stress  distribution  for  two  extreme  cases  of  complete  slip  and 
complete  adhesion,  but  his  results  are  too  lengthy  to  be  inclu- 
ded here.  Applying  his  equations  to  the  problem  of  a glass  sur- 
face indented  with  a steel  cylinder,  it  can  be  shown  that  tensile 
stresses  exist  at  the  edges  of  the  contact  zone.  He  attempted 
to  reduce  the  friction  effect  by  applying  high  pressure  lubri- 
cants. He  observed  that  the  failure  load  using  dry  platens  was 
approximately  35%  higher  than  that  obtained  using  greased  pla- 
tens, a result  which  is  in  general  agreement  with  the  theorv. 

Hooper  also  suggested  that  since  a cylinder  has  a finite 
length,  the  stress  distribution  is  actually  three-dimensional, 
and  the  plane  strain  distribution  is  not  adequate  to  describe 
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the  problem.  The  solution  for  the  elliptical  pressure  distri- 
bution in  an  elastic  half-space  due  to  a rectangular  contact 
area  was  given  by  Kunert  (1961).  His  analysis  predicted  the 
existence  of  tensile  stresses  along  the  edges  of  the  contact 
zone  if  the  length  of  the  cylinder  is  less  than  approximately 
40  times  the  contact  half-width. 


1.3  STRENGTH  DISTRIBUTIONS  FROM  ACOUSTIC  EMISSION 

Acoustic  Emission  and  Brittle  Materials 

Acoustic  emissions  are  the  stress  waves  originating 
from  strain  energy  release  as  a result  of  various  physical 
and  structural  changes  within  the  material  caused  by  the 
applied  stress. 

The  strain  energy  released  in  the  formation  or 
growth  of  micro-cracks  prior  to  maximum  load,  observed, 
for  example,  by  Evans  (1974) , is  thought  to  be  the  cause 
of  acoustic  emission  in  many  polycrystalline  ceramics  and 
possibly  other  materials  that  exhibit  barely  any  dis- 
location motion  at  room  temperature.  Also,  as  would  be 
expected,  no  acoustic  emission  is  observed  from  amorphous 
brittle  materials,  except  at  the  onset  of  catastrophic 
failure  (Noone,  1973).  Others,  such  as  Graham  and  Alers 
(1973),  have  observed  a relationship  between  detected 
acoustical  events  and  fractured  grains  and  for  Lucalox 
found  a one-to-one  correspondence.  Noone  (1973)  in  the 
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same  material  similiarly  attributed  a single  acoustical 
event  with  the  parting  of  two  grain  surfaces. 

The  above  strongly  suggests  that  for  many  brittle 
materials  any  detected  acoustic  emission  is  due  to  the 
fonr.ation  of  stable  micro-cracks  and  not  due  to  other 
mechanisms  such  as  dislocation  motion. 

Obtaining  the  Statistical  Parameters 

The  formation  of  stable  micro-cracks  at  stresses 
below  the  ultimate  fracture  load  suggests  that  the  general 
method  of  McClintock  and  Zaverl  (1976)  can  be  used  to 
obtain  the  statistical  parameters  m,  and  that 
govern  crack  initiation  in  grain  boundaries,  idealized  as 
an  extreme  value  distribution  of  the  third  asymptotic  kind. 

The  approach  to  be  taken  is  similiar  to  the  analysis  of 

Evans  etal  (1974).  However,  it  is  shown  here  that  the  number 

of  micro-cracks  is  largely  a function  of  the  applied  stress 

level  and  that  time-dependent  effects  are  usually  only  of 

minor  importance.  The  acousv.lc  emission  data  of  Evans  et  al 

(1974)  is  shown  in  Fig.  1.5  plotted  in  terms  of  'count-rate' 

against  applied  stress  for  polycrystalline  alum.ina  ('Lucalox') 

for  a test  conducted  at  constant  stress  rate  S.  From  the 

data^with  r the  gradient  of  the  straight  line,  the  count- 

rate  ^ is  given  in  terms  of  the  stress  S and  constant 
dt 

C by 
s ■' 
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Fig.  1.5  Acoustic  emission  count  rate  against  applied 

stress  for  Lucalox  by  Evans,  Russel  and  Linzer, 
(1974)  . 
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(1.12) 


r log  S - r log  C^. 

The  observations  of  Graham  and  Alers  (1973)  and 
Noone  (1973)  suggest  that  the  number  of  counts  n is 
related  to  the  number  of  micro-cracks  in  terms  of  a 
constant  D so  that 


dn  _ 
dt  “ 


(1.13) 


It  turns  out  that  the  slight  stress  dependence  of  D, 
as  proposed  by  Tetelman  (1973),  affects  the  calculated 
r umber  of  micro-cracks  by  only  10  - 20%,  and  may  be  ne- 
glected. Combininq  Eqs„1.12  and  1.13  leads  to 


d\ 

d? 


1 

C D 
s 


(1.14) 


Integrating  Eq.  1.14,  taking  the  lov/er  limit  as  zero 


(J'‘  vcs 


1 S 

t D(r+1)  ’ 

s 

The  analytical  theory  of  McClintock  and  Zaverl 

given  an  equation  of  similiar  form  for  = 0 

L 

size  N grain  boundary  segments, 


(1.15) 

(1975) 
and  part- 


(1.16) 


From  Eqs.  1.15  and  1.16  the  extreme  value  exponent  m is 
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in 


r+1 


and  the  parameter  Sq  is 

, ^l/(r+l) 

S„  = C^[N(r+l)DJ 

For  the  unpolished  specimen  m=7  . The  effect  of 
polishing  was  to  significantly  reduce  the  amount  of  acoustic 
activity  as  well  as  increasing  m to  about  25  . This 
suggests  that  many  of  the  low  strength  grains,  damaged 
during  manufacturing  and  machining^  have  been  removed  from 
the  surface  by  polishing.  Removal  of  the  low  strength 
cracks  during  polishing  would  account  for  the  low  acoustic 
activity  from  optically  finished  samples  of  zinc  selenide 
as  observed  by  Evans  (1976).  In  fact  the  amount  of  precur- 
sor acoustic  emission,  being  first  detected  at  90%  of 
the  fracture  stress,  was  of  such  small  size  that  applica- 
tion of  the  above  m-^thod  was  not  possible.  It  is 
recommended  that  numerical  modelling  of  unpolished  sur- 
faces be  carried  out  using  the  surface  model,  as  described 
in  Chapter  3,  with  prexisting  cracks. 
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J.  FLAW  DETECTION  BY  SCATTERED  LIGHT 
J.l  INTRODUCTION 

A flaw  detection  technique  is  useful  when  fracture 
occurs  through  the  propagation  of  preexisting  cracks. 
Attempts  can  then  be  made  to  detect  these  cracks  at  test 
stresses  low  enough  to  avoid  premature  component  failure. 

It  was  shown  by  McClintock  .and  Hahn  (1975)  that  an 
impractical  number  of  tests  would  be  necessary  to  accurately 
(jG'termine  the  tail  of  the  flaw  distribution  curve,  where 
low  strength  cracks  cause  failure  to  occur  well  below  the 
mean  strength.  Consequently,  failure  rate  predictions 
cannot  be  made  using  empirical  strength  distribution  curves. 
A method  of  non-destructive  evaluation  of  laser  windows  is 
required  which  will  increase  confidence  rn  reliability 
estimates  v/ithout  impairing  the  serviceability  of  the 
inspected  component.  The  possibility  of  using  the  scatter- 
ing of  light  from  cracks  is  explored  by  dark  field  viewing 
of  pre-cracked  glass  test  specimens. 
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J.2  THEORETICAL  ANALYSIS 


In  order  to  evaluate  a possible  flaw  detection 
technique  the  sizes  of  typical  strength  impairing  cracks 
must  be  found.  The  length  2c  of  cracks  ir  brittle  laser 
window  materials  can  be  obtained  from  the  fracture  toughness 
and  the  nominal  critical  fracture  strength 


K_  = a /ttc 
Ic  c 


(J-1) 


The  corresponding  total  crack  opening  displacement 
2u  of  an  internal  crack  under  a test  stress  c5'  is 

y ^ 

given  for  plane  strain  by  Green  and  Zerna  (1973)  in  terms 
of  the  distance  x from  the  center  of  the  crack,  Poisson's 
ratio  0 , and  the  shear  modulus  G by 

2<51  (1-;)) 


2u  = 

y 


, 2 2,1/2 
(c  -X  ) ' 


(J.2) 
3/2 


For  example,  in  zinc  selenide  with  = 0.75  MN/m' 
(McClintock  and  Hahn,  1975)  a crack  causing  failure  at 
6500  psi  would  have  a half-length  of  90  pm  and  a maximum 
center  opening,  at  a third  of  the  fracture  stress,  of 
0.08  pm.  Equation  J.2  is  plotted  in  Fig.  J.l  showing 
the  maximum  crack  opening  displacemtns  of  cracks  of 
different  strengths  under  different  test  loads  that  would 
be  expected  in  zinc  selenide,  A flaw  detection  technique 
must  be  capable  of  detecting  cracks  of  magnitude  which 
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MAX.  CRACK  OPENING  (;jm) 


Fig . 


J.l  Predicted  maximum  openings  of  cracks  having 

different  strengths  under  different  test  loads 
for  zinc  selenide. 


are  deemed  critical,  typically  with  openings  between  0.02 
yim  and  0.10  yn.  depending  on  the  material,  in  order  to  be 
applied  to  such  high  strength  laser  windows. 


For  scattering  a surface  crack  can  be  considered 
to  be  equivalent  to  a thin  gap  between  two  homogeneous  media. 
Bloembergen  (1975)  showed  that,  for  cracks  with  openings 
small  compared  with  the  wavelength  > of  the  electromagnetic 
radiation  passing  through  the  material  having  refractive 
index  n^  fraction  of  scattered  light  1^  emitted  for  normal 
incidence  is  given  by 


n 


(J.3) 


Such  scattering  of  light  from  cracks  suggests  a possible 
method  of  detecting  the  more  open  low  strength  cracks  and 
even  the  shorter  cracks  in  higher  strength  laser  windows. 


Applying  a small  load,  of  say  10  to  30%  of  the  de- 
sign stress,  will  tend  to  open  any  preexisting  surface 
microcracks,  which  might  be  closed  or  disguised  by 
smearing-over  due  to  polishing.  If  light,  or  whatever 
short  wavelength  electromagnetic  radiation  can  be  trans- 
mitted, is  directed  through  the  edge  of  the  window  running 
parallel  with  and  immediately  below  the  surface  in  tension, 
any  cracks  present  could  be  located  by  the  scattered 
radiation. 
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J.3  EXPERIMENTAL  OBSERVATION  AND  DISCUSSION 

The  light  scattered  from  preorarice':  glass  test 
specimens,  whilst  under  under  carbon-arc  illumination,  was 
viewed  against  a dark  field.  The  precise  crack  opening  was 
subsequently  established  using  a scanning  electron  microscope. 
Cracks  with  openings  of  about  0.1  pm  and  greater  have  been 
observed  by  this  method,  visual  detection  of  any  scattered 
light  was  found  to  be  the  most  effective,  whereas  photographi< 

detection  using  high  speed,  3000  ASA,  film  required  exposure 
tiiiiGs  of  about  120  secs. 


Minimum  observed  crack  openings  are  between  one  and 
five  times  greater  than  maximum  openings  of  typical  median 
strength  impairing  flaws  under  loading  of  10-30%  of  the 
design  stress.  However,  the  depths  of  such  flaws  would  be 
only  about  a tenth  of  the  cracks  viewed  in  glass.  Furthermore, 
some  laser  window  materials,  zinc  selenide  for  example,  are 
only  transparent  to  wavelengths  longer  than  visible.  Con- 
sequently, scattered  intensities  from  cracks  in  such  high 
strength  material  would  be  reduced  by  a factor  of  about  50  for 
calcium  fluoride  to  about  1000  for  zinc  selenide. 

It  IS  concluded  that  until  scattered  intensities 
can  be  Increased  or  detection  methods  Improved,  proof 
testing  must  be  used  to  increase  brittle  structure 
reliability.  However,  a scattered  light  technique  could 
prove  useful  in  deteoting  the  critical,  but  relatively 
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rare,  low  strength  cracks  responsible  for  the  uncertainty  in 

the  tail  of  experimentally  determined  strength  distribution 

curves.  it  is  therefore  recommended  that  further  tests  be 

carried  out  on  the  optically  transparent  brittle  laser  window 

materials,  to  establish  the  extent  to  which  such  cracks  can  be 
detected . 
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K.  SUMMARY 


The  space  charge  theory  for  ionic  materials  has  been 
modified  to  include  a secondary  interaction  force  between 
an  aliovalent  ion  and  a grain  boundary  or  surface  which  acts 
in  addition  to  the  usual  electrostatic  interaction  forces. 

In  a purely  ionic  system  such  a force  would  likely  arise 
from  an  ion  size  misfit,  i.e.,  a strain  energy  which  is 
relieved  near  the  boundary.  The  space  charge  problem  has 
been  reformulated  to  include  such  a force.  Numerical 
solutions  show  that  the  electrostatic  field  and  the  dis- 
tribution of  solute  ions  and  defects  near  the  boundary  are 
substantially  modified  by  the  coupling  of  the  two  inter- 
actions through  the  required  satisfaction  of  Poisson's 
equation  by  all  charged  species.  Solutions  were  also 
obtained  for  the  case  of  two  solute  ions  with  different 
strain  interactioi  forces.  The  distributions  of  the  two 
ions  are  found  to  be  interdependent. 

Boundary  mobility  calculations  were  performed  taking  the 
secondary  interaction  and  the  modifying  effect  on  the  space 
charge  cloud  into  account.  Because  of  uncertainty  about 
the  form  of  the  strain  field  and  the  velocity  dependence  of 
the  electrostatic  field  the  resulting  calculations  have  more 
uncertainties  than  do  the  calculations  for  the  purely 
electrostatic  case.  Nevertheless  some  interesting  results 
are  indicated.  The  mobility  is  affected  most  strongly  at 
temperatures  near  the  isoelectric  temperature.  Here  the 
strain  interactions  induce  much  stronger  electrostatic  forces 
than  would  otherwise  exist.  These  are  apparently  helpful  in 
preventing  breakaway  of  the  boundary  from  the  solute  could. 
The  theory  predicts  that  the  effectiveness  of  divalent 
solutes  in  inhibiting  boundary  motion  should  increase  in 
the  order  of  Ba^^,  Sr^^,  to  Ca''^^. 

In  many  cases  it  is  physically  possible  for  a boundary 
to  breakaway  from  its  solute  cloud  and  move  at  the  much 
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higher  intrinsic  mobility.  At  intermediate  driving  forces 
this  apparently  happens  by  a spatial  fluctuation  of  the 
boundary  which  induces  a nonuniform  distribution  of  solute 
along  the  boundary.  The  more  rapidly  moving  parts  of  the 
boundary  accelerate  and  shed  the  solute  cloud.  This  break- 
away region  can  spread  laterally  and  carry  the  entire 
boundary  into  the  high  velocity  mode.  Such  rapidly  moving 
boundaries  can  then  cause  much  faster  recrystallization  or 
discontinuous  grain  growth.  Thus,  breakaway  of  a boundary 
from  the  solute  cloud  can  be  an  important  step  in  controlling 
recrystallization  in  forged  materials  as  well  as  in  prevent- 
ing discontinuous  grain  growth  and  pore  entrapment  during 
sintering  of  oxides.  A nonlinear  perturbation  analysis  was 

used  to  indicate  the  conditions  under  which  such  breakaway 
can  occur. 


Foi  typical  oxides  four  possible  mechanisms  may  control 
the  grain  boundary  mobility.  These  are  the  intrinsic 
mobility,  solute  drag,  pore  or  particle  drag,  and  liquid 
phase  control.  These  have  been  evaluated  to  consider  the 
relative  importance  and  possible  mechanisms  for  transitions 
among  the  regimes  suggested.  The  particular  mechanism 
controlling  the  mobility  can  affect  the  kinetics  of  grain 
growth  and  the  resulting  grain  morphologies.  A survey  of 

the  data  for  many  ceramics  indicates  qualitative  agreement 
with  theory. 

Single  crysta  l a a^f  CaF2  were  hot  forged  at  temperatures 

between  0.46  (47b'>C)  and  0 . 9 j (1230  = C)  of  the  melting 

point.  The  flow  stress  temperature  curve  was  found  to  have 

a mild  orientation  dependence  for  T < 0.65  Ttt,  with  a > 

_ _ m wirn  > 

<110>  '^<111> 


The  oubcfcil  sii..;,*  jjas  found  to  be  inversely 


proportional  to  the  forging  t fr-eiss.  .^\t  T > 0.65  T^ 
recrystallization  occurred  during  forging.  The  nature  of 
the  microstructure  developed  depends  upon  the  crystal 
orientation.  In  general  <100>  crystals  have  more  uniform. 
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equiaxed  subgrains,  compared  to  <111>  crystals  which  tended 
to  be  more  nonuniform  with  elongated  subgrains.  At  lower 
temperatures  the  microstructures  were  less  uniform, 
especially  in  the  <111>  crystals  in  which  deformation  bands 
were  common.  Crystals  forged  in  the  <100>  direction  showed 
less  tendency  to  crack  during  forging  or  cooling  than  <111> 
crystals , 

At  lower  temperatures  the  subgrain  boundaries  are  less 
Well  knit  and  contain  considerable  excess  or  extrinsic 
dislocations.  At  higher  forging  temperatures  the  subgrain 
boundaries  become  equilibrated  and  there  are  fewer  if  any 
extrinsic  dislocations  found  within  them.  The  specific 
dislocation  reactions  involved  in  sub-boundary  equilibration 
were  identified.  The  misorientation  angles  of  the  sub- 
boundaries are  of  the  order  of  a degree  although  there  is  an 
indication  of  a bimodal  distribution  in  which  lower  angle 
subceMs  form  within  higher  angle  subgrains. 

Models  were  developed  for  calculating  the  effects  of 
microcracking  and  microcrack  coalescence  on  fracture 
behavior  and  statistics.  The  problem  is  to  determine  the 
maximum  stress  at  which  catastrophic  crack  propogation  occurs 
in  a body  containing  a distribution  of  subcritical  flaws. 

From  the  models  the  stresses  for  first  cracking  and  for 
failure  were  determined. 

The  effects  of  principal  stress  ratio  on  brittle  crack 
statistics  were  studied  for  homogeneous  arrays  of  regular 
hexagons,  having  grain  boundary  strengths  correspono.ing  Lo 
an  extreme- value  distribution.  In  a numerical  nodej.-  ;ie 
cracks  were  simulated  by  dislocation  pairs,  choker.  >o  satisfy 
the  stress-free  conditions  on  each  cracked  grain  segments. 
Cracking  was  assumed  to  occur  in  the  grain  boundai.y  t,h  fhe 
lowest  strength  relative  to  the  local  stress;  or  stress 
intensity  factor  on  it  due  to  the  applit ' strera  acting 
through  the  current  crack  pattern.  An  analyt.c  apt  oximal , 
was  developed,  based  on  the  density  of  grains  croi;i«;ed  so  far 
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and  the  probability  that  they  would  aggregate  into  a crack 
of  critical  length  under  the  current  applied  stress. 

For  a typical  strength  distribution,  the  stress  at 
first  cracking  rises  by  50%  in  going  from  biaxial  to  uniaxial 
applied  stress.  The  maximum  stress  shows  less  size  effect 
than  predicted  from  extreme-value  statistics,  and  a some- 
what smaller  rise  in  strength  in  going  from  biaxial  to 
uniaxial  stress.  Crack  closure  and  friction  begin  at  a 
compressive  stress  of  about  half  the  applied  tensile  stress. 

The  problem  of  calculating  the  progress  of  micro- 
cracking  in  brittle  materials  towards  the  configuration 
giving  maximum  load,  by  repeated  stress  analysis  of  the 
developing  micro-crack  distribution,  had  been  treated  initially 
by  considering  grain  boundary  cracking  in  a two-dimensional 
array  of  regular  hexagon  grains  embedded  within  an  infinite 
body.  Since  brittle  materials  often  fracture  by  the 
propagation  of  surface  cracks,  the  numerical  model  was 
extended,  using  near  surface  dislocation  dipoles,  to  study 
the  effects  on  brittle  crack  statistics  of  surface  or  near 
surface  initiated  fracture. 

Part  sizes  up  to  1000  grain  segments  were  studied 
using  the  numerical  model  for  two  different  variabilities. 

It  was  found  that  the  maximum  strengths  for  surface 
initiated  fracture  are  10-20%,  depending  on  orientation  and 
variability,  below  those  for  embedded  fracture.  The 
relatively  small  drop  in  maximum  strengths  suggest  that  the 
statistical  treatment  of  brittle  materials  is  somewhat 
independent  of  whether  fracture  is  initiated  in  the  bulk  or 
at  the  surface  of  the  material  as  long  as  the  relevant  dis- 
tribution of  grain  boundary  strengths  is  used. 

The  method  of  determining  the  flaw  density  distributions 
of  brittle  materials  from  hardness  indentation  was  summarized. 
The  method,  developed  earlier  in  this  project,  neglected  the 
interfacial  friction  between  the  indenter  and  the  specimen 
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surface.  Thus,  the  apparent  strength  calculated  was  too 
high.  It  was  suggested  that  the  effect  of  friction  could  be 
included  in  the  analysis  by  modifying  the  stress  distribution 
function . 

A plane  strain  cylindrical  indenter  was  proposed  for 
studying  sub-surface  flaws.  Such  an  indenter  gives  no  stress 
on  the  specimen  surface  and  onl.y  compressive  stress  in  the 
interior.  Friction  again  influences  the  location  of  maximum 
stress . 

The  possibility  of  determining  strength  distributions 
from  acoustic  emission  was  studied.  Acoustic  emissions  were 
used  to  obtain  the  statistical  parameters  that  govern  crack 
initiation  in  grain  boundaries,  idealized  as  an  extreme 
value  distribution  of  the  third  asymptotic  kind.  It  is  found 
that  polishing  significantly  reduces  the  amount  of  acoustic 
activity  as  well  as  increasing  the  value  of  the  statistical 
parameter  m. 

A flaw  detection  technique  must  be  capable  of  detecting 
cracks  of  magnitude  which  are  deemed  critical,  typically 
with  openings  between  0.02pm  and  0.10pm  depending  on  the 
material.  The  possibility  of  using  the  scattering  of  light 
from  cracks  was  explored  by  dark  field  viewing  of  pre-cracked 
glass  test  specimens.  Cracks  with  openings  of  about  0.10pm 
and  greater  have  been  observed.  Therefore,  until  scattered 
intensities  can  be  increased  o."  detection  methods  improved, 
proof  testing  must  be  used  to  increased  brittle  structure 
reliability. 
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